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PREFACE 


This book follows closely a lecture course I gave entitled ‘Astronomical 
Spectroscopy' to third-year undergraduate students at University College 
London between 1998 and 2003. The students who attended had done 
a prior introductory course on Quantum Mechanics which covered the 
hydrogen atom but no further atomic physics or spectroscopy. A similar 
level of prior knowledge is assumed in the current work. 

There are many people whose help have been essential for the comple- 
tion of this book. First I must thank Bill Somerville who inaugurated the 
course Astronomical Spectroscopy and taught it for two years before me. 
He selflessly shared his lecture notes and other materials with me. I would 
like to thank Ceinwen Sanderson for turning my hand-scrawled lecture 
notes into ATEX, and my colleagues Tony Lynas-Gray, Bill Somerville, Peter 
Storey and Jeremy Yates for their extensive comments on the draft of the 
book. I owe a debt of gratitude to my graduate students Bob Barber and 
Natasha Doss who checked all the problems and found many errors. I 
thank all of them for the corrections; any errors that remain are all mine. 

I must also thank the students who attended my Astronomical Spec- 
troscopy course. It was great fun to teach, not least because the latest 
developments in astrophysics often fed straight into the lectures. Partic- 
ular thanks are due to the class of 2003 who made a number of helpful 
comments and suggestions on the contents of the book. 

A book on spectroscopy thrives on good illustrations and I have 
shamelessly plundered the literature and other sources for spectra to illus- 
trate this one. I must thank Xiaowei Liu for help with digitising many 
of the published spectra, my student Iryna Rozum, my son Matthew, 
and especially David Rage for their help with the other illustrations. 
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I thank the journal publishers and the many authors who greeted my 
requests to reproduce their work with prompt enthusiasm, especially 
those authors who adapted figures at my request. Each journal and author 
is individually acknowledged in the figure captions. 

Finally I must acknowledge the UCL astronomers of the past and 
present who have answered my many questions on astrophysics with a 
patience their frequent stupidity probably did not deserve. Particularly 
high on this list are Pete Storey and Mike Barlow, but the rest of the varied 
lunch crew should not be forgotten. Without you my knowledge of things 
astronomical would be the same as it was the day I arrived at UCL — 
nothing. 


Jonathan Tennyson 
London, 
July 2004 
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CHAPTER ONE 


WHY RECORD SPECTRA OF 
ASTRONOMICAL OBJECTS? 


“We will never know how to study by any means the chemical 
composition (of stars), or their mineralogical structure 


— Auguste Comte (1835) 


1.1 A Historical Introduction 


In the first part of the 19th century, astronomers began to make paral- 
lax measurements which revealed for the first time how distant even the 
closest stars are from us. Since travel to the stars was, and still is, impossi- 
ble with foreseeable technology, many scientists believed that the compo- 
sition and character of the stars would forever remain a mystery. This view 
is pithily summarised by the quote from the positivist French philosopher 
Auguste Comte (1798-1857) given above. 

Today, the composition of stars, and indeed of the diffuse material in 
the large spaces in between the stars, is well known. How did this situation 
come about? In fact the first steps to finding the solution to the problem 
had been taken even before Comte began writing. 

In 1814, Joseph von Fraunhofer (1787-1826) used one of the high- 
quality prisms he had manufactured to diffract a beam of sunlight, taken 
from a slit in his shutters, onto a whitewashed wall. Besides the character- 
istic colours of the rainbow, which had been observed in this fashion since 
Newton, he saw many dark lines (see Fig. 1.1). He meticulously catalogued 
the exact wavelength of each dark line — which are still known today as 
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Fig. 1.1. The solar spectrum as recorded by Fraunhofer. 


Fraunhofer lines — and labelled the strongest of them with letters. Many 
of these labels, such as the sodium D lines (see Sec. 6.4) are still used 
today. Fraunhofer not only recorded the first astronomical spectrum, he 
recorded the first-ever high-resolution spectrum. Fraunhofer's spectrum 
was the first to resolve discrete line transitions. 

Fraunhofer did not know what caused the dark lines he observed. 
However he performed a similar experiment using light from the nearby 
red-star Betelgeuse and found that the pattern of dark lines he observed 
changed significantly. Fraunhofer concluded correctly that most of those 
features were somehow related to the composition of the object he was 
observing. In fact some of the lines were due to the Earth's atmosphere, 
the so-called telluric lines. For example, the features Fraunhofer marked 
A and B in his solar spectrum are actually due to molecular oxygen in our 
own atmosphere. 

The first real step in understanding Fraunhofer's observations came in 
the middle of the 19th century with the experiments of Gustav Kirchhoff 
(1824-1887) and Robert Bunsen (1811-1899). These scientists studied the 
colour of the light emitted when metals were burnt in flames. They found 
that in certain cases the wavelength of the emitted light gave an exact 
match with the Fraunhofer lines. The sodium D lines, which give sodium 
street lights their characteristic orange colour, were one such example. 
These experiments demonstrated that the Fraunhofer lines were a direct 
consequence of the atomic composition of the Sun. 

Any understanding of how these lines came about had to wait 
until the arrival of the 20th century with the revolution of scientific the- 
ory represented by quantum mechanics. The developments of quantum 
mechanics and spectroscopy have always been closely linked. As it is 
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through the study of spectra that we have learnt of many of the riches 
in the Universe around us, the development of astrophysics has also been 
closely linked to that of spectroscopy and quantum mechanics. This book 
aims to give an introduction to the spectroscopy of atoms and molecules 
that are important for astrophysics. This book is not a text on quantum 
mechanics, and indeed, some basic knowledge of quantum mechanics is 
assumed, for it is not possible to understand or interpret spectra without 
some understanding of quantum mechanics. 

Hearnshaw (1986) gives a fascinating historical view of the relation- 
ship between astronomy, spectroscopy and the technical developments in 
both fields (see further reading). 


12 What One Can Learn from Studying Spectra 


Essentially all information about astronomical objects outside the solar 
system comes through the study of electromagnetic radiation (light) as 
it reaches us. This light can contain much detailed information which 
is only obtained by careful analysis. Generally speaking, one can clas- 
sify the information obtained by observing light according to the spec- 
tral resolution; that is the degree of sensitivity to different wavelengths, 
used to make the observation. One can classify such observations using 
the following general categories. 

When one looks at the night sky with the naked eye, most astronom- 
ical bodies appear white. White light is actually light that is composed 
of many wavelengths which are not resolved into their different colours. 
Monitoring white light gives the positions of objects in the night sky. It can 
be used to construct maps of stars and galaxies. It can also be used to plot 
the movements of heavenly bodies such as comets through the night sky. 

If one looks carefully at some celestial objects, such as the planets Mars 
and Jupiter, or stars such as Betelgeuse, one can see that these objects are 
tinged with a certain colour. Using instruments with low resolving power, 
it is possible to separate the light arriving at Earth into broad band colours. 
Observing colours tells us something about temperatures. For example, 
blue stars are hotter than red ones; objects that emit X-rays, such as the 
solar corona, are very hot, whereas cold objects may only emit light of 
very long wavelengths such as radio waves. 

The most detailed astrophysical information is only obtained from 
high-resolution studies which involve detecting the light arriving at the 
earth as a function of its component wavelengths. This allows detailed 
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spectroscopic features to be identified separately from broad band features 
such as colour. At the highest resolution, such studies not only yield the 
central wavelength of any feature, often referred to as a line, but also the 
shape of the feature. Such studies can yield significant extra information 
and this book is largely devoted to the physical basis of this information 
and how it can be interpreted. 

To interpret an astronomical spectrum, one needs considerable knowl- 
edge of atomic and molecular physics. This knowledge usually comes 
from laboratory studies which provide the basic physical parameters 
necessary for understanding the astronomical spectrum. There is a direct 
relationship between these physical parameters and the astronomical 
information that can be obtained by observing spectra. Thus for any line 
Observed in an astronomical spectrum, one can potentially use laboratory 
data to extract the following information. 

The composition of the object being observed can be inferred 
by knowing which atom (or ion or molecule) produces the observed 
transition. 

The temperature and other physical conditions can be deduced from 
assigning the actual transition being observed to precise energy levels in 
the atom. Transitions take place between many different states in a partic- 
ular atom. Knowing which states are involved gives direct information on 
the degree of excitation of the system. This can be used to determine the 
physical conditions, such as the temperature or density of the environment 
local to the system. 

The abundance of the species undergoing the transition can only be 
determined if the intrinsic strength of the transition being observed is 
known. Line strengths can be hard to determine in the laboratory. Astro- 
nomically, the strength of a transition is directly related to the number of 
atoms undergoing the transition under suitable conditions of optical depth 
(see below). Knowledge of the intensity of transitions is therefore impor- 
tant for determining the abundance of any species. 

Motions of the species being observed relative to the earth, or indeed 
the whole region containing the species, lead to a shift in the wavelength 
of the line; this shift is known as the Doppler shift. The Doppler shift is the 
change in the line position from the position measured in the laboratory. 
This shift is given by the Doppler formula, 


vo AX 


Why Record Spectra of Astronomical Objects? 5 


where v is the velocity of the source in a direction away from us, 
c = 2.99792458 x 105 m - s^! is the speed of light, A is the rest wavelength 
of the transition and A) is the change in wavelength, known as the Doppler 
shift. Application of this formula requires laboratory measurement of the 
rest wavelength to high accuracy. Formula (1.1) is for non-relativistic 
motions. When an object is moving towards us, the transition is shifted 
to shorter wavelengths (‘blue-shifted’), and when the object is moving 
away from us, it is shifted to longer wavelengths (‘red-shifted’). It was 
through the monitoring of Doppler shifts of spectra of hydrogen atoms 
that allowed Edwin Hubble (1889-1953) to show in 1929 that our universe 
is uniformly expanding and so started from a single point or Big Bang. 

The pressure or density of the environment local to the species under- 
going the transitions can be monitored by observing the line profile. 
Such observations require particularly high resolutions. Spectral lines are 
broadened by collisions between species; the more frequent these col- 
lisions are, the greater the broadening. This process is called 'pressure 
broadening’. Lines are also broadened by the thermal motions according to 
the Doppler formula. Doppler broadening arises because hot species move 
about faster than cold ones. Both of these reveal information about the 
physical environment of the species being observed. However, the com- 
bined effects of pressure and temperature on the line profile can only be 
resolved using ultrahigh resolution observations. 

Any magnetic field present can be monitored as certain spectral lines 
will be split into more than one component. Energy levels of states which 
possess angular momentum are split in the presence of a magnetic field. 
The result is that a single transition can become two or more distinct tran- 
sitions. The degree of separation between these component lines depends 
directly on the strength of the local magnetic field. Such splittings, if 
Observed, can therefore provide a measurement of this field. 

The information obtained from such observations is the key to most 
astronomical knowledge. However, to interpret any astronomical spec- 
tra requires detailed information about the intrinsic properties of atomic 
spectra. For each atom or ion or molecule being observed, one needs to 
know: 


(1) Its important spectral lines: these are often summarised using figures 
called Grotrian diagrams (see Sec. 5.4). 
(2) Its energy level structure: also summarised on Grotrian diagrams. 
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(3) The intrinsic line strength of the transition(s) being observed. 
(4) The precise rest (ie. laboratory) wavelength of any transition 
observed. 


Additional information is required to interpret pressure broadening of 
spectral lines and splitting in magnetic fields, however these topics will 
not be pursued in this book. Understanding and use of all this detailed 
spectroscopic information requires considerable knowledge of quantum 
mechanics. 

At all wavelengths there are observed spectral lines which have yet 
to be identified (see for example Figs. 6.8 and 7.6). A particularly long- 
running current example are the diffuse interstellar bands or DIBs. This 
means that laboratory astrophysics, the study of astrophysical processes 
in the laboratory, based on either experiment or theory (or both), remains 
an active area of research. 


Problems 
Answers to problems are given at the end of the book. 


1.1 While observing stars in a distant galaxy, Edwin Hubble observed dis- 
crete line emissions at 411.54 nm, 435.50 nm, 487.75 nm and 658.47 nm. 
There are H-atom transitions with rest wavelengths of 410.17 nm, 
434.05 nm, 486.13nm and 656.28nm. Verify that these lines are all 
Doppler-shifted by the same amount. What is the speed of the distant 
star relative to earth? Is it moving towards us or away from us? 


CHAPTER TWO 


THE NATURE OF SPECTRA 


“We all know what light is; 


but it is not easy to tell what it is.’ 


— Samuel Johnson, quoted by Boswell (1776) 


2.1 Transitions 


All atoms have a series of discrete, quantised energy levels. In fact there 
are infinitely many of them. A spectral line can be obtained when a 
jump between two of these levels with different energies occurs with 
light of the correct wavelength. This wavelength, A, corresponds to the 
exact energy difference, E, between the energy levels via the Planck 
relationship 
he 

E=hv= A! (2.1) 
where v is the frequency of the light and c is the speed of light. h is Planck's 
constant, 6.626068 x 10 ?* J.s, and hv gives the energy carried by each 
particle of light, known as a photon. 

Each atom absorbs light at a series of characteristic wavelengths. 
While individual transitions belonging to different species may acciden- 
tally coincide, the whole series of lines is unique for each atom, and indeed, 
every ionisation stage of each atom. It is often said that the spectrum of an 
atom gives a unique fingerprint for that atom which distinguishes it from 
all others. A closer analogy is that of the barcode used by supermarkets 
and elsewhere to uniquely tag their products. The spectrum of an atom is 
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a barcode which, if we can read it, yields much detailed information about 
the atom. 


2.2 Absorption and Emission 


An atom can either absorb light, jumping to a higher-lying energy level, or 
emit light, dropping to a lower energy level (see Fig. 2.1). Not all tran- 
sitions are equally likely; each transition is therefore characterised not 
only by a precise wavelength but also by a probability of that transition 
occurring. 

Any level in an atom can absorb light, but it requires light of the 
correct wavelength from another source to make the atom jump to a 
higher energy level and for the photon to be absorbed. This is rela- 
tively easy to arrange in the laboratory but will not occur in all astro- 
nomical environments. A typical situation where an absorption spectra 
is observed comes from the atmospheres of stars, where the core of 
the star provides a continuum light source (see Fig. 2.2). This light 
source is approximately a black body curve with the temperature of the 
star. Species in the photosphere of the star are observed in absorption 
against this curve. This is the nature of Fraunhofer's solar spectrum (see 
Fig. 1.1). Absorption in the interstellar medium against a more distant star 
can also be observed given a suitable arrangement of the astronomical 
bodies. 


Absorption Emission 
Level 3 
Er N. > 
Level 2 
hv hv 
DIN PIN 
Level 1 








Fig. 2.1. Emission and absorption in a schematic three-level system. 
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7 
Pd Photosphere 
(cooler) 
Hot stellar — 
Interior — 
~a 
— 
N 
N 
N 


Fig. 2.2. The spectrum of a star: spectral lines appears as absorption against a 
continuous black body curve. 


Emission requires that the atom involved starts in an excited state. 
Under such circumstances the atom can spontaneously emit a photon. The 
intensity of an emission spectrum from level j to level i is given by 

NjhvjA ji , (2.2) 
where N; is the number of atoms in state j per unit volume, and hv; gives 
the energy difference between the levels. A j; is called the Einstein A coef- 
ficient for spontaneous emission and it gives the number of transitions 
per second of an atom from excited state j to state i. This is an impor- 
tant quantity for measuring the likelihood of a particular transition. The 
Einstein A coefficient, which is usually given in units of s !, will be used 
throughout this book to establish both transition strengths and timescales 
for transitions. 

The intensity with which an atom in level i absorbs light to jump to 
level j is given by 


NihvjiBijp, , (2.3) 


where N; is the number of atoms in state i, and p, is the density of radiation 
with frequency v. Bj; is the transition probability for absorption, called the 
Einstein B coefficient. 

Einstein proved that 





(2.4) 
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where g; and g; are the statistical weights, also known as degeneracies, of 
states j and i respectively. This proof can be found in standard textbooks 
such as Bransden and Joachain (2003) (see further reading). 


2.3 Other Measures of Transition Probabilities 


The Einstein Aj; and Bj; coefficients can be calculated using the wave- 
functions of states i and j. Although Einstein A coefficients will be used 
throughout this book to quantify transition probabilities, there are a num- 
ber of other ways of doing this. 

One often-used quantity is the oscillator strength: 
p 4T EoMe 
where me is the mass of an electron and e is the charge on an electron. 
The factor 47 €o, where c, is the permittivity of a vacuum, enters when the 
expression is given in SI units. 

If g; is the degeneracy of level i, it can be shown that 


&ifij ^ —gifj (2.6) 


where the convention that f is positive for absorption and negative for 
emission has been used. f has the advantage of being dimensionless, so 
it can be tabulated unambiguously. Furthermore, oscillator strengths also 
have interesting sum rule properties which can be used, for example, to 
quantify missing absorption. For more information see Woodgate (1983) in 
further reading. Laboratory data for atomic systems are usually presented 
as f-values. 





pm hvjiBij, (2.5) 


2.44 Stimulated Emission 


Figure 2.3 shows a third way in which an atom can undergo a transition. 
This is a stimulated rather than a spontaneous emission. To get a stimu- 
lated emission, a photon of the correct wavelength is required to initiate 
the emission of the photon. The two outgoing photons are coherent and 
travel in the same direction. If there are a significant number of atoms in 
a particular excited state, then stimulated emission can provoke a cascade 
of photons which can result in laser (Light Amplification by Stimulated 
Emission of Radiation) action. 
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Level 3 Level 3 


hv 
LUN 
hv UA ay 
AUAN 
NAVES hv hv 
Level 2 US Level 2 ATTN 
Level 1 Level 1 


Fig. 2.3. Stimulated emission leading to laser action in a schematic three-level 
system. 


To achieve laser action it is necessary to not only get stimulated emis- 
sion but to also get population inversion: a greater population in the 
upper state than the lower state. The original laboratory lasers, reported 
by Townes and Schawlov in 1954, worked in the microwave. Maser 
(Microwave Amplification by Stimulated Emission of Radiation) action is 
also important astronomically. Specific astronomical examples will be dis- 
cussed in Secs. 3.10 and 10.5. 


2.5 Optical Depth 


If many, many atoms of a certain species lie in the line of sight of a par- 
ticular observation, then all the atoms may not be required to absorb all 
the light of the wavelength under study. To quantify this situation it is 
usual to define a quantity called the optical depth, T, which is a measure 
of transparency of the medium. The optical depth is defined by 


r= ns Q.7) 


where the integral runs over the path of the light being absorbed. The 
extinction coefficient, K, represents the product of the number density 
of the atoms and their opacity at the wavelength in question, where the 
opacity is a measure of a material's ability to absorb or block photons. 

If the optical depth is large (T >> 1), then a transition is classed as opti- 
cally thick. Under these circumstances, it will not require all the atoms to 
absorb all the light at that wavelength. In this case the total absorption 
Observed is not directly related to the number of absorbers, usually called 
the column density. Conversely, in the limit of only a few photons being 
absorbed, the optical depth is small (7 « 1). In this case the line is optically 
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thin and the intensity of the observed spectrum is directly proportional to 
the number of absorbers. Under these circumstances, the spectrum can be 
used to measure the column density directly. 

To obtain column densities and other related information from spec- 
tra, it is important to know the optical depth of the observed transitions. 
Sometimes it is readily apparent that a particular transition is saturated, or 
optically thick (see the central feature in Fig. 3.9 for an example). However, 
it is often not obvious for a given observation whether a particular line is 
optically thin or not. In some cases a knowledge of atomic physics can be 
used to resolve this problem. An important example of this is discussed 
in Sec. 6.4. 

The transport of radiation through astronomical objects is a compli- 
cated issue. It is dealt with more thoroughly by Emerson (1996) (see fur- 
ther reading). 


2.6 Critical Density 


From a spectroscopic point of view, whether a particular emission behaves 
as if the emitting species is in thermodynamic equilibrium with its local 
environment depends on something called the critical density. For den- 
sities above the critical density, emissions are thermal in origin, whereas 
below this density, each collision leads to an emission. The critical den- 
sity n, thus depends on the ratio of the timescale for decay by emission 
to collisional de-excitation. Specifically, the critical density for level j is 
given by 


—OXEiejAji 


Yizidji 


C 


where Ajj is the Einstein A coefficient and the sum runs over all pos- 
sible emissions, and qj; is the rate for collisional de-population of level 
j, summed over all possible processes. This expression often simplifies 
to the ratio of two numbers, since for many systems and environments 
there is a single important route for emission and a dominant collisional 
de-excitation process. 

The critical density is the density of collision partners, often elec- 
trons, above which the collisional de-excitation from the upper level 
occurs more quickly than the radiative de-excitation. In this case the 
emissions behave as if they are thermal and can be used as a measure of 


The Nature of Spectra 13 


temperature. Below the critical density, each collisional excitation leads to 
an emission and, given knowledge of both the line strength and the colli- 
sion cross section, the emissions can provide information on the density of 
the collision partners. 

Much of the matter in the Universe is not in local thermodynamic 
equilibrium (LTE), as determined by the critical density. In non-LTE condi- 
tions, the concept of temperature is of limited use, although the free elec- 
trons are still usually thermalised to some temperature; such media are 
often characterised by the electron temperature Ty. 

This definition of critical density is distinct from the cosmologi- 
cal critical density which determines whether the Universe is open or 
closed. 


24 Wavelength or Frequency? 


Measurements of spectra obtained in the laboratory are usually reported 
in frequency units such as Hz (which are s !) or the related wavenum- 
ber unit of cm !. Frequencies and wavenumbers are directly propor- 
tional to the energy jump of the transition. They are therefore much 
more amenable to physical understanding and interpretation. For this 
reason, most of the discussion presented below will be given in the 
energy / frequency / wavenumber domain. 

Conversely, astronomers, at least in the infrared, visible and ultravi- 
olet, tend to use wavelengths, typically um, nm or A (equal to 1076, 107? 
and 10^? m, respectively). Spectrographs on telescopes work naturally in 
wavelengths. Indeed, the spectral resolution of an instrument is given by 
its resolving power R, as 


R=5, (2.9) 


where AA is the smallest wavelength difference that can be resolved. The 
ratio 5 gives the velocity resolution. Thus a resolving power of R = 30000 
can resolve velocities in excess 10 km - s71. 

The dichotomy between working in the wavelength or frequency 
domain means frequent application of the standard formula 


\= (2.10) 


C 
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When working in wavenumbers I find it helpful to remember that 


10000 
A(um) = 


TE Q.11) 


At some wavelengths, particularly at radio frequencies, it is stan- 
dard practice for astronomers to present spectra as Doppler shifts from 
the (known) rest (that is laboratory) wavelength Ao. This works via the 
Doppler formula: 


U,—C—, (2.12) 


where A, represents a shift in wavelength. c is the speed of light and v, is 
the velocity of the source away from us; both usually given in km s !. This 
method provides a convenient way of showing the velocity structure of the 
object(s) under study. Indeed, using velocity shifts means that several tran- 
sitions from the same object can be shown on the same frequency scale; 
Figures 3.17, 3.22 and 10.19 are three of several examples given in this 
book. However, one note of caution must be exercised. This representa- 
tion assumes that all transitions in a particular observational window have 
the same rest frequency Ao. This means that transitions that are nearby in 
frequency but due to physically distinct processes will appear in the spec- 
trum as velocity shifts. Examples of this are quite common: see for instance 
the top panel in Fig. 3.17 where the transition labelled H509 is not given 
relative to its rest wavelength. 


28 The Electromagnetic Spectrum 


Figure 2.4 gives an overview of the electromagnetic spectrum. The most 
tightly defined and smallest region in the electromagnetic spectrum is the 
visible region, which covers wavelengths between 4000 À and 7000 Á. For 
astronomical studies, it is important to have an appreciation of the differ- 
ent regions of the electromagnetic spectrum for a number of reasons. 

At the theoretical level, different physical processes occur at different 
energies, which can then be directly associated with particular regions of 
the spectrum. Astronomically there is also a correlation between wave- 
length and temperature, particularly for emission spectra. Thus, for exam- 
ple, very hot regions are bright at short, X-ray wavelengths, whereas 
cold interstellar clouds are extensively studied at long, microwave wave- 
lengths, often referred to as radio waves. 
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Fig. 2.4. The electromagnetic spectrum. 


From an astronomical perspective there are also practical reasons for 
understanding the electromagnetic spectrum. All ground-based observa- 
tions are necessarily made through the earth's atmosphere. The atmosphere 
is largely transparent at visible and radio wavelengths. Observations are 
possible through several infrared ^windows' in the atmosphere. These win- 
dows lie at wavelengths between regions where water vapour, and other 
atmospheric species absorb much of the incoming light. However in all 
wavelength regions there are atmospheric species which cause isolated 
absorptions known as 'telluric' lines. One can seek to reduce the effects 
of the atmosphere by the choice of good observing sites, such as the high, 
dry, near-equatorial sites in Hawaii and Chile, which are both home to 
clusters of telescopes. However there are many wavelength regions where 
ground-based observations are simply impossible (see Fig. 2.5). 
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Fig. 2.5. Absorption of light of different wavelengths in the Earth's atmosphere. 
Gray areas show strong telluric absorption. 


The development of satellite technology over the last quarter-century 
means that there have now been space missions flown which study all 
wavelengths. Observing from space is very expensive and observing 
time is therefore necessarily limited. Furthermore, state-of-the-art satellite 
observations usually yield much poorer spectral and spatial resolutions 
than is obtainable on the ground because it is not practical to fly such large 
telescopes or such good spectrometers. Large telescopes and long integra- 
tion times are also needed to observe faint objects. 


Problem 


2.1 The top panel of Fig. 3.17 in Chapter 3 shows two emission features 
labelled H40a and H509. The velocity scale is relative to the rest fre- 
quency of H40a, which is at 99.023 GHz. Estimate the rest frequency 
of the H508 transition. 


CHAPTER THREE 


ATOMIC HYDROGEN 


‘The two most common things in the Universe are Hydrogen and Stupidity.’ 


— Harlan Ellison (1934) 


3.1 Overview 


Spectra of atomic hydrogen, H, are of paramount astronomical impor- 
tance. This is because approximately 90% of atomic matter by number is 
hydrogen. This occurs in a variety of forms: H* or protons, H atoms, H» 
molecules and indeed the molecular ions H7 and Hj. The spectral lines of 
hydrogen are prominent in a great variety of astronomical objects and are 
much studied. All aspects of hydrogen spectroscopy therefore need to be 
considered in detail. 

The spectrum of atomic hydrogen also plays an important role in the 
theory of quantum mechanics. Hydrogen is the simplest atom, comprising 
a single electron and a proton. It is the only atom for which exact quantum 
mechanical solutions can be found for its energy levels and wavefunctions. 
These solutions will be used extensively here but will not be derived. Such 
derivations are a standard part of most introductory courses in quantum 
physics [see Rae (2002) in further reading]. 


3.2 TheSchródinger Equation of 
Hydrogen-Like Atoms 


Any atom comprising a single electron orbiting a nucleus of charge Z can 
be described as hydrogen-like. The Hamiltonian operator for this system 
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can be written 


H= Hg _ Ze ] 
2u 4TEor 





(8.1) 


where r is a vector, of length r, linking the electron to the nucleus. This 
Hamiltonian comprises two terms. The first term, given by the Laplacian 
operator, is the kinetic energy operator for the electron. The second term 
is the potential energy term, in this case the Coulomb attraction between 
the electron, charge —e, and the nucleus, charge + Ze. 

The Hamiltonian (3.1) contains a number of constants. It is often given 
in terms of atomic units, which scale the dimensions of the problem to the 
atomic scale. In atomic units (see Sec. 3.4 for details), the charge on an 
electron e, the mass of an electron m,, h = x and 47r€9 all equal one. This 
simplifies the Hamiltonian to: 

H= age: es (8.2) 
2u r 

In general the Schrödinger equation for a system of energy E and with 

wavefunction w is written 


fy = Ey. (3.3) 


For the particular case of the hydrogen-like atom in atomic units, the 
explicit form of the Schrédinger equation is therefore 


1 Z 
a -Z= E| yA =0. (3.4) 


3.3 Reduced Mass 


The hydrogen atom is actually a two-particle problem, concerning the 
motion of an electron of mass m, and a nucleus of mass M. To solve for the 
internal motion of this system it is necessary to rewrite the equations-of- 
motion of the system into one equation representing the overall translation 
of the whole system, with mass M + me, in space, and a second equa- 
tion representing the internal motions. The Schródinger equation (3.4) is 
the quantum mechanical equation representing this internal motion of the 
hydrogen atom once the equations governing the translational motion of 
the whole system through space have been separated. In this Schródinger 
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equation, the effective mass of the reduced system is denoted u and is 
known as the reduced mass. 
For a system with two particles of mass m, and mp, 





PENES Lo M (3.5) 
mı + m» 
Clearly, if m = m; = m», then u = 5. 
For the case of the hydrogen-like atom, 
mM 
= 3.6 
a (3.6) 


In the limit of an infinite nuclear mass, i.e. M = co, u = me. In practice, 
M is very much bigger than m,. For example, for H itself, M is about 
1836 me. Under these circumstances the reduced mass is very close to me. 
However, as discussed in Secs. 3.7 and 3.11, the small shift implied by a 
correct treatment of the reduced mass is important for astronomical obser- 
vations as it means that hydrogen and deuterium give distinct spectra that 
can be distinguished at moderate resolution. 


3.4 Atomic Units 


Atomic units provide a complete integrated unit system in the manner 
of SI units but with quantities scaled to the dimensions of the atom. In 
atomic units: 
Unit of mass is the electron mass, me = 1.6605402 x 1077 kg, 
Unit of electric charge is the electron charge, e = 1.602188 x 10°C, 
Unit of length, the Bohr radius ao = ine = 5.29177249 x 107" m, 
Planck's constant divided by 27, A = 1a.u. = 1.05457 x 10 ?*J.s, 
Similarly 47€, = 1. 
In this unit system, the atomic unit of energy is known as the Hartree and 
is denoted E;. 1E; = 2Rx = 27.2113661 eV = 4.3597482 x 10 5 J, where 
Ræ is the Rydberg constant (see Sec. 3.7). The atomic unit of time is 
2.41884 x 1077 s and the speed of light c is 137.03599 a.u. 
Atomic units are often denoted a.u. and must be carefully distin- 


guished from the somewhat larger Astronomical Unit, AU, or indeed the 
less well-defined ‘arbitrary units’. 
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3.5 Wavefunctions for Hydrogen 


The Schródinger equation (3.4) can be solved analytically by working in 
spherical polar coordinates, i.e. r = (r, 0, à). In these coordinates the wave- 
function is separable into radial and angular solutions 


vr, 0, à) = Ril) Xin, Q). (3.7) 


The radial solutions, R,;j, can be expressed analytically in terms of 
Laguerre polynomials [see Rae (2002) in further reading]. Figure 3.1 shows 
both the wavefunctions and the probability distribution of the lowest few 
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Fig. 3.1. Wavefunctions (left) and probability distribution (right) for the radial 
coordinate of the hydrogen atom. (T.S. Monteiro, private communication.) 
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Fig.3.2. Wavefunctions for the angular motions of the hydrogen atom which are 
given in terms of spherical harmonics, Y, ,,(0, $). Plots are for the x-z plane except 
for the | = 2, m = 2 plot for which the x-y plane is shown. The plots give the abso- 
lute value of the spherical harmonic as a distance from the origin for each angle; the 
signs indicate the sign of the wavefunction in each region. (S.A. Morgan, private 
communication.) 


radial solutions. An important feature of these wavefunctions is the nodes 
which are the points where the wavefunctions pass through zero. 

The angular solutions for the hydrogen atoms, Yim, are called spherical 
harmonics. These functions are complex for m # 0. Figure 3.2 depicts the 
first few of these spherical harmonics. Readers familiar with chemistry 
will recognise these shapes as the s, p and d orbitals used to represent 
atomic structure and chemical bonds. 


3.6 Energy Levels and Quantum Numbers 


For bound states, the solutions of the Schródinger equation (3.4) occur 
with discrete energies given by the formula 


254 2 
Bee URS (3.8) 


En =~ aab Wa 
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where the zero of energy is taken to be a completely separated electron and 
nucleus with zero kinetic energy. The constant R in Eq. (3.8) is the Rydberg 
constant discussed in the next section. 

Each bound state of the H atom is usually characterised by a set of 
four quantum numbers (n, l, m, sz). Each of these is defined below, along 
with a fifth quantum number s. 


n is the principal quantum number. It takes the values n = 1, 2,3, ..., 00. 
n determines the energy of the atom according to Eq. (3.8). 

l is the electron orbital angular momentum quantum number. The 
actual angular moment is given by A[/(/ + 1)]2. I can take the values 
0,1,2, ..., n — 1. By convention, the values of / are usually designated 
by letters (see Table 3.1). 

m is the magnetic quantum number, so called because it determines the 
behaviour of the energy levels in the presence of a magnetic field. mh 
is the projection of the electron orbital angular momentum, given by l, 
along the z-axis of the system. It can take (2/ + 1) values —l, —1 + 
d cir sco b LE 

s is the electron spin quantum number. The electron spin angular 
momentum is given by A[s(s + 1)]? which, for a one-electron system, 
equals Bh, since an electron always has spin one-half. 

sz gives the projection of the electron spin angular momentum, given 
by s, along the z-axis of the system. This projection is actually Asz. In 
general, s; can take (2s + 1) values given by —s,—s+1,...,s—1,s. 
For a one-electron system, this means s; can take one of two values: 
-i or +4. 


The simplest notations for the various states of H is to denote each state 
by its nl quantum numbers. Thus the ground state is denoted 1s; the first 
excited states are 2s and 2p; the n = 3 states are 3s, 3p and 3d. (See Fig. 3.3.) 
These notations leave the m and s, quantum numbers unspecified since 
these quantum numbers are really only significant for H in the presence 
of an external field. This means that each nl configuration is 2(2/ + 1)-fold 


Table 3.1. Letter designations for orbital 
angular momentum quantum number 1. 


0 1 2 3 4 5 6 7 8 
s pd f gh i kl 
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Fig. 3.3. Schematic energy levels of the hydrogen atom with various series iden- 
tified. [Adapted from P.W. Merrill, Lines of the Chemical Elements in Astronomical 
Spectra (Carnegie Institute of Washington Publications, 1956).] 


degenerate. Such configurations are used to build up basic atomic struc- 
ture (see Sec. 4.4). More complete and complicated spectroscopic notations 
used to denote states of H and other atoms are discussed in Sec. 4.8. 


3.7 H-Atom Discrete Spectra 


The spectrum of the hydrogen atom comes from electrons jumping 
between different levels in the atom. Given that the energy levels depend 
only on the quantum number n, this means that the electronic spectrum 
of the H atom comes from changes in n. The wavelengths, A, for these 
transitions are given by the Rydberg formula: 


1 1 1 1 

X = go Em 2r En) = Ru( sa =e z) ny < na . (3.9) 
This formula was constructed by the Swedish physicist Johannes Rydberg 
(1854-1919) to model the experimental observations of Balmer and others. 
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Of course, the formula agrees perfectly with the energy differences that 
one obtains using the quantum mechanical expression for the energy 
levels given in Eq. (3.8). In this expression, R is the Rydberg constant. 
This constant takes slightly different values for different H-like atoms 
since it incorporates the reduced mass u (see Sec. 3.3). For hydrogen itself, 
Rg is particularly well-determined experimentally; it takes the value 
109677.581 cm !. 

The Rydberg constant for a system of infinite nuclear mass is denoted 
by Ræ and is related to Ry by 


u My 
euge (ema ce ae 1 
Ru = ERs (a) (3.10) 


where R4, = 109737.31 cm~!. Rydberg constants for other species such as 
deuterium, can be derived using the formula above with the appropriate 
reduced mass. 

The spectrum of H is divided into a number of series linking different 
upper levels n with a single lower level n; value. Each series is denoted 
according to its nı value and is named after its discoverer. Table 3.2 sum- 
marises the six main H-atom series. 

The range of each H-atom series is given by the lowest frequency tran- 
sition, between the levels nı and n; = n4 + 1, and the series limit, which is 
the transition between nı and n» = oo. As discussed in Sec. 3.8.1, the series 
limit is not always observable. Table 3.2 gives the spectral region in which 
each series is observed. As the Balmer series lies in the visible region, it 
is particularly easy to observe from Earth. As a result, Balmer lines have 
been particularly important in the study of H-atom spectra. 


Table 3.2. Spectral series of the H atom. Each series comprises the transitions 
ny — n1, where n1 < n» < oo. 


Range/cm ! 


ni Name Symbol Spectral region 15 — n4 41 Ny = oo 
1 Lyman Ly ultraviolet 82257 - 109677 
2  Balmer H visible 15237 - 27427 
3 Paschen P infrared 5532 - 12186 
4 Brackett Br infrared 2468 - 6855 
5 Pfund Pf infrared 1340 - 4387 
6 Humphreys Hu infrared 808 - 3047 
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Within a given series, individual transitions are labelled by Greek 
letters. These letters denote the change in n or An. In this notation: 


An = lisa, 


An = 2 is p, 
An —3is y, 
An = 4is ô, 
An = bis c. 


Thus Lya is the transition between nı = 1 and m = 2, and Hy is that 
between nı = 2 and nz = 5. Greek letters are usually only used for the 
most important transitions with low An. Transitions with high An are com- 
monly labelled by the number n2. Thus, H15 is the Balmer series transition 
between nı = 2 and nz = 15. 

The wavelength of each individual transition can be predicted using 
the Rydberg formula. For example Lya lies at 


a 1\ 30 E 
x- Ru(1 E i) = [Ru = 82258.2cm™, (3.11) 


meaning that 
A = 1.21568 x 10? cm = 1215.68 A = 121.168 nm. 


All hydrogen series transitions between bound states are described as 
bound-bound transitions. Figures 3.4, 3.5, 3.6 and 3.7 give sample H-atom 
spectra recorded in very different astronomical environments. Figure 3.4 
shows an optical spectrum of the B-type star O! B Ori. Absorption by 
the Balmer series up to H14 (n = 14) is clearly visible. Figure 3.5 
shows Lyman series absorption in a shell of gas expanding about a hot, 
Wolf-Rayet star. Figure 3.6 shows higher members of the Paschen series 
recorded in absorption in a spectrum from a soft X-ray transient binary 
star. Figure 3.7 shows infrared emissions due to Pfund and Brackett lines 
recorded in the gas ejected from supernova 1987a. 

The reduced mass factor means that deuterium spectral lines are 
shifted with respect to H-atom spectra. In principle, one can detect D and 
H together (see Fig. 3.8). However, the cosmic abundance of deuterium is 
about 2 x 10~° that of hydrogen. This makes it difficult to measure abun- 
dance ratios directly using neighbouring transitions as the H transitions 
are likely to be optically thick if those of D can be observed. Figure 3.9 
shows high-resolution simultaneous observations of heavily red-shifted 
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Fig. 3.4. Spectrum of B-type star ©! B Ori showing Balmer series absorptions 
recorded at the Anglo-Australian Telescope. (X-W. Liu, private communication.) 
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Fig. 3.5. Satellite spectrum recorded from shell nebula GRB 021004. The doublet 
structure in the Lyman series absorptions is caused by Doppler effects in the shell. 
[Adapted from N. Mirabel et al., Astrophys. J. 595, 935 (2003).] 


H and D Lya lines. These data need to be interpreted very carefully since 
H absorptions with a different Doppler shift to the main peak, such as 
the peak on the left of the figure, can be confused with absorptions by 
D. This particular spectrum was taken as part of a systematic study of 
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Fig. 3.6. Spectrum soft X-ray transient (SRT) binary star GRO J1655-40 recorded 
with the Anglo-Australian Telescope, showing higher members of the Paschen 
series. P13, P15 and P16 are blended with the Call triplet transitions discussed 


in Sec. 6.5. [Adapted from R. Soria, K. Wu and R.W. Hunsted, Astrophys. J. 539, 445 
(2000).] 
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Fig.3.7. Infrared spectrum of supernova 1987a taken using the Anglo-Australian 
Telescope 192 days after the initial explosion. Strong Brackett and Pfund emissions 
from atomic hydrogen are clearly visible. This spectrum gave the first-ever detec- 
tion of the molecule H7 outside the solar system; its spectrum is modelled by the 
solid line. The structure in the spectrum about 3.2 um is due to telluric effects in 


this region. [Analysis reported by S. Miller et al., Nature 355, 420 (1992).] 
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Fig. 3.8. Lyman absorption spectrum for hydrogen and deuterium recorded 
using the FUSE satellite along a line of sight towards B-type star HD 191877. 
[Adapted from C.G. Hoopes et al., Astrophys. J. 586, 1094 (2003).] 


the absorption of light from distant quasars which has firmly established 
the primordal D to H ratio, one of the fundamental parameters of the Big 
Bang. It should be noted that since deuterium is heavier, its lines should 
have smaller Doppler width. 


Worked Example: H-atom Ha emission occurs at 15237 cm~!. At what 
wavenumber is the corresponding transition for the D atom? 
The Rydberg formula for Ha gives: 


1 1 
15237 cm™! = hi — ;) (8.12) 


The Rydberg constant for deuterium, Rp, is given by 


puc BP (3.13) 
HH 
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Fig.3.9. Hubble Space Telescope spectrum of Lya absorption against quasar QSO 
1009+2956 which lies at redshift z = 2.5. The continuous lines give the results of 
models which fit several Ly series absorptions. The vertical lines give line centres; 
the three lines near 4258 A are due to D; other absorption features are all due to H. 
[Adapted from S. Burles and D. Tytler, Astrophys. J. 507, 732 (1998).] 


where uy and up are the reduced masses of the hydrogen and deuterium 
atom, respectively. 


HH  Mudcm, Mom, — (My + me)M 


D 
ee ES ED 00 3.14 
uo — Mum, Mp-4m, (Mp4 m)Mg (3.14) 





where the masses My = 1836.1 m, and Mp = 3670.4 m, have been 
used. Using these numbers gives Ha for D as 15233 cm !, compared to 
15237 cm ! for H. 

Although the Rydberg constant is known very precisely, the wave- 
lengths of various transitions in hydrogen do not coincide completely 
with the values predicted by the Rydberg formula. It should be noted 
that this formula is only correct within certain assumptions which neglect 
small effects due to relativity and the finite size of the hydrogen nucleus. 
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For this reason, the measured wavelengths should be used if high accuracy 
is desired. 


3.8 H-Atom Spectra in Different Locations 
3.8.1 Balmer series 


Balmer series lines are the most studied H-atom lines since they are in the 
visible. Indeed, Ho emissions at 6563 À can be clearly seen with the naked 
eye as the red light surrounding the Sun during a total eclipse. In fact, 
Ha and H lines were labelled C and F respectively in Fraunhofer's solar 
spectrum (see Fig. 1.1). 

Stellar envelopes are fairly high density environments. This means 
that the population of the different atomic energy levels is thermal and 
given by the Boltzmann distribution: 


P= Ee(-55:). (3.15) 


where P; is the population of the ith level given as a proportion of the 
atoms in level i. P; can therefore take values between zero, meaning no 
atoms in level i, and one, meaning all the atoms are in level i. In Eq. (3.15), 
gi is the degeneracy (statistical weight) of level i, and the energy above the 
ground state is AE;. Q is the partition function which ensures that the sum 
over all populations is unity. k is Boltzmann's constant. When analysing 
spectra, it is useful to remember that k = 0.695 cm ! K 1. 

However the temperature, T, varies between stars, which results in 
very different spectral characteristics. Since the centre of the star generates 
bright radiation, spectral lines are seen in absorption. The strength of these 
absorptions depends on the spectral type of the star. 

Hydrogen Balmer lines are strongest for A0 stars which have a tem- 
perature of about 10000 K. Few other transitions are seen in these stars. 
Cooler stars have less population in the n = 2 level of H, so the Balmer 
line absorption diminishes. Indeed, in the coolest stars (T « 5000K), 
molecules form and little atomic H remains. For stars significantly hot- 
ter than 10000 K, H atoms become increasing ionised and the strength of 
the Balmer series drops again. 

This discussion illustrates a general point. Atoms lose electrons to 
become ions. As the environment gets hotter, the degree of ionisation 
increases. Since this is essentially a thermal effect, for every atom, there is a 
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particular ionisation stage which is dominant at a particular temperature. 
At this temperature, the spectral lines of this ion will be at their strongest. 
Furthermore, successive ions of particular atoms have spectra quite unlike 
each other; for example H ions are protons which have no spectrum at all! 

Figure 3.10 shows the temperature distribution of iron ions. The 
Roman numerals in this figure are used to designate ionisation stages. 
Thus, Fe I corresponds to neutral iron Fe; Fe II to singly ionised iron Fe*; 
Ferr to Fe?*, and so forth. Fe xxvI denotes Fet which is H-like, i.e. it 
has one electron. Fe XXVII is a bare iron nucleus. Strictly speaking, Roman 
numerals label a spectrum, but they are often used to indicate an ion. 

For hydrogen, it is usual to talk of H1 regions where H is neutral. 
Often, no optical H-atom spectrum is seen from these regions. Conversely, 
in H II regions, the H atoms are ionised. As discussed below, H I recombi- 
nation spectra are seen from H II regions. 

Spectra are sensitive to pressure effects. The width of spectral lines 
depends on pressure, and information on the luminosity class of a star can 
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Fig. 3.10. Abundance of iron ions with different ionisation stages as a function 
of temperature. These results are taken from the coronal model of J.M. Schull. 
[Adapted from A. Dalgarno and D. Layzer (eds.), Spectroscopy of Astrophysical Plas- 
mas (Cambridge University Press, 1987).] 
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be obtained by measuring these line-widths. If the pressure is low, higher 
members of the Balmer series can be seen (see Fig. 3.4). However, the com- 
plete series is never observed in stars and the spectrum shows a discon- 
tinuity known as the Balmer jump (see Fig. 3.11). The jump occurs 
at the point where continuous (bound-free) absorption switches on. The 
position of the Balmer discontinuity shifts to longer wavelengths at higher 
pressure. 

To understand why, it is necessary to consider the physical size of the 
H atom. The approximate radius of the electron orbit in the Bohr atom, rn, 
is given by the formula: 


Aregh? n? 
2 0 
=n 27 =7 ao, (3.16) 





where for H(1s), rı = 1 and ap = 0.529 x 10! m. At high pressure, the 
mean atom-atom distance will be smaller than some r;, and higher orbits 
will not be found since they will be destroyed by atom-atom collisions. 
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Fig. 3.11. Model spectrum of an A5-type star showing both the Balmer and 
Paschen discontinuities. (R.J. Sylvester, private communication.) 
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Worked Example: The Sun has a number density, N, of about 10" cm ?. 
What is the highest H-atom 7 level that one would expect to find? 
The approximate radius of an electron orbit in the H atom is: 


ry = n? x 0.529 x 1079 m. 


Therefore the approximate volume of an H atom in state n is 


Vn & er ns Enn x 1.48 x 107 md, 
3 3 
On the surface of the Sun, N = 10!” cm? = 109 m ?. This means that each 
atom occupies a maximum volume of about 107? m°. 
Assuming that the highest n value corresponds to the maximum volume 
allowed for each atom gives 


4 
37" x 148 x 10?! = 1072, 


which gives n 2:16 as the highest level that one would expect on the Sun. 


3.82 Lyman series 


The Lyman series is expected to be strong in absorption spectra of hot 
stars which have significant ultraviolet continuum. However, there are a 
number of technical issues that need to be considered as it is not possible 
to observe such ultraviolet transitions from the ground. 

So far, only a few satellites with the capability of recording spectra this 
far into the ultraviolet have been launched. The Hubble Space Telescope, 
for example, does not possess this capability; its best region for spec- 
troscopy is given by \ > 1300 A, whereas Lya lies at 1216 A. The highly 
successful International Ultraviolet Explorer (IUE) satellite, which oper- 
ated for 19 years starting from 1978, covered Lyo but no higher Lyman 
lines. However, this situation changed with the launch of NASA's FUSE 
(Far Ultraviolet Spectroscopic Explorer) satellite in June 1999 (see Fig. 3.8). 
FUSE covers A — 900-1200 À at good resolution and sensitivity. FUSE's 
range thus extends below the Ly series limit at 912 À. 

Lyman absorption comes from the H-atom ground state. It is thus 
dominated by a strong interstellar component. See Fig. 3.9 for an exam- 
ple. Interstellar H atoms are all in their ground states, and are thus only 
sensitive to Ly wavelengths. This complicates the interpretation of any 
Ly emission spectrum. 
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Fig.3.12. Schematic Lyman a emission line showing geocoronal emission and, in 
the wings, absorption against starlight from the interstellar medium. 


Spectra taken in the locality of the earth are contaminated by strong 
geocoronal Lyo emission (see Fig. 3.12). Such emissions are present in IUE 
spectra, which were taken in geosynchronous orbit, and also some Hubble 
Space Telescope spectra, for which Fig. 7.7 gives an example. Figure 10.13, 
which was recorded with FUSE, shows the clear signature of geocoronal 
Ly emission. 

Lya can be used as a probe of the Universe at earlier epochs. Quasars 
emit ultraviolet light with a quasi-continuous distribution. Absorption of 
this light by H atoms via the Lya transition can be monitored with differ- 
ent epochs being distinguished by their differing redshifts. 

Figure 3.13 shows a high-resolution spectrum of a high redshift 
quasar. The strong and broad Lyman a emission line is observed at a 
wavelength of 5622 Á, indicating that the quasar is at a redshift z — 3.625. 
The spectrum of the quasar at shorter wavelengths is eaten away by a 
great number of narrow absorption lines. Most of these are single Lya lines 
formed in gas located between the quasar and us; they therefore appear at 
many discrete redshifts between z — 3.625 and 0. These absorption lines 
are collectively known as the Lyman a forest. They are a powerful probe of 
physical conditions in galaxies and the intergalactic medium at early times 
in the evolution of the universe. 


3.8.8 Infrared lines 


The H-atom series higher than Balmer absorb and emit in the infrared. 
These can be seen in a number of locations. For example, Fig. 3.7 shows the 
spectrum of the outer expanding shell of supernova 1987a taken 192 days 
after the initial supernova explosion. This spectrum was used to report the 
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Fig. 3.13. Spectrum of the high redshift QSO O14224231, obtained with the 
High Resolution Echelle Spectrograph on the Keck telescope in Hawaii showing 
numerous Lya absorptions by the intervening interstellar medium. These absorp- 
tions are known collectively as the Lyo forest. (M. Pettini, private communication.) 


discovery of the H5 molecular ion in the expanding gas, however it also 
shows strong Bro emission and several Pfund series lines. 


3.9 H-Atom Continuum Spectra 
3.9.1 Processes 


The continuum of a proton and an electron is not quantised. Therefore an 
H atom in its 1s state can be ionised by any photon with A < 912A, i.e. 
beyond the Lyman series limit. 


H(1s) +hv H* +e. (3.17) 


This process is called photoionisation; it is a bound-free transition. Simi- 
larly, the Balmer continuum is observed for \ < 3646 A (see Fig. 3.11). The 
probability of photoionisation is greatest when near threshold and drops 
away as A decreases. See Fig. 3.14 which illustrates bound-bound and 
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Fig. 3.14. Schematic stellar spectrum showing the hydrogen Balmer lines and 
continuum. There is a large departure from a black body spectrum at ultraviolet 
wavelengths below the Balmer series limit. 


bound-free (photoionisation) spectra for Balmer lines absorbing against 
the black body radiation curve of a star. 
The reverse process, 


Ht +e > H(nl) thy, (3.18) 


is called radiative recombination; it is a free-bound transition. The probabil- 
ity that an electron will emit a photon as it passes the proton is low, so this 
is a very unlikely process. However, the chances increase if the electron is 
travelling slowly, i.e. near threshold. 

Radiative recombination is essential, if inefficient, for cooling. For 
example the recombination era, which started about 300000 years after the 
Big Bang and lasted for about 1 million years, was the period in which 
neutral H and He formed in the early Universe. This is the earliest epoch 
of our Universe for which it is theoretically possible to record spectra. 


3.9.2 H-atom emission in HII regions 


Figure 3.15 shows many higher Balmer lines recorded in Orion. These 
emissions come from H II regions — ionised regions around hot stars where 
there are strong ultraviolet fluxes and gases with temperatures of about 
10000 K. Under these conditions any H atom will be rapidly photoionised. 

Although the recombination of an electron and a proton to produce 
an H atom is an inefficient process, it still happens continually in H 1 
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Fig. 3.15. Spectrum of the Orion nebula (M42) recorded using the ESO 1.52m 
telescope in Chile. The vertical scale gives the observed flux in units of 10-16 ergs- 
cm2. s71. AT! - arcsec™?. Also shown are are polynomial fits to the continuum 
bluewards and redwards of the Balmer jump at 3646 A. These can be used to 
determine the temperature of the nebula. [Adapted from X.-W. Liu et al., Astrophys. 
J. 450, L59 (1995).] 
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Fig. 3.16. Recombination spectra start with the emission of a continuum photon 
followed by a cascade of emission lines known as recombination lines. 


regions. This produces continuum emissions and a small proportion of 
neutral H atoms. Typically, about 176 of the protons are in the form of neu- 
tral atomic H. Recombination can occur at different levels of the atom. If 
the atom is formed in an excited level, it can decay by a series of emis- 
sion lines (see Fig. 3.16). Since the emission lines follow in a series, the 
process is sometimes referred to as a cascade. Different routes are possi- 
ble and which is preferred depends on the relative values of the different 
Einstein A coefficients. The ratio of the A coefficients, which determine the 
relative importance of the competing pathways, are known as branching 
ratios. 
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H-atom lines from H II regions occur as a result of recombination from 
the continuum. They are called recombination lines. 

Collisional excitation of H in lower energy levels is not important in 
H 1 regions — a ground state atom is much more likely to be photoionised. 
Populations of different levels are therefore determined by radiative pro- 
cesses and not a Boltzmann distribution. 

Recombination and collisional excitation thus form competing mech- 
anisms for driving emission spectra. In H II regions, the most common 
situation is that while the H-atom spectra are dominated by recombination 
lines, the spectra of the other atomic species present are driven by collisional 
excitation. This means that by careful study of both types of emission spec- 
tra, one can obtain a measure of both the density and the temperature of 
the nebula, as well as the strength of the radiation field which drives the 
ionisation process. This means that by studying spectra one can get infor- 
mation on both the nebula and the central star that powers it. 


3.10 Radio Recombination Lines 


As discussed in Sec. 3.8.1, the Balmer series is truncated in the atmo- 
spheres of stars by pressure effects. This means that the long wavelength 
transitions between higher levels cannot be observed. Densities in plan- 
etary nebulae are much lower, about 10?—10* cm ?, compared to about 
10 cm ? at the Sun's surface region, or about 10? cm ? in the Earth's 
atmosphere at sea level. At very low densities, bound states with large 
n can exist. The Bohr radius for an H atom with n = 137 is approximately 
1 um, meaning that atoms in this state can only survive at densities signif- 
icantly less than 10? particles - cm ?. 

Recombination lines, transitions between neighbouring high values 
of n, can be detected at radio wavelengths. These transitions are labelled 
using a high n version of the series notation described in Sec. 3.7. 


Transitions n — n + 1 are called ‘Hno’, 





Transitions n — n + 2 are called ‘Hng’, 


and so forth. It should be noted that ‘H’ here stands for hydrogen and not 
Balmer. The transition to state n is strongest when An = 1, i.e. the Hna tran- 
sition (see Fig. 3.17). Table 3.3 gives wavelengths and lifetimes for some 
radio-frequency transitions. 
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Fig. 3.17. Millimeter wavelength hydrogen radio recombination lines observed 
toward the very luminous, peculiar star 7 Carinae using the SEST telescope. Note 
that the intensity is given by the antennae temperature, T4, which is conventional 
in radio astronomy. [Reproduced from P. Cox et al., Astron. Astrophys. 295, L39 
(1995).] 


Transitions with very high n have been observed, such as H766a, 
where the size of the H atom involved approaches 0.1 mm. Certain lines 
are observed more than others. For example, the transitions H109a and 
H1375 both have wavelengths near 6 cm and can be observed at the same 
time. The H166a line is near 21 cm, where many detectors are available for 
observing the so-called 21 cm transition (see Sec. 3.14). 
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Table 3.3. Wavelengths, A, and Einstein A coeffi- 
cients for some typical radio frequency hydrogen 
atom transitions. 


H50a H508 H50y H150a 


À 59mm 30mm 21mm 15.5cm 
A/s 18 9.4 6.1 0.008 


Radio recombination of H atoms can lead to population inver- 
sion. Population inversion is one of the conditions necessary for maser 
(microwave amplified stimulated emission of radiation) action (see 
Sec. 2.4), a process physically similar to that of a laser. In fact, H atoms 
can mase in dense HII regions. Such masing has been observed from 
MCW 349 via the na lines. In this location, maser action has been observed 
for all transitions with 7 € n < 90. However, molecular masers are much 
more common and important than those seen in atomic sources. Molecular 
masers will be discussed in Sec. 10.5. 


3.11 Radio Recombination Lines for Other Atoms 


An important feature of spectroscopy is that any atom, ion or molecule has 
a unique spectrum by which it can be readily identified. An exception to 
this general rule are the radio recombination lines. When a single electron 
is promoted to a state with a very high principal quantum number, n, the 
electron experiences a potential due to the ion core which feels like that 
due to a single point charge such as a proton. This is because this electron 
is so far away from the nucleus that the nucleus and the other electrons 
appear to it as if they occupy a single point. 

Under these circumstances the energy levels of the outer electron 
satisfy the H-atom formula, Eq. (3.8), with an effective nuclear charge, 
Ze, equal to one. Note that a more sophisticated treatment uses quantum 
defect theory and replaces n? in Eq. (3.8) with (n — u?, where u is known 
as the quantum defect. Quantum defect theory is discussed in Sec. 6.1. 

Within the confines of the H-atom formula, the only difference 
between the high n energy levels, and hence the radio recombination spec- 
trum, of different atoms arises from the different nuclear masses. These 
different masses give rise to reduced masses, and Rydberg constants, 
which depend on the atom in question in a fairly simple fashion. The 
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reduced-mass factor differs because the nuclear masses, M, differ between 
atoms, giving a scaled Rydberg formula 


Ro M 


Debo Te, 
^ n? Mm, 


(3.19) 

Figure 3.18 gives a sample radio recombination spectrum showing 
lines from hydrogen, helium and carbon. Radio recombination lines of 
H and He atoms have been observed from many ionised locations. Lines 
of heavier atoms such as C have also been detected. Traces of heavier ele- 
ments, which could be S or Mg, can also be detected sometimes. However 
as the atoms get heavier, the frequencies of the recombination lines con- 
verge (see Table 3.4), making it very difficult to conclusively identify the 
species involved. Note that in Fig. 3.18 the relative strengths and profiles 
of He and C lines differ. This is because recombination and hence the emis- 
sion lines arise from different regions. Helium is only ionised in an inner 
zone of the H II region which has energetic photons, so it is hot. As a result, 
the line is strongly Doppler broadened. Conversely, carbon is ionised in 
cooler regions where even the H is neutral (see Fig. 3.19), so the line is 
much narrower. 

How does the energy balance work in the H II regions from which the 
recombination emissions are observed? The heating comes from incoming 
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Fig. 3.18. Radio recombination lines of hydrogen, helium and carbon observed 
in the gaseous nebula W3. [Adapted from A. Dalgarno and D. Layzer (eds.), 
Spectroscopy of Astrophysical Plasmas (Cambridge University Press, 1987).] 
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Table 3.4. Frequencies, f, for the 50a radio recom- 
bination lines of hydrogen, helium, carbon and an 
atom with an infinite nuclear mass. The frequencies 
of the heavier atoms converge on this value. 


H50a He50a C509 M = œ 50a 
f/MHz 51,072 51,092 51,097 51,099 
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Fig.3.19. Layers of ionised gas about a central hot star. The figures in eV give the 
energy required to maintain the ions indicated, and hence the maximum energy of 
the ultraviolet radiation, found outside each sphere. 


radiation, usually from a central star, which is ‘captured’ by photoionisa- 
tion in a bound-free transition: 


H+hv>H*+e. (3.20) 


This is balanced by outgoing radiation, generally at longer wavelengths. 
This comes in the form of continuum photons formed by recombination, 
the reverse process to Eq. (3.20), and the resulting cascade emissions. It 
also comes from emissions from heavier atoms and from ‘bremsstrahlung’. 
Bremsstrahlung, a German word meaning 'braking radiation', is gener- 
ated by electrons changing velocities as they pass close to charged nuclei 
(see Fig. 3.20). Bremsstrahlung photons are not quantised; they can thus 
form a continuous distribution in frequency. Such transitions are classified 
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Fig. 3.20. Bremsstrahlung (‘braking radiation’) is emitted when a fast electron is 
slowed down by passing a charged nucleus. 


as being free-free transitions as both the upper and lower states are not 
bound. 


3.12 Angular Momentum Coupling in the 
Hydrogen Atom 


One-electron atoms such as hydrogen contain several sources of angular 
momentum: 


e Electron orbital angular momentum /; 
e Electron spin angular momentum s; 
e Nuclear spin angular momentum i. 


The nuclear spin arises from the spin coupling of the various nucleons: 
protons and neutrons both have an intrinsic spin of a half. For atomic spec- 
tra it is not necessary to consider the spin of individual nucleons, just their 
total spin represented by the quantum number i. For many-electron atoms 
however, one has to consider the / and s quantum number for each electron 
in the system. This situation will be considered in Sec. 4.7. 

As in classical mechanics, only the total angular momentum is a con- 
served quantity. It is therefore necessary to combine angular momenta 
together. This is best done by adding two angular momenta at a time. The 
order in which this is done is referred to as a coupling scheme. The choice 
of coupling scheme usually reflects the strength of the actual couplings: 
the strongest couplings are considered first. 

For hydrogen, the usual coupling scheme is to combine / and s to give 
the total electron angular momentum j. These are added vectorially as 


l+s=j. (3.21) 
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One then combines the total electron and nuclear spin angular momenta 
to give the final angular momentum f: 


jti-f. (3.22) 


To do this one needs to know the rules for addition of angular momenta. 
This is based on the use of vector addition. Note that angular momentum 
is a vector as it has both a magnitude and an orientation. 

In classical mechanics, adding vector a and vector b gives a vector c, 
whose length must lie in the range 


la—bl<c<a+b, (3.23) 


where a, b and c are the lengths of their respective vectors. This is some- 
times known as the triangulation condition since the lengths of the vectors 
must allow them to form a triangle. 

In quantum mechanics a similar rule applies except that the results 
are quantised. The allowed values of the quantised angular momentum, 
c, which arise from adding a and b, span the range from the sum to the 
difference of a and b in steps of one: 


c — |a —b|,|a - b| -1,...,a-b — 1,a b. (3.24) 


For example, add the two angular momenta l; = 2 and I, = 3 together to 
give L. This is usually written in vector form 


L-h-tb, (8.25) 
and the result is: 


L= lly — DL lh — 5| 1,...,h-b - 1h +h, 
—1,2,3,4,5. 


3.13 The Fine Structure of Hydrogen 


Electron spin arises as part of a relativistic treatment of quantum mechan- 
ics. Relativistic effects couple electron orbital angular momentum, l, and 
electron spin, s, to give the so-called fine structure in the energy lev- 
els which are split according to the value of the total electron angular 
momentum j. 
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Table 3.5. Fine structure effects in the hydrogen atom: splitting of 
the nl orbitals due to fine structure effect for | = 0,1,2,3. The result- 
ing levels are labelled using H atom, and the more general spectro- 
scopic notation of terms and levels (see Sec. 4.8). 


Configuration l s J Hatom Term Level 
ns 0 l l nsi n?S n?$, 
2 2 

np 1 i L3 mp1, npa n?po  g?po,g?po 

2 2 

nd 2 l 3,3 nds, nds n2D n?D;, n?D; 

1 5 7 2go 2go 2go 

nf 5 3,5 nfs,nf n*F n*F2,n*^F 
3 5 55 prem 3 z 





For hydrogen, s = 4 so that, except for the | = 0 case, j = 1 + 5 (see 
Table 3.5). This table labels the resulting levels with the common H-atom 
notation nl;, where I is given by its letter designations, s, p, d, etc., and 
by spectroscopic notation for which labels of the €5^U[; are used. A full 
discussion of spectroscopic notation can be found in Sec. 4.8. 

Table 3.5 shows the fine structure levels of the H atom. This table 
shows that the states with principal quantum number n = 2 give rise to 
three fine-structure levels. In spectroscopic notation, these levels are 2 ?S p 
2 Pi and 2 "PS. 

So far the discussion on H-atom levels has assumed that all those with 
the same principal quantum number, n, have the same energy. In other 
words, the energy does not depend on l or j. This is not correct: inclusion 
of relativistic (or magnetic) effects split these levels according to the total 
angular momentum quantum number j. This splitting, called 'fine struc- 
ture’, has been well-studied in the laboratory. An even more subtle effect 
called the Lamb shift, which is due to quantum electrodynamics, can also be 
observed. Values of these splittings for the n = 2 levels are given in Fig. 3.21. 
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Fig. 3.21. Splitting in the n = 2 levels of atomic hydrogen. The larger splitting is 
the fine structure and the smaller one the Lamb shift. 
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For hydrogen, fine-structure and Lamb-shift splittings are too small 
to be important for most astronomical applications. The fine structure is, 
however, of great importance for complex atoms and will be discussed 
further in Chapters 4 and 6. 


3.14 Hyperfine Structure in the H Atom 


There is one more source of angular momentum in the H atom which has 
not yet been included. This is the nuclear spin, i; for H, i — 1. Coupling i 
to the total electron angular momentum, j, gives the final angular momen- 
tum, f [see Eq. (3.22)]. For H this means 





f-jct I ; (3.26) 


The ground state of H is 1s; or "8; and has j — 5. This means that 
nuclear spin coupling can split this state into two leve with f = 0 or 1. 
There is a very small, 6 x 10? eV, splitting between the lower f = 0 
and higher f — 1 levels of H caused by magnetic effects. The f — 0-1 
transition between these levels has a frequency of 1420.406 MHz which 
corresponds to a wavelength of 21 cm. The 21 cm line is probably the sin- 
gle most important line in astronomy. It is used to map H-atom densities 
throughout the ISM (see Fig. 3.22 for example). 
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Fig. 3.22. 21-cm-line profiles observed with the Very Large Array for two galac- 
tic lines of sight recorded as part of a study which constructed a face-on galactic 
map of the H II region complexes. The vertical axis of the spectra represents the 
observed line-to-continuum intensity, which is equivalent to exP(— Tine). [Adapted 
from M.A. Kolpak et al., Astrophys. J. 582, 756 (2003).] 
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The 21 cm line is a powerful tool because: 


(1) Displacement of the line gives the line-of-sight velocity; 

(2) Its intensity gives the number of atoms. Note that the line is very weak, 
its Einstein A coefficient is 2.9 x 10? s! which corresponds to a life- 
time of 10 million years, so it is always optically thin; 

(3) Theline profile can be used to deduce the temperature of the gas. Thus, 
for example, Fig. 3.22 suggests that location C23.437-0.207 is much 
cooler that G43.1724-0.006, where the lines appear to be significantly 
Doppler broadened. 

(4) The Zeeman splitting of the transition can be used to measure the mag- 
netic fields. 


Note that the line is seen in both emission and absorption and its analysis 
can be complicated by the presence of several clouds along any particular 
line of sight. 


3.15 Allowed Transitions 


For hydrogen, transitions which correspond to any change in the principal 
quantum number, An, are allowed. However transitions are not observed 
between all states of the H atom or indeed complex atoms. Transitions 
are governed by selection rules which determine whether they can occur. 
The selection rules can be derived rigorously using the rules of quan- 
tum mechanics [see Schiff (1969) in further reading], but will simply be 
stated here. 

Strong transitions are driven by electric dipoles. Weaker transitions, 
driven by both electric quadrupoles and magnetic dipoles, are observed 
astronomically (see Sec. 5.2), but are not important for hydrogen. Rigorous 
quantum mechanical analysis shows that for electric dipole transitions in 
hydrogen: 





As = 0, (for H this is always satisfied as s = i for all states); 
Aj =0 pels 1; 
Ami = 0, t1. 








The last selection rule on mj is only important in a magnetic field. 
Thus if one considers the Ha transition which corresponds to n = 2—3, 
the Al = +1 condition shows that not all transitions between sub-levels 
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occur. Specifically, the transitions 2s -3p, 2p -3s and 2p -3d are allowed 
whereas the transitions 2s —3s, 2p - 3p and 2s -3d are not allowed. 

If fine structure effects are considered, then the selection rules can 
give further constraints. Considering only the Ho transitions designated 
allowed above, the selection rule on j shows that 


2s 17 3p; is allowed; 
-3p; is allowed; 
2p; -3d; is not allowed; 
- 3s; is allowed; 
- 3d; is allowed; 
2p; -3s, is allowed; 
- 3d; is allowed; 


- 3d; is allowed . 


3.16 Hydrogen in Nebulae 


Hydrogen atom emissions in HII regions and planetary nebulae are very 
similar but the latter are generally brighter, which means that more weak 
line emissions can be observed. In particular, lines belonging to the Balmer 
series are often seen strongly in emission. Indeed, the characteristic red 
colour seen in many nebulae comes from Ho. 

Balmer or other spectral series are obtained from excited atoms spon- 
taneously emitting photons. Every excited state has a half-life 7, similar to 
that encountered in radioactive decay, which is related to the strength of 
emission. Thus excited states which decay only by weak line emission are 
long-lived and those which decay via strong transitions are short-lived. 
However, most excited states can emit to more than one other state. 

The lifetime of excited state i is given by 


-1 
Ti = (z A) , (3.27) 
J 


where A;; is the Einstein A coefficient (see Sec. 2.2). 

Lifetimes for allowed atomic transitions are short, perhaps a few times 
10^? s. Table 3.6 gives some examples for the H atom. A glaring exception 
in Table 3.6 is the lifetime of the 2s level of H. This state has a lifetime of 
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Table 3.6. Lifetimes, 7r, for decay by spontaneous emission for 
low-lying excited states of the hydrogen atom. 


Level 2s 2p 3s 3p 3d 
T/S 014 16x10? 16x10" 54x10? 23x107 


2s 





ls 


Fig. 3.23. Decay of the metastable 2s state of hydrogen giving two continuum 
photons. 


~0.14 s, i.e. it lives 10? times longer than the 2p state. This is because the 
transition 2s — 1s is strongly forbidden. The 2s state is metastable which 
means that on the atomic scale, it is long-lived. 

So how does the 2s state decay? By the process of two-photon emis- 
sion, which is an inefficient process and in this case has an Einstein A 
coefficient of 7s~' which can be compared to A(2p — 1s) = 6.3 x 105 s 1. 
The combined energy of the photons emitted must correspond to the 
energy difference E(2s) — E(1s) but the photons themselves can take any 
energy within this constraint (see Fig. 3.23). The photons thus appear as 
continuous emission radiation. Indeed the two-photon decay of the H 2s 
state is responsible for approximately one half the continuum emission 
observed from H II regions. 


Problems 


3.1 Give an expression for the energy levels of the hydrogen atom in terms 
of the Rydberg constant Ry. Assuming a value Ry = 109677.58 cm 1, 
derive a wavenumber for the Lya transition of atomic hydrogen 
in cm. Explain why the Rydberg constant, R4, = 109737.31 cm 1, 
is more appropriate than Ry for a heavy one-electron atom. Hence 
obtain an estimate for the wavenumber of the Lya transition of 
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hydrogen-like iron, Fe?*. From what astronomical environments 
would such transitions occur and how might they be observed? 

A proton and an electron recombine to form atomic hydrogen in its 
4p state. At what wavelengths will recombination lines be observed? 
Label each wavelength by its standard series notation. How would the 
observed emissions differ if the atoms had recombined to the 4s level? 
Hydrogen Ha has a rest wavelength of 6564.71 A. At what rest wave- 
length would you expect to observe deuterium Ha? What telescope 
resolution (R = à) would be required to resolve this difference? 
How does this compare to the resolution required to resolve H and 
D Lya transitions? Given a high enough resolution telescope, what 
other problems do you anticipate in obtaining the B abundance by 
observing a single spectral line such as Ha? The mass of H is 1836.1 
and D is 3670.4 in atomic units. In these units the mass of an electron 
is one. 

A typical star has a number density, N, of about 1016 cm~’, while in 
an H II region the number density is more typically 10* cm ?. In each 
case what is the highest level of atomic hydrogen that is likely to be 
occupied? State any assumptions made in obtaining this answer. The 
partition function for atomic hydrogen can be written as: 


lore) 5 Ry 
Z= 2 exp ET . 
This series cannot be summed as it diverges. Can you suggest why, in 
practice, the partition function for the H atom is finite? 

Use the Rydberg formula to obtain the wavelength of the 80a radio 
line for an atom of infinite mass. Hence, taking the mass of the hydro- 
gen nucleus to be 1836.1 electron masses, obtain the frequency of the 
80a transition of atomic hydrogen. What resolving power would be 
required to separate the two transitions? 


CHAPTER FOUR 


COMPLEX ATOMS 


‘The truth is rarely pure and never simple.’ 


— Oscar Wilde, The Importance of Being Earnest (1895) 


4.1 General Considerations 


Consider an atom with N electrons and nuclear charge (atomic number) Z 
It is straightforward to write down the non-relativistic Schródinger equa- 
tion for this system: 


BEL accen 


EACUS 4T €ofj ur urn 4neo|r; — r| ril 
X Witt, E ,tn) =0, (4.1) 


where 1; is the coordinate of the ith electron, with its origin at the nucleus. 

The first sum in Eq. (4.1) contains a kinetic energy operator for the 
motion of each electron and the Coulomb attraction between that electron 
and the nucleus. The second summation contains the electron-electron 
Coulomb repulsion term. The Coulomb repulsion between pairs of elec- 
trons means the above equation is not analytically soluble, even for the 
simplest case, the helium atom, for which N — 2. This is how all atoms 
with two or more electrons justify their label ‘complex’. To make progress 
on understanding these systems it is therefore necessary to introduce 
approximations. 
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4.2 Central Field Model 


The easiest simplification of the complex atom problem is to try and regain 
a single particle problem, which is similar in spirit to the hydrogen atom 
problem and for which relatively easy solutions can be found. In general 
such solutions depend on allowing the electrons to move in a potential 
which does not depend on their angular position about the nucleus. Such 
potentials, which include the hydrogen atom potential, generate a central 
field since the force acting on each electron only depends on its distance 
from the nucleus at the centre. 

Let us assume that each electron moves in its own angle-independent 
(central) potential given by V;(r;) for electron i. Within this model the 
Schródinger equation (4.1) can be separated into N single electron equa- 
tions. This gives a simplified Schródinger equation for the motion of the 
ith electron: 


hz 
"om Vi + Viri) oir) = Eier). (4.2) 


Using this model, the total energy of the system is given by the sum of 
single electron energies 


E= 2. E;. (4.3) 


The solutions of Eq. (4.2), ¢i(r;), are known as orbitals. Since much 
of the discussion on atomic structure is done in terms of such orbitals, it 
is important to emphasise that the exact solution to complex atom prob- 
lems cannot be written as products of orbitals and that orbitals only exist 
with a central field or independent particle model. For this reason, this 
model is also often known as the orbital approximation. 

Within the orbital approximation, each atomic orbital can be written 
as the product of a radial and an angular function, similar to hydrogen 
[see Eq. (3.7)]. 


oir) = Rna; (ri)Ynn(0;,0;) ; (4.4) 


The angular part of each wavefunction is independent of the other elec- 
trons and is therefore simply a spherical harmonic Y;,, (0;,9;). 

So far no attempt has been made to specify the form of the central 
potential, V;(r;). The standard choice is to consider the motion of each 
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electron in the average field of all the other electrons. With this choice 


— Ze? e 
Vir) = —————À), 4.5 
iri) Arcor; F > ul = rj ) ( ) 


where the notation (...) represents an average. With this form of the cen- 
tral potential, the radial part, R,,; (rj), of orbital (4.4) depends on all the 
other electrons. It therefore cannot be obtained in analytic form and has 
to be evaluated numerically using computers. As V also depends on the 
orbitals of the other electrons, the usual method of solution is iterative. 
An initial guess for the orbitals is used to generate the potential, solving 
Eq. (4.2) gives improved orbitals which then give a new potential, and so 
on, until self-consistency is reached. See Bransden and Joachain (2003) in 
further reading for a detailed discussion on the self-consistent field (SCF) 
problem. 

It is standard to use the hydrogen atom orbital labels, n, | and m, to 
label the orbitals of other atoms. The angular behaviour, given by ! and m, 
is indeed the same, but this is not so for the radial functions. In this case, 
orbitals designated by a certain n, | and m do not have the same radial 
structure as those given in Fig. 3.1 for hydrogen-like atoms. However they 
do have the same number of nodes. 

Within the orbital approximation presented above, the total wavefunc- 
tion for the system would be written 


Wty, ro, +--+ EN) = Pi(ty) G2(12), «++, ONGEN). (4.6) 


However this expression ignores the fact that one cannot distinguish 
between electron i and electron j. Before considering the indistinguisha- 
bility of electrons, it is necessary to explicitly consider the role of electron 
spin. 

To consider spin it is necessary to generalise the definition of an orbital 
to that of a spin-orbital: 


$i(j) = biltj, o), (4.7) 


where 6;(j) means that spin-orbital i is a function of the space-spin coordi- 
nates of electron j. In Eq. (4.7), c; is a spin coordinate variable; it will not be 
necessary to explicitly define these. The four dimensions of the space-spin 
variable of a single particle are similar to the four vectors used to represent 
relativistic behaviour of a single body in classical systems. 
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The space-spin wavefunction for the total system can now be written: 
Y(1,2, ..., N) = B)(1)@2(2)---By(N). (4.8) 


For non-relativistic treatments, such as that given by Eq. (4.1), the energy 
given by this wavefunction is the same as that given by w in Eq. (4.6). 


43 Indistinguishable Particles 


Consider a system with two identical particles. These particles can be any 
microscopic particles such as electrons, protons, neutrons, and so forth. If 
the wavefunction of these particles is V(1, 2) and the particles are indistin- 
guishable, then what property must this wavefunction have? 

Of course it is not the wavefunction but the probability distribution, 
P], which is physically observable. If the particles are truly indistinguish- 
able, this distribution cannot be altered by interchanging the particles. This 
means that 


24, 2)? = NOTIS (4.9) 


i.e. the probability distribution is unaltered by the interchange of particle 
1 and particle 2. Equation (4.9) has two possible solutions. There is the 
symmetric solution 


(1, 2) = HVQ, 1) (4.10) 

or the antisymmetric solution 
"(L2 = —¥(2, 1). (4.11) 
To explain the observed electronic structure of atoms, the German 


physicist Wolfgang Pauli (1900-1958) postulated that: 


Wavefunctions are antisymmetric with respect to interchange of 
identical Fermions. 


Fermions are any particles with half-integer spin such as electrons, protons 
or neutrons. This statement is known as the Pauli Principle and should 
not be confused with the less general Pauli exclusion principle discussed 
below. The Pauli Principle means that for any many-electron system the 
wavefunction must satisfy Eq. (4.11) for each pair of electrons. 
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Within the orbital approximation, a two-electron wavefunction which 
obeys the Pauli Principle can be written 


(1, 2) = = [@,(1)®,(2) — 6,(2)5,(1)] = —¥(2, 1), (4.12) 
V2 
where the factor of is included to keep the wavefunction normalised. 
The Pauli exclusion principle arises naturally from this expression for 
a two-electron wavefunction. If the two spin-orbitals, ®, and ®,, are the 
same, i.e. D, = ®,, then the total wavefunction, (1, 2), is zero. This solu- 
tion is not allowed as it cannot be normalised. Hence solutions which 
have the two particles occupying the same spin-orbital are inadmissible 
or excluded. The Pauli exclusion principle is often summarised as: 


No two electrons can occupy the same spin-orbital. 


This exclusion is of course the key to atomic structure and accounts 
naturally for the shell structures of atoms, and indeed nuclei. It also pro- 
vides the degeneracy pressure which holds up the gravitational collapse 
of white dwarfs and neutron stars. 


44 Electron Configurations 


For a hydrogen-like atom, the energy of the individual orbitals is deter- 
mined only by principal quantum number n. The energy ordering is 


E(1s) < E(2s) = E(2p) < E(3s) = E(3p) = E(3d) < E(49) --- 


For complex atoms the situation is not so simple. The degeneracy on 
the orbital angular momentum quantum number] is lifted. This is because 
electrons in low l orbits ‘penetrate’, i.e. get inside orbitals with lower 
n-values. Penetration by the low l electrons means that they spend some of 
their time nearer the nucleus experiencing an enhanced Coulomb attrac- 
tion. This lowers their energy relative to higher l orbitals which penetrate 
less or not at all. Figure 3.1 illustrates the greater penetration of s orbitals. 
As a result of this effect, the orbitals of complex atoms follow a revised 
energy ordering: 


E(1s) < E(2s) < EQp) < E(3s) < E(3p) < E(3d) ~ E(49) --- 


Following the Pauli exclusion principle, each orbital labelled nl actu- 
ally consists of orbitals with 21 + 1 different m values, each with two possi- 
ble values of sz. Thus the nl orbital can hold a maximum 2(21 + 1) electrons. 
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As an example, a p orbital, for which I = 1, can hold up to six electrons, 
since the following combinations of magnetic and spin quantum numbers 
are possible: 


1 1 1 
(ns - 5); (m=+1,8.=-3); (m=0,s.= +5); 
1 1 1 
(m=0,s.=-3); (m=-1,8.=+5); (m=-1,8.=-5). 


An atomic configuration is given by distributing electrons amongst 
the orbitals. The lowest energy or ground state configuration involves fill- 
ing the atomic orbitals in energy order from the lowest energy orbitals 
upwards. 

For example, carbon has six electrons which give the following ground 
state configuration: 


1s?2s?2p? ; 


where the superscript on the orbital gives the number of electrons in that 
orbital. No superscript is often used to denote a single electron. In car- 
bon, the 2p orbital contains only two electrons, which means that it is not 
full. Partially-filled shells usually give rise to several states with the same 
distribution of electrons between orbitals or configuration. This ‘open 
shell’ structure introduces complications which will be discussed further 
in Sec. 4.10. 

The neon atom has ten electrons and its ground state configuration is 


15:2s2p* ; 


All the occupied orbitals are full and this is known as a closed shell con- 
figuration. A closed shell or sub-shell (such as 2p°) makes no contribution 
to the total orbital or spin angular momentum ie. L or S. This property 
greatly simplifies the calculation of total angular momenta for complex 
atoms. 

Atomic ions which have the same number of electrons form what are 
called isoelectronic series. These ions have the same ground state con- 
figuration. Thus, for example, all ions with ten electrons are described as 
neon-like. Hence Fe XVII (or Fe!**) is called neon-like iron. 
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4.5 The Periodic Table 


The periodic structure of the elements was originally proposed by the 
Russian chemist Dmitry Ivanovich Mendeleyev (1834-1907) on the basis 
of careful study of the chemical properties of each element. The structure 
of the periodic table can be understood in terms of the configurations of 
the individual atoms. 

Table 4.1 gives the configurations of the first twenty elements in the 
periodic table. In this table, closed shells have be designated using the 
standard shell notation which labels the n = 1 orbitals as K, n = 2 as L, 
n —3as M, n = 4 as N, and so forth. Note that for potassium and calcium 
the M shell is not full since the 3d orbital is empty. The elements immedi- 
ately after calcium, known as the transition metals, begin to progressively 


Table 4.1. Atomic configurations of the first 20 elements 
in the periodic table. Z is the atomic number which corre- 
sponds to the charge on the nucleus. 


Atom Z Configuration 
hydrogen H 1 1s 
helium He 2 is^ 
lithium Li 3 K 2s 
beryllium Be 4 K 2g2 
boron B 5 K 2s2p 
carbon C 6 K  2s?2p? 
nitrogen N 7 K 2s2p? 
oxygen O 8 K 2s2p! 
fluorine E 9 K 2s2p 
neon Ne 10 K dep 
sodium Na 11 K L 3s 
magnesium Mg 12 K L 8s? 
aluminium Al 13 K L 3s?3p 
silicon Si 14 K L 8s?3p? 
phosphorus P 15 K L 3s?3p? 
sulphur S 16 K L 3s?3p! 
chlorine Cl 17 K L 8s?3p? 
argon Ar 18 K L 3s?3p$ 
potassium K 19 K L 3s?3pÓ 4s 
calcium Ca 20 K L 3s?3p$ 4s? 
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fill up the 3d orbital. These elements, which include the astronomically 
important iron, have particularly complicated electronic structures and 
associated optical properties. 

Atoms with the same electron configuration outside a closed shell 
share similar chemical and optical properties. Indeed it was the observa- 
tion of patterns in chemical behaviour which led to Mendeleev's original 
proposal of the periodic structure of the elements in 1871. 

For example, lithium, sodium and potassium all have a single elec- 
tron outside a closed shell. These atoms, known as the alkali metals, all 
have very similar optical properties. The spectra of alkali metals will be 
discussed in detail in Chapter 6. 

Electronically-excited states of atoms usually arise when one of the 
outermost electrons jumps to a higher orbital. These excited states can be 
written as configurations. Such excited states for helium might include 
1s2s, 1s2p or 1s3s. It should be noted that each of the configurations actu- 
ally gives rise to more than one excited state. 

States with two electrons simultaneously excited are possible but are 
less important. For many systems, all of these states are unstable. They 
have sufficient energy to autoionise by spontaneously ejecting an electron. 
For example, the lowest two-electron excited state of helium has the con- 
figuration 2s?. This state spontaneously decays to the 1s ground state of 
the He* ion and a free electron. 


4.6 Ions 


A neutral atom has the same number of electrons as the charge on the 
nucleus in atomic units. If an atom is not charge neutral it is called an ion. 
Under these circumstances the number of electrons, N, is not equal to the 
nuclear charge, Z, and the system carries a net charge of Z — N. 

Positive ions have N « Z. Ions with the same number of electrons, N, 
and different nuclear charge, Z, are called isoelectronic (see Fig. 4.1). Iso- 
electronic ions have similar structure and spectra, and are often referred 
to using the neutral atom at the head of the series, such as the example 
neon-like iron for Fe!®* discussed above. The spectrum of members of the 
sequence is modified by the different effective charge Ze. For example, a 
single outer electron may feel an effective charge of Zef = Z — N + 1. The 
energies of individual levels, and hence the transition frequencies between 
the levels, shift in proportion to Ze 
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Fig. 4.1. Isoelectronic series for low Z atoms. 


In the one-electron case, the ‘hydrogen-like’ ion transition frequencies 
scale closely with Z?. For complex ions there are considerable differences 
in detail between spectra of ions belonging to the same isoelectronic 
sequence. 

Negative ions have N > Z, although in practice only N = Z + 1 is 
important. Even then not all atoms can bind an extra electron to form a 
stable negative ion. Hydrogen, carbon and oxygen can bind an electron 
while helium and nitrogen cannot. Furthermore, most negative ions have 
only one stable level, and so possess no line ('bound-bound") spectrum. 
This means that the only possible transitions are continous bound-free 
absorption, also known as photoionisation. 

In general, negative ions are much less astronomically important than 
positive ones. However the hydrogen anion, H7, is a major source of opac- 
ity in cool stars including our Sun. The second electron is only bound by 
0.75 eV, in contrast to the first electron which is bound by 13.6 eV. It can 
therefore be photoionised by all light with wavelengths less than 1.6 jum, 
which lies in the infrared. 


47 Angular Momentum in Complex Atoms 


Complex atoms contain more than one electron and thus have several 
sources of angular momentum. In particular, the ith electron has orbital 
angular momentum /;, and spin angular momentum s;, which equals 
1. There is only one conserved angular momentum in the atom. Ignor- 
ing nuclear spin effects, this angular momentum is the total (spin plus 
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orbital) angular momentum for all electrons, J. Note that above and else- 
where the convention is followed that single electron angular momenta 
are represented by lower case letters, I, s, j, etc., and many electron angu- 
lar momenta are represented by upper case letters, L, S, J, etc. 

There are two coupling schemes or ways of summing the individual 
electron angular momenta to give the total angular momentum. 


47.1 L-S or Russell-Saunders coupling 


In Russell-Saunders coupling, the orbital and spin angular momenta 
of the electrons are added separately to give the total orbital angular 
momentum, L, 


L-Yl, (4.13) 
and the total electron spin angular momentum, S, 


S= > 3 (4.14) 


These are then added to give J 


J=Lt+S. (4.15) 
It is useful to remember that, as a result of the Pauli Principle, closed shells 
and sub-shells, such as ls? or 2p°, have both L = 0 and S = 0. This means 
that it is only necessary to consider ‘active’ electrons, those in open or 
partially-filled shells. In most cases this means only one or two electrons. 
When more than two angular momenta need to be added together, they 
should be added in pairs. The result is independent of the order in which 
the addition is performed. 


Worked Example: Consider O III with the configuration: Is?2s?2p3d. 
1s? and 2s? are closed, so contribute no angular momentum. 

For the 2p electron /; = 1 and sı = $i 
for the 3d electron l = 2 and s» = i. 
L=LĠ +hĠ > L=1,2,3; 

S= s +s, > 5—0,1. 
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Combining these using all possible combinations of L and 5, and the rules 
of vector addition, gives: 


J Level 

1 Tp? 

01,2 Spp DUIS 
2 

1 

3 


J=L+S5> 


1D9 

, 2, 3 AD, "Dope 
AES 

3 1 2,354 3P9,3F9,3F9, 


OU l2 l2 o3 n [m 


Each state of an atom or ion is characterised by a unique combination 
of L, S and J, known as a ‘level’, the notation for which is explained in 
Sec. 4.8. Thus twelve levels arise from the configuration 1s?2s?2p3d. Note 
that although some values of J appear several times, they all correspond 
to distinct states of the ion and it is important to retain them all. 


4.7.2 j-j coupling 


An alternative scheme for coupling angular momenta is to consider the 
total angular momentum, jj, for the ith electron by combining l; and s;: 


j-hts, (4.16) 
and then coupling these j’s together to give the total angular momentum. 


]J-Yji. (4.17) 


This scheme is known as j-j coupling. Again J = 0 for closed shells and 
sub-shells. 


Worked Example: Again consider 2 III with configuration pu ae 
For the 2p electron /; = 1 and s; = 3, j, = 4h + s, giving ji = 3,3; 
for the 3d electron lz = 2 and s2 = 5, j, = h + $9, giving j = 3,3. 


Combining these gives: 


A hR J 
J=i,ti> 3 f 12 

3 3 0,123 

i 3 2,8 

3 3 1,2,3,4 


62 Astronomical Spectroscopy 


These correspond to exactly the same twelve J values, or levels, obtained 
with L-S coupling. 


4.7.5 Why two coupling schemes? 


Given that L—S and j-j coupling schemes give the same results for J, why 
is it necessary to have two different schemes? The answer is that the two 
methods are used under different circumstances. 

In the non-relativistic formulation of atomic problem, as represented 
by the Schródinger equation (4.1), all states with same values for L and 
S have the same energy. In practice, relativistic effects split this degener- 
acy and these effects are usually introduced as a perturbation of the non- 
relativistic treatment. 

For light atoms, which generally mean atoms lighter than iron, rela- 
tivistic effects are weak. Under these circumstances the values of L and S 
are approximately conserved quantities, and the L-5 coupling scheme is 
the most appropriate. 

For heavy atoms, those beyond iron for example, relativistic effects 
are much stronger. Under these circumstances, L and 5 are no longer con- 
served quantities and j-j coupling is more appropriate. 

Light atoms are much more prevalent in astronomical spectra, so 
only L-S coupling is considered in any detail in this book. However it 
is important to remember that L and S are only approximately conserved 
quantities. 


48 Spectroscopic Notation 


Atoms with many active electrons can have a number of different energy 
levels arising from a single configuration. The splitting of configurations 
according to their L, 5 and J values requires a new method of denoting 
the different states. The standard notation is called spectroscopic notation 
and works within L-S coupling. Using this notation, states can be labelled 
as either ‘terms’ or ‘levels’. 

A ‘term’ is a state of a configuration with a specific value of L and S. It 
is denoted 


25-1r (0) y 


The leading superscript gives the spin multiplicity, i.e. the degeneracy of 
the spin state. A state with S = 0 is a ‘singlet’ as 25 +1 = 1; a state with 
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S= i is a ‘doublet’; one with 5 = 1 is a ‘triplet’, and so forth. The value of 
the orbital angular momentum, L, is given using a capital letter using the 
standard letter designation (see Table 3.1). Thus S represents a state with 
L —0, P a state with L = 1, D means L = 2, and so forth. Historically, this 
notation arises from analysis of the sodium spectrum (see Sec. 6.1). 

The trailing superscripted 'o' in the term means 'odd parity'. Even 
parity terms are usually written without a superscript. The parity of a par- 
ticular term, which is the same for all terms arising from a particular con- 
figuration, is defined in the following section. 


Worked Example: The O ii ion with configuration 1s?2s?2p3d can have 
L — 1,2 or 3 and S = 0 or 1 (see Sec. 47.1). Taking all combinations of L 
and S, this gives rise to terms 


lpo, pe, Ip, Spe lpo and 3pe . 


It should be noted that the splitting of configurations into terms with 
different energies arises even in the non-relativistic (Schródinger) formula- 
tion. However, inclusion of relativistic effects splits these terms into levels 
according to their J value. A level is denoted 


2541 Le 


where the only difference from the term notation is the addition of the 
subscript, J, which represents the total electron angular momentum. 

Thus, for example, a 3F° term, such as the one arising from the 
1s?2s?2p3d configuration of O I1, splits into three levels with J = 2,3 and 
4. These levels are designated ?F$, ?F? and 3F9. 

For completeness it is worth noting that for each level, characterised 
by a particular value of J, there are 2J + 1 sub-levels. These are called 
states and are designated by the total magnetic quantum number 
My which takes values 


Myj--J,-Jc*M...,] -MJ. 


These states are degenerate in the absence of an external field. The splitting 
of levels into states in a magnetic field is generally known as the Zeeman 


effect. 
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49 Parity of the Wavefunction 


The parity of the wavefunction is determined by how the wavefunc- 
tion behaves upon inversion. Inversion is the operation of reflecting the 
wavefunction through the origin, here the atomic nucleus, and is equiv- 
alent to replacing vector r with —r. Given the symmetry of the atom, the 
square of the wavefunction, i.e. the probability distribution of the elec- 
trons, must be unchanged by this operation. Neglecting spin, this means 


vn, Vor seer rn) = cwn, — Eye, T YN) 3 (4.18) 


Even parity states are given by +7 and odd parity states are given by —v. 

In practice the parity of all terms and levels arising from a particular 
configuration can be determined simply by summing the orbital angular 
momentum quantum numbers for each of the electrons. With this simple 
rule, the parity is given by 





(1), (4.19) 


As closed shells and sub-shells have an even number of electrons, it is 
again only necessary to explicitly consider the active electrons. 

Thus for the O III configuration Is?2s?2p3d, it is only necessary to con- 
sider the sum /(2p) = 1 and /(3d) = 2. This gives (—1)!*? = —1, which 
explains why all the terms and levels arising from this configuration were 
all labelled odd above. 

The parity of a configuration is important since it leads to a rigorous 
dipole selection rule known as the Laporte rule. The Laporte rule states: 


All electric dipole transitions connect states of opposite parity. 


In other words (strong) transitions can only link configurations with even 
to those with odd parity, and vice versa. 


4.10 Terms and Levels in Complex Atoms 

A single configuration can lead to several terms. These terms have differ- 
ent energies. It is worth considering a few examples. 

Example 1: The helium atom. 


(1) The ground state is 1s?. 
This is a closed shell, with L = 0 and S = 0, hence it gives rise to a 
single, even parity term !S, and level 'So. 
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(2) The first excited configuration is 1s2s. 
This has l, = l, = 0 and hence L = 0, 
but sı = s) = 1 giving both S = 0 (singlet) or S = 1 (triplet) states. The 
energy ordering of atomics states is given by Hund's rules. Hund's 
first rule governs ordering of terms with different spin multiplicities: 


For a given configuration, the state with the maximum spin multiplicity 
is lowest in energy. 


So the °S term (°S; level) is lower in energy than the !S term ('S level). 
In practice the splitting between these terms is 0.80 eV. 


(3) The next excited configuration is 1s2p, which has odd parity. 
This has |; = 0 and l = 1, giving L = 1; 
again $1 = $5 = L giving both S = 0 and S = 1 terms. 
Following the rule above, the 3P? term is lower than the !P? term, 
in this case by 0.25 eV. The ?P? is also split into three levels: ?P9, ?P? 
and ?P9. 
Figure 5.2 depicts the energy levels of helium in a form known as a Gro- 
trian diagram. The layout of these diagrams is discussed in Sec. 5.4. 


Example 2: The carbon atom. 


Start by considering the excited state configuration 1s?2s?2p3p. 
It is only necessary to consider the outer two electrons for which: 


h = 1, $1 = 

h = 1, $5 — 
These give L = 0, 1, 2 and S = 0,1, which give rise to the following terms, 
all with even parity: 1S, 3S,!P, ?P, !D and ?D. 


Now consider the ground state configuration of carbon Is?2s?2p?. 

This configuration also has J, = 1, s1 = i and h = 1, s2 = i. 

However the Pauli Principle restricts which terms are allowed. For exam- 
ple the term ?D includes the state: 


(L =1,m, = +41, sı —1, sız = 41) 
(Io = 1, m» = +1, $2 = 1, S2; = +3) 


which is allowed when 14 = 2, 1; = 3, but is forbidden when nı = 2,n2 = 2 
by the Pauli exclusion principle since both electrons have precisely the 
same quantum numbers. 
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There are general methods of determining which terms are allowed 
for a configuration with a multiply occupied open shell [see Bransden and 
Joachain (2003) in further reading]. However there is a rule of thumb which 
suffices for present purposes. It turns out that for systems with equivalent 
electrons, that is, electrons which have the same n and ! values, then the 
sum L + S for these electrons must be even for the Pauli Principle to be 
satisfied. 

The ground state configuration of carbon, C I, thus gives terms lg 3P 
and !D. The?P term has the highest spin and is thus the ground state term. 
The other two terms have however, the same spin multiplicity, so which is 
lower in energy? Hund's second rule states: 


For a given configuration and spin multiplicity, the state with the maximum 
orbital angular momentum is the lowest in energy. 


In the case of the ground state configuration of carbon, the ID state 
lies 1.42 eV lower in energy than the !S state, but is 1.26 eV above the 
?P state. 

The examples above have only considered terms and have neglected 
splitting according to J where it arises. This fine-structure splitting 
becomes increasingly important for high Z (i.e. heavy) atoms. 

The ?P ground state of carbon has L = 1 and S = 1, which lead to 
J — 0, 1, 2. Allowed levels are thus ?Po, ?P, and ?P;. 

The energy order of these is given by Hund's third rule: 


The lowest energy is obtained for lowest value of J in the normal case and 
for highest J value in the inverted case. 


The normal case is a shell which is less than half filled, for example 2p? as 
in carbon. The inverted case is a shell which is more than half full such as 
the 2p* found in the ground state of atomic oxygen, which also has a ?P 
ground state. 

Thus for carbon one gets the energy order: 


3Po < 3P} < 3P,, 
whereas for oxygen, one gets 
3P, < 3P} < ?Po p 


It should be noted that this situation does not arise for configurations 
with half-filled shells since the lowest energy term, the only one for which 
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Hund's rules apply rigorously, always has L = 0 so that there is only a 
single level. 
To summarise the energy ordering of levels given by Hund's rules: 


(1) For a given configuration, the term with maximum spin multiplicity 
lies lowest in energy; 

(2) For a given configuration and spin multiplicity, the term with the 
largest value of L lies lowest in energy; 

(3) For atoms with less than half-filled shells, the level with the lowest 
value of J lies lowest in energy; 

(4) For atoms with more than half-filled shells, the level with the highest 
value of J lies lowest in energy. 


It should be noted that Hund's rules are only applicable within L-S cou- 
pling. Furthermore, they are only rigorous for ground states. However 
they are in practice almost always followed for all states of atoms and 
ions. They are therefore also useful for determining the energy ordering of 
excited states. 


Problems 


4.1 Whatis the ground state configuration, term and level of the beryllium 
atom, Be? One of the outer electrons in Be is promoted to the 3d orbital. 
What terms and levels can this configuration have? 

4.2 An excited helium atom has the configuration 3d?. What values of the 
total orbital angular momentum quantum number, L, and total spin 
angular momentum quantum number, 5, are allowed? Use spectro- 
scopic notation to give the terms which arise from the combinations 
of L and S allowed under the Pauli exclusion principle. Use Hund's 
rules to suggest the energy order in which these terms are likely to 
occur. For each term deduce the allowed values of the total angular 
momentum quantum number, J. Hence give the full designation of 
each level, including parity, using spectroscopic notation. Which level 
has the lowest energy? 

4.3 The symbol for a particular level is quoted as ‘F3. What are the values 
of L, S and J for this level? How many states does it have? What are 
the other levels for this term? 
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Use the value of S to determine the minimum number of electrons that 
could give rise to this term. Suggest a possible configuration that could 
give this term. 

4.4 Symbols for particular levels of three different atoms are written as 
1D,,°D 3 and Ap 3- Explain in each case why the symbol must be wrong. 

4.5 The lithium atom Li has three electrons. For each of the following 

configurations, what terms will be present: (a) 1s?2p, (b) 1s2s3s, and 
(c) 1s2p3p? 

4.6 The 1s*2s*2p°3s3p°3d5p configuration of Ca is formed during recom- 
bination. Derive the terms that arise from this configuration and 
explain which you would expect to be lowest in energy. What levels 
can arise from the lowest energy term? 


CHAPTER FIVE 


HELIUM SPECTRA 


‘Born of the Sun.’ 


— Stephen Spender, 7 Think Continually (1933) 


5.1 Herand Hell Spectra 


Helium is the second-most abundant element after hydrogen compris- 
ing more than 25% of the Universe’s atomic matter by weight. Helium 
can exist in its atomic form, Hel, as singly ionised Hell, or as doubly 
ionised He 1. He III is of course the base helium nucleus He^*, which is 
also the a particle formed during radioactive decay, and which has no 
spectrum. 

The He* ion is a one-electron system. It therefore has a hydrogen-like 
spectrum except that the binding of the energy levels and the transition 
frequencies are scaled by a factor of Z?, where for helium, Z = 2. 

Thus for example, the n = 2-1 Lya for hydrogen is observed at 
a wavelength of 1216 À. The corresponding transition in Hell, which 
is also referred to as Lya, is observed at 304À. This transition is 
observable in absorption in the interstellar medium (ISM) where it is 
difficult to measure helium abundances by other means. Such obser- 
vations have been made from the extreme ultraviolet explorer (EUVE) 
satellite. 

As discussed in Chapter 3, hydrogen lines occur at 


m E) - 
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For the Balmer series, nı = 2, n; = 3,4,5,.... For He Il, similar lines in the 
visible occur at 


1 1 1 1 1 
x = Rue (n) fe —— 31 (5.2) 
i a E s 


where Rye is the Rydberg constant for Het, which is intermediate in value 
between Ry and R» (see Sec. 3.7). Lines with nı = 4 and nz = 6,8,10,... 
are therefore almost coincident with the Balmer series. The shift due to the 
reduced mass factor in the Rydberg constant can only be seen at high reso- 
lution. However, intermediate lines with 11; = 5,7,9, ... are also observed. 
It was originally thought that this spectrum, which is known as the 
Pickering series, was due to atomic hydrogen with half-integer quantum 
numbers! The Pickering series can be observed in the spectra of O-stars. 
Indeed the classification of O-stars depends on the relative strength of the 
absorption by He II and Her in their atmosphere (see Fig. 5.1). 

Emission lines of Hell are also observed among the recombination 
line series in a nebulae. However these emissions only come from the 
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Fig. 5.1. Spectra of early (O4 star HDE 269698), mid (O7 star AzV 232) and late 
(O9 star Sk-66° 169) type O-supergiant stars showing the sensitivity of the subtypes 
to the ratio of He 1 4471 A to He 11 4542 A absorptions. The spectra were recorded 
using the Very Large Telescope (VLT) by P.A. Crowther et al. [see Astrophys. J. 579, 
774 (2002)]. (P.A. Crowther, private communication.) 
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region of nebula close to the central star where there are sufficient high- 
energy, ultraviolet photons to ionise He*. He* requires 54.4 eV, compared 
to 13.6eV which is needed to ionise H. HeI recombination spectra also 
come from a smaller physical region than the H II region as the ionisation 
potential of He is 24.6 eV. 

Strong He I emissions are observed from nebulae including the lines: 


4471 A 1s2p - 1s4d ?P? -?D, 
5875 A 1s2p - 1s3d ?P? - °D, 
6678 A 1s2p - 1s3d !P? - !D. 


The 5875À line lies in the yellow and absorbs in the Sun. Observa- 
tion of this previously unknown line in the Sun during the solar eclipse 
of 1868 led Sir James Lockyer (1836-1920) to postulate the presence of 
a new element: helium. He was widely ridiculed for his boldness but 
was vindicated when Sir William Ramsay (1852-1916) isolated helium 
in his laboratory at University College London more than twenty years 
later. 

Like hydrogen, transitions involving the ground state of He lie far 
in the ultraviolet. For example, the He 1 1s?- 1s2p 1S—1p° transition lies 
at 584 À. This transition is known as the He! resonance line. The phrase 
'resonance line' is used to denote the longest wavelength, dipole-allowed 
transition arising from the ground state of a particular atom or ion. 


5.2 Selection Rules for Complex Atoms 


Strong transitions are driven by electric dipoles. Electric dipole selection 
rules are of two types: rigorous rules which must always be obeyed, 
and propensity rules which, when violated, lead to weaker (forbidden) 
transitions. 

The rigorous selection rules which all electric dipole transitions, even 
so-called forbidden ones, must satisfy are: 


(1) AJ must be 0 or +1 with J = 0 — 0 forbidden. 
(2) AM; =0, +1. 
(3) Parity changes i.e. even — odd. This is the Laporte rule (Sec. 4.9). 








All electric dipole transitions must obey these rules. There is an additional, 
and more complex set of rules, which is satisfied by systems with a single 
electron but not rigorously by complex atoms. Failure to satisfy these rules, 
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which I will call propensity rules, does not completely rule out a transition 
but generally makes it much weaker. 

The propensity rules for electric dipole transitions, which lead to 
stronger transitions if satisfied, are: 


(4) The spin multiplicity is unchanged, AS — 0. 

(5) Only one electron jumps, i.e. the configuration of the two states must 
differ by only the movement of a single electron. This movement is 
governed by the rules An any and Al = +1, thus: 





2s? & 2s2p is allowed; 
2s? — 283d or 2s? +> 3p? is forbidden. 


Configuration interaction (CI) weakens this rule since CI means that 
the real wavefunction is a (small) mixture of different configurations. 
For example, the ground state of Be, whose configuration is usually 
written 1s?2s?, is more precisely represented by 95% 1s?2s? mixed with 
about 5% of the same symmetry configuration 1s?2p?. 

(6) AL —0, +1, L = 0 o 0 forbidden, thus: 





15 1p? or ĉD © ?P? are allowed; 


16 «415° or 9S à ?D? are forbidden. 


Note that it is possible to change / and not L, for example, 2p3p 1P — 
2p3d !P° is allowed. Such a transition is observed for CI. 


A full list of selection rules are given in Table 5.1. 

Electric dipole transitions which satisfy all the rigorous selection rules 
as well as the propensity rules are referred to as allowed transitions. These 
transitions are generally strong and have Einstein A coefficients which are 
typically bigger than 106 s 1. 

Photons do not change spin, so transitions usually occur between 
terms with the same spin state, as expressed by the rule AS — 0. However, 
relativistic effects mix spin states, particularly for high Z atoms or ions. 
As a result of relativistic effects one can get (weak) spin changing transi- 
tions; these are called intercombination lines. Intercombination lines are 
denoted by one square bracket, for example: 


Cir] 2s? 1S - 2s2p ?P? at 1908.7 À . 


This transition is important because the C?* 2s2p ?P? state is metastable, 
i.e. it has no allowed radiative decay so that this transition determines the 
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Table 5.1. Selection rules for atomic spectra. Rules 1, 2 and 3 must always 
be obeyed. For electric dipole transitions, intercombination lines violate 
rule 4 and forbidden lines violate rule 5 and/or 6. Electric quadrupole and 
magnetic dipole transitions are also described as forbidden. 



























































Electric dipole Electric quadrupole Magnetic dipole 
1. AJ=0,+1 AJ —0, «1, : AJ 20, +1 
Not J 20—0 Not J 20— 0,5 -4,0—1 Not J 20—0 
2. AM;=0,+1 AM, = 0, AG 2 AM, —0, +1 
3. Parity changes Parity unchanged Parity unchanged 
4. AS=0 AS =0 AS =0 
5. Oneelectronjumps One or no electron jumps No electron jumps 
An any An any An=0 
Al = +1 Al=0, +2 Al = 0 
6. AL=0,+1 AL=0,+1,+2 AL — 0 
Not L=0-—0 Not L=0-—0,0-1 


lifetime of this state. Actually, the situation is more subtle than this. The ?P? 
term splits into three levels: ?P9, ?P? and °P9. The electric dipole intercom- 
bination line at 1908.7 A is actually ‘Sp -?P?. It has an A value of 114 s 1. 

The transition 1S) -?P9, which occurs at 1906.7 A, is completely for- 
bidden by dipole selection rules as AJ — 2. It only occurs via a very weak 
magnet quadrupole transition. The 1906.7 À line is 10? times weaker than 
the already-weak line at 1908.7 A; it has an A value of 0.0052 s^!. These 
two lines can be used to give information on the electron density, as dis- 
cussed in Sec. 7.1. Finally the transition !S; -?P9 is a J = 0-0 transition, 
which is completely forbidden by both dipole and quadrupole selection 
rules. This transition is not observed. 

Electric dipole transitions which violate the propensity rules 5 and/or 
6 are called forbidden transitions. These are labelled by square brackets. 
For example, 


1906.7 A [C 111] 2s? 'Sp—2s2p ?P5; 
322.57 A [C ni] 2s? !S; -2p3s ! P5; 


are both forbidden lines of C?*. The former is a magnetic transition while 
the latter is an electric dipole transition involving the movement of two 
electrons. Forbidden transitions are generally weaker than intercombina- 
tion lines. 
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It is also possible to get transitions driven by higher electric multi- 
poles or magnetic moments. The only important ones of these are electric 
quadrupole and magnetic dipole transitions. The selection rules for these 
transitions are also given in Table 5.1. Even when all the rules are satisfied, 
electric quadrupole and magnetic dipole transitions are both much weaker 
than the allowed electric dipole transitions. They are thus also referred to 
as forbidden transitions. 

Typical lifetimes, that is inverse Einstein A coefficients, for allowed 
decays via each mechanism are 


-8 -3 
Tdipole ^" 10 "s, Tmagnetic ~ 10s, Tquadrupole ^" 1s. 


These timescales mean that states only decay by forbidden transitions 
when there are no decay routes via allowed transitions. 

Finally it should be noted that even the rigorous selection rules given 
above can be modified when nuclear spin effects are taken into considera- 
tion. These result in rigorous selection rules for electric dipole transitions 
based on the final angular momentum. In particular: 





AF must be 0 or € 1 with F = 0 e 0 forbidden. 


It is only very rarely necessary to consider this. 


5.3 Observing Forbidden Lines 


States decaying only via forbidden lines live for a long time on an 
atomic, if not an astronomical, timescale. Such states are called metastable 
states. 

Forbidden lines are often difficult to study in the laboratory as 
collision-free conditions are needed to observe metastable states. In this 
context it must be remembered that laboratory ultrahigh vacuums are sig- 
nificantly denser than so-called dense interstellar molecular clouds. Astro- 
physically, low density environments are common. In these environments 
the time between collisions is very long and an atom in an excited state 
has time to radiate even when it is metastable. 

Emissions due to forbidden lines are important in hot, low density 
regions such as HII regions, planetary nebulae, the solar corona and the 
Earth’s aurora. Observing them gives direct information on the popula- 
tions of excited levels. As the transitions are weak they have low optical 
depth and therefore give reliable population information. 
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Forbidden transitions are usually only important when an excited 
state cannot decay via an allowed transition, ie. when this state is 
metastable. For this reason, forbidden lines are normally only important 
for low-lying states since higher states nearly always have possible radia- 
tive decay routes via allowed transitions. For neutral atoms or ones with 
low ionisation, this means that forbidden transitions are often observed in 
the infrared. Of course, this wavelength shifts with the level of ionisation 
so that forbidden lines for ions occur throughout the visible, ultraviolet 
and even at X-ray wavelengths (see Sec. 5.6). 

One common source of forbidden transitions in the infrared arises 
from the relaxation of excited terms within the ground state configuration 
and levels within the ground state term. Such transitions are completely 
dipole-forbidden by the Laporte rule and hence undergo weak magnetic 
transitions instead. 


5.4 Grotrian Diagrams 


The German astrophysicist Walter Grotrian (1890-1954) invented a simple 
diagrammatic means of representing the many states and transitions of 
atoms and ions. These are called Grotrian diagrams and several of them 
can be found in this book. 

The standard structure of a Grotrian diagram is as follows: 


(1) The vertical scale is energy. It starts from the ground state at zero, 
and extends to the first ionisation limit. Sometime the binding energy, 
expressed relative to the first ionisation limit, is given as the right-hand 
vertical scale. Terms (levels) are represented by horizontal lines. 

(2) States with the same term, or sometimes the same level if fine structure 
effects are large, are stacked vertically and labelled by the principal 
quantum number n of the outer electron. For example, He 1 singlet 
states (see Fig. 5.2), with configuration 1snd appear in the column 
headed !D and are labelled by n = 3,4,5,.... 

(3) Terms are grouped by spin multiplicity. 

(4) States are linked by observed transitions with numbers giving the 
wavelength of the transition, usually as an integer, in À. Thicker lines 
denote stronger transitions and forbidden transitions are given by 
dashed lines. 


There are some variations on this structure of the diagram accord- 
ing to the composer and the system being considered. Grotrian diagrams 


76 Astronomical Spectroscopy 



































| 
F | F cm! 
| 
aot 4 > 18697 
| 
| 
| 
l -]- 50,000 
| 
| 
| 
z^ 
538.9 L 
ri n y 1 
| ++ 100,000 
592.4 l 
d l 
| 
| He I 
l Ionisation level 
! 198310.8 cm”! -]- 150,000 
| (24.5876 eV) 
| 
| 
I 
| 
ott | 
| 


— 200,000 


Fig. 5.2. Grotrian diagram for Het. The running numbers denote the principal 
quantum number of the active electron. The left-hand side of the figure is for ‘para’ 
singlet helium and the right-hand side is for ‘ortho’ triplet helium. 


are useful for getting a feel for a particular spectrum but not suitable for 
detailed or high accuracy work. A particularly comprehensive and useful 
series of Grotrian diagrams for each element and ionisation stage has been 
compiled by Bashkin and Stoner (see further reading). 

Figure 5.2 gives the Grotrian diagram for Hel. This diagram splits 
into two: the left-hand side gives singlet (‘Para’) He, while the right corre- 
sponds to triplet (‘Ortho’) He. The terms ortho and para, short for ‘ortho- 
dox’ and ‘paradox’, arise from the spin multiplicity of the atom with the 
ortho or orthodox states being three, 25 + 1 = 3, times as likely as the 
25 +1 = 1 para or paradox states. This terminology is not used for other 
atomic systems but has been adopted for molecular spectra where it is the 
nuclear rather than the electron spin multiplicities which give rise to dif- 
fering weights in the spectrum. 

The ortho and para states of Her are linked by very weak intercom- 
bination lines, two of which are shown in Fig. 5.2. In particular, the 2 35 
state can only decay to the 1 !S ground state. This transition is dipole-for- 
bidden by the Laporte rule. It decays via a magnetic dipole transition. The 
transition is weaker than an allowed magnetic dipole transition since it is 
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both spin-forbidden, as AS — 1, and involves an electron jump. It has A — 
1.2x10-4^s^! which means the helium 1s2s ?S state has a lifetime of 8000 s. 


5.5 Potential Felt by Electrons in Complex Atoms 


In the orbital approximation, each electron moves in an effective potential, 
Vi(ri) (see Sec. 4.2). V; is complicated but aspects of its general form can be 
understood by considering its behaviour as r; — 0 and as r; — oo. 
Consider an atom or ion with N electrons and nuclear charge -- Ze. 
For simplicity assume that N < Z. As the coordinate of the ith electron 
gets very large (r; — oo), there are N — 1 electrons left near the nucleus. In 
this regime the potential felt by this outer electron is given by: 
>% Z (N-1 Z-N+1 
o e ee ee (5.3) 
Yi Yi Yi 
For a neutral atom N = Z, and Vj(rj) ^ -ri. Conversely as the ith electron 
gets very close to the nucleus (r; — 0), the electron moves inside all the 
other electrons and therefore feels the full nuclear charge. 
s0 Z 
vir) == =. (5.4) 
i 
This analysis does not solve the problem of the form of V;(r;) but it 
does give its limits: 


- (5.5) 


-Z < V(r) < — (=) y 
r 

This is shown schematically in Fig. 5.3. It is tempting to interpolate 

between the two limits shown on the figure but this would be wrong. In 

practice, the shell structure of the atom introduces fluctuations into the 

potential felt by the individual electrons at intermediate values of r. 

If the outermost electron of an atom or ion lay entirely outside the 
other electrons it would be quasi-hydrogenic and its binding energy 
would be given by the expression 

2 


Z 
Ent ~ -Ro -5 , Zæ=Z-N+1. (5.6) 


However penetration makes the potential more attractive. As discussed 
in Sec. 4.4 and illustrated by Fig. 3.1, penetration is the ability of an outer 
electron to penetrate the charge cloud of the inner electrons. This can be 
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Fig. 5.3. Approximate shape of the potential felt by a single electron in a complex 
atom with N electrons. The vertical axis represents rV (r). 


thought of as raising the value of Zep felt by the outer electron, but it is bet- 
ter to think of penetration as lowering the principal quantum number n: 


eff f (5.7) 





where v is known as the effective quantum number. It is useful to define 
the effective quantum number in terms of n, and something called the 
quantum defect, Hı, via the expression 


Vy] =N— Hn. (5.8) 


As s orbitals penetrate more than p, p orbitals penetrate more than d, and 
so forth, one finds in general that 


Eno < Eg € Ego ..., 
which means that 
Hno > [nt > Hn2.... 


Quantum defects are discussed in more detail in Sec. 6.1. 


5.6 Emissions of Helium-Like Ions 


Any ion with only two electrons is known as a helium-like ion. It turns 
out that spectra of helium-like or K shell ions have become an important 
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diagnostic because their characteristic spectra can give important physical 
information about a hot environment. If a helium-like ion has one electron 
excited to its n = 2 level, then the various states, which are all fairly close 
in energy, can decay by the following routes: 


W 1s2p- 1s? 1p? —16,, resonance line, 
X 1s2p - 1s? ?P2 -!S,, magnetic quadrupole line, 
Y 1s2p- 1s? 2p? — 150, intercombination line, 


Z 1s2s—1s* 361 — 190, forbidden transition . 


The letters are the standard labels for these transitions, an example of 
which is given for O VII in Fig. 5.4; note that lines X and Y are close together 
and are usually not separately resolved. 

Careful monitoring of the relative intensities of these lines can give 
information both on the temperature and density of the environment and 
the ionisation mechanism involved. In particular the critical density n, (see 
Sec. 2.6) of the transitions increases strongly with the nuclear charge, Z. 
So, for example, n; of Si XIII is about 10000 times greater than n, of the 
isoelectronic C v. This means that a large range of electron densities can be 
studied by monitoring a variety of helium-like ions. 

Figure 5.5 gives a simple Grotrian diagram for O VII which is helium- 
like, and one of the species commonly observed. Its spectrum is also 
observed in solar flares. The O VII resonance line lies at 21.6 A, which is 
in the X-ray region of the spectrum. At these wavelengths the transitions 
depicted in the figure appear as a triplet as the intercombination doublet is 
usually not resolved. These transitions lie close enough together for their 
spectra to be recorded simultaneously by instruments on board satellites 
such as the XMM-Newton. X-ray spectra are discussed in Chapter 8 where 
Fig. 8.3 shows emission from He-like iron or Fe xxv. 
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Fig. 5.4. Helium-like oxygen, O VII, triplet X-ray spectra from the coronae of the 
stars Procyon (left) and Capella (right). See the text for an explanation of the labels 
for each line. Spectra were recorded using the satellite Chandra. [Adapted from 
S.M. Kahn et al., Philos. Trans. R. Soc. Lond., Ser A 360, 1923 (2002).] 









80 Astronomical Spectroscopy 
2lpo 
23p° 
1630À 
Intercombination 23s 
22.1À 
Resonance : 
line 21.6À _” Forbidden 
' — 228À 
1!S 

Fig. 5.5. Term diagram showing for helium-like oxygen, O VII, showing transi- 
tions from the 1s2/ states. 
Problems 
5.1 A transition of hydrogen-like helium, “Het, is observed close to 


5.2 


5.3 


5.4 


5.5 


hydrogen Ha. Between which states of He* is this transition? Estimate 
the wavenumber of this transition, stating any assumptions made. 
Emission of Hell are observed at 1640 A. To what transitions does 
this radiation correspond to? What other emission lines must also be 
present in the spectrum? 

A nebula is composed only of hydrogen and helium, and consists of 
four shells around a hot, central star. In shell A, closest to the star, 
all atoms are fully ionised; in shell B, helium and hydrogen are both 
singly ionised; shell C is an H II region with neutral helium; shell D, 
which is furthest from the central star, is an H I region. Explain what 
(if any) recombination line spectra you would expect to observe from 
each region. 

The helium-like ion N VI is formed in the hot intergalactic medium 
with configurations 1s2s and 1s2p. What levels can be formed from 
these configurations? By what mechanism would you expect each 
level to emit to the 1s? ground state? Order each of these transitions 
according to its approximate strength. 

The lithium atom, Li, has three electrons. Consider the following con- 
figurations of Li: (a) 1s?2p, (b) 1s2s3s, (c) 1s2p3p. By considering the 
configurations only, state which of the three sets of transitions between 
the configurations (a), (b) and (c) are electric dipole allowed and elec- 
tric dipole forbidden transitions? 


CHAPTER SIX 


ALKALI ATOMS 


Lithium, sodium, potassium and rubidium all have ground state electronic 
structures which consist of one electron in an s orbital outside a closed 
shell. This single ‘optically active’ electron gives these atoms, the alkali 
metals, similar chemical behaviour and fairly simple spectra. Even so, the 
presence of the inner or core electrons lead to a number of complications 
which are not present in the spectrum of simple one-electron atoms. 


6.1 Sodium 


Sodium, Na, has Z = 11 and a ground state configuration of 1s?2s?2p93s!. 
If the outer 3s electron was completely screened then it would feel an effec- 
tive nuclear charge Ze, = 11 — 10 = 1 and its energy levels would obey the 
Rydberg formula. In practice, the 3s electron penetrates and reduces the 
effective quantum number of the electron giving a revised formula: 
Zn 

d (n ES ma i p 
where un is the quantum defect. The quantum defect formula was origi- 
nally proposed by Rydberg. Unlike his hydrogenic formula, Eq. (3.9), this 
depends on / as the degree of electron penetration is |-dependent: in par- 
ticular, electrons with low / penetrate more and hence have lower energy. 
An advantage of this form is that usually j1,,; only depends weakly on n, 
allowing values of the quantum defect to be transfered between states with 
different n but the same 1. Table 6.1 gives values of the quantum defect for 
low-lying states of sodium. Note the u43 ~ 0 means that there is effectively 
no penetration for f electrons; this is also true for electrons with | greater 
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Table 6.1. Quantum defects 1; for sodium. 








I n=3 n=4 n=5 n=6 n=% 
0 s 1.373 1.357 1.353 1.351 1.348 
1 p 0.883 0.867 0.862 0.857 0.855 
2 d 0012 0013 0.014 0.014 0.015 
3 f — 0.000 0.000 0.000 0.000 
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Fig. 6.1. Grotrian diagram for sodium. Transitions are labelled with more than 
one wavelength due to the effects of fine structure. 


than 3. Use of the quantum defects for sodium gives a hydrogen-like level 
structure with splitting on / (and J) see Fig. 6.1. 


Worked Example: The spectrum of S VI, sodium-like sulphur, shows a 
series of ns—3p transitions with the series limit at 604310 cm™!. The first 
two transitions in this series lie at 257109 cm~! and 398238 cm~!. Estimate 
where the third transition in the series should lie. 
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Fig. 62. Simplified Grotrian diagram for S VI showing transitions discussed in 
the worked example. 


This data is best shown diagrammatically as can be seen from Fig. 6.2. 
Using this figure, the energy of the 4s level, E(4s) is 257109 — 604310 = 
—347201 cm! and E(5s) = 398238 — 604310 = —206072 cm~!. For the outer 
electron in S?*, Ze = 6, so the quantum defects are given by the expression 


Rx 


p(ns) =n—6 Ets): 





Using the values of E(4s) and E(5s) determined above gives u(4s) = 0.627 
and u(5s) = 0.622. Assuming u(6s) ~ 0.62 gives E(6s) = —136487cm~! and 
hence the 6s—3p transition at 467823cm™!. In practice, this transition is 
observed at 467949 cm '!. Note that E(ns) is negative since it is binding 
energy which goes to zero at the ionisation limit. 

The Grotrian diagram for sodium, Fig. 6.1, shows that transitions with 
Al = +1 and AL = +1 dominate the spectrum. However, as implied by the 
multiple wavelengths given for most transitions, these lines are actually 
split into components. 
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Table 6.2. Spectral series of sodium. 


Series name Transitions n values Multiplicity 
Sharp n ?$; —3?pP9, n=4,5,6,... doublets 
»2 
Principal n?po , 5 3?9, n = 3,4,5,... doublets 
23 2 
Diffuse n "D; 373 2Po, =n =3,4,5,... triplets 
: 23 


Fundamental n?p9 , > 32D n=4,5,6,... triplets 
2/3 


5 
2 


MIO 


" 


Consider the spectrum of sodium and, in particular, emissions into 
the n = 3 shell of sodium. The spectral series, most of whose transitions 
lie in the visible, are given in Table 6.2. These series are named after their 
spectral appearance rather than their discoverers. Note the initial letters of 
the names of each series; these gave rise to the use of the notations s, p, 
d and f for orbital angular momentum states / = 0,1,2 and 3 respectively, 
which of course, are the corresponding orbital angular momenta of the 
emitting states. Orbitals with / > 3 are simply labelled alphabetically, with 
the omission of the letter j. 

As noted in Table 6.2, each member of each series is split into more 
than one line. These splittings are large enough to be resolved with a high- 
resolution spectrograph and, as discussed later in this chapter, are very 
important observationally. However before discussing their astronomical 
importance it is necessary to understand the physics that gives rise to the 
different multiplets, and in particular why some series give two lines (dou- 
blets) and others three lines (triplets). 


6.2 Spin-Orbit Interactions 


Both the orbital angular momentum L and spin angular momentum S 
give rise to an internal magnetic field within the atom. It is the magnetic 
interactions between the spin and orbital motions of the electrons which 
split terms into levels. 

The most important interaction is between the spin magnetic moment 
of an electron which is given by 


p= —2778, (6.2) 


and the magnetic field due to the orbital motion of the electron. This is the 
only source of spin-orbit interaction considered here. 
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Equation (6.2) defines the spin magnetic moment in terms of the Bohr 
magneton ug, where 


eh 1 
2m, 2 





a.u. = 9.27 x 10 ^ JT *. 








HB 


The Bohr magneton is the basic unit, or quantum, of magnetic moment. 
The magnetic field due to the orbital motion of an electron is 
WoxrdVv 
B= y 6.3 
= 2 cr dr (63) 
where V(r) is the potential experienced by the electron. As the orbital 
angular momentum can be defined in terms of the velocity v by the 
expression 





l= Mer X V, (6.4) 
B l dv 
B-- 2 merc? dr ` (65) 


The factor of c ? means that the resulting magnetic field, and hence the 
energy shift, is small. 

The energy shift due to interactions is given by — u : B. The spin-orbit 
interaction contribution to the Hamiltonian operator can be considered as 
a perturbation to the standard atomic Hamiltonian used in Eq. (4.1). It can 
be written as 


fso = +f (nf £. (6.6) 


For a many-electron atom one can write 


Hso = L.$, (6.7) 





where L and $ are the total orbital and spin angular momentum operators. 
A(L, S) is a constant for a given term, that is a given configuration and 
values of L and S. The spin-orbit coupling is still proportional to r1 4*, 
which becomes r^? for a pure Coulomb potential. This means that A is 
largest for states with small n. 

The energy shifts which result from spin-orbit interaction can be eval- 


uated by considering the expectation value of the spin-orbit Hamiltonian: 


AEso = [ wr Asotir, (6.8) 
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where the integral runs over all space and spin coordinates of the 
electronic wavefunction V. The wavefunction 'V for a complex atom is 
usually not known. However within a non-relativistic framework, the 
quantum numbers E, L, S, J and My, are constants of motion and can 
therefore be taken as known. In particular the angular momentum quan- 
tum numbers are obtained as eigenvalues of the following operator 
equations: 


[^v = L(L4- 1)g?w 
Sw = S(S - 1)? w, (6.9) 
f*-Jg ng. 





Furthermore, by definition the angular momentum operators obey the 
relationship 





]-L-$. (6.10) 
This means that 
P=L+8 L+S) 
= 2484288 
-L-S= (f - 2-8). (6.11) 


One can use this operator relationship to evaluate the effect of L- $ acting 
on the wavefunction: 


I> 


M = s(P -Dp-$ 


e 
= Sg: 1) - L(L4-1) - S(S + 1)”. (6.12) 





The expression for the energy shift due to spin-orbit coupling is therefore 


Al S) 





AEso = ————[J(J 4-1) — L(L - 1) —- S(S + 1)]. (6.13) 


This derivation follows the work of Alfred Landé (1888-1976). 

Note that Hund’s third rule (see Sec. 4.10), means that A(L, 5) is posi- 
tive for an atom with a less than half-full shell, such as sodium, and nega- 
tive for an atom whose shell is more than half full. Use of Eq. (6.13) is best 
illustrated by example. 
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3?P45 — J-30, L-1, S-1/2 AE = +A/2 
3 P 17.2 cm! 






2 J=1/2, L=1, S=1/2 AE =-A 
3 Pin 


Term Level 


Fig. 6.3. Spin-orbit effects in the 3 2P term of sodium. AE gives the shift of the 
level relative to the term. 


Worked Example 1: The Sodium D Lines. 
The transitions 3 PS — 375, in Nat lie in the orange part of the visible 


spectrum. The 3p G 2p) term in sodium is split by spin-orbit interaction, 
as given in Fig. 6.3. 

The sodium D lines are so called because they were so labelled by 
Fraunhofer in his original solar spectrum (see Fig. 1.1). However, the D line 
is actually a doublet and the components are usually labelled: 


D» 5890 A 3p -3s 3 °P; -378;, 
D; 5896 A 3p -3s 3 *P: -3 "5, 


Worked Example 2: The ground state of carbon. 

The ground state of carbon has the term ?P (see Sec. 4.10), and levels given 
by J = 0,1,2. The spin-orbit terms are evaluated in Table 6.3. As A is pos- 
itive for CI, this gives the energy ordering 3p, > ?P, > 3Po. Figure 6.4 
illustrates the observed splittings. 

It can be seen from Fig. 6.4 that the splittings are not exactly in the 2:1 
ratio implied by the values in Table 6.3. This is because the treatment given 
above is highly simplified. There are many other small (magnetic) interac- 
tions which need to be considered in a full treatment. However, for low 
Z atoms and ions, the splitting between the levels approximately follow 
the intervals given by Eq. (6.13) or what is called the Landé interval rule. 


Table 6.3. Spin-orbit interaction terms in the ground 
state ?P term of the carbon atom. 





Level L S J liJg«n-rLL-410-5(541) 


3P, 1 1 2 +1 
3P} 1 1 -1 
3Po 1 1 0 -2 


= 
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Fig. 6.4. Spin-orbit effects in the ground 3P term of carbon. AE gives the shift of 
the level relative to the term as predicted by Landé theory; splittings measured in 
the laboratory are given for comparison. 


The value of A increases rapidly with ion charge, approximately as 
ES Thus, for example, the ground state of C1, considered above, has 
A = 16cm !, whereas the isoelectronic ONI ion has A ~ 95cm ^! in its 
?P ground state. 


6.3 Fine Structure Transitions 


The different transitions between components of spin-orbit states are 
called 'fine structure transitions'. The important selection rule for these 
transitions is AJ = 0, + 1 with the exception that J = 0 — 0 transitions are 
not allowed. 

Following this rule, the transition Na 3?D; 3 — 3 pu 1 has three com- 





ponents because: 
?D por "PS, 2D Pc ^P and D; = P are all allowed but 
"D; = ey is not. 


The three allowed transitions give rise to a triplet structure (see Fig. 6.5). 
The strength of the individual fine structure transitions is given by a 
degeneracy factor times the line strength. Within a non-relativistic treat- 
ment, the line strength is determined completely by the terms involved 
and is therefore the same for all the fine structure transitions. Thus the 
degeneracy factor gives the relative strength of the fine structure individ- 
ual components. The actual formula for these degeneracy factors depends 
on the angular momentum quantum numbers of the upper and lower lev- 
els. A general formulation is fairly complicated and only specific answers 
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Fig. 6.5. The sodium 3 ?2D -3 ?P? triplet, g = 2J +1 gives the statistical weight of 
each level. 


will be quoted here. Tables giving the ratios for all cases can be found in 
Allen's Astrophysical Quantities (see further reading). 

Using these tables the relative strength of the triplet ^D; ; — 3 "PP 1 
transition discussed above is 

*Ds —?P$ : °D; —?P$ : ^D; — P1 of 10:4: 2 
2 2 2 2 z 2 

This corresponds to a ratio of 5:2:1. Any spectrum observed in optically- 
thin conditions should show these ratios. If the spectrum is optically 
thick then a ratio closer to 1:1:1 will be observed. Thus the intensity ratio 
between these transitions give direct information on the optical depth of 
the spectrum. 


6.4 Astronomical Sodium Spectra 


The sodium resonance line, NaI 3s—3p, is prominent in absorption in 
the solar spectrum (see Fig. 6.6), and is known as the sodium D spec- 
trum. For an S-P doublet the intrinsic ratio of line intensities is always 
2:1. When unsaturated, or optically thin, the strength of the absorption 
by the doublet that make up this line should be such that the D; line is 
twice as strong as the D, line (see Fig. 6.7). This ratio is approximately 
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reproduced in Fig. 6.6 which implies that the lines are optically thin in the 
Sun. The sodium D lines are also observed in absorption against starlight 
in the interstellar medium. The D lines are usually very saturated in such 


spectra. 





1. 


Normalised Flux 
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Fig. 6.6. A solar spectrum reflected from the Moon just before a lunar 
eclipse taken at the University of London Observatory. (SJ. Boyle, private 
communication.) 


3p, | d. | g- 4 
¥P ip gm 
3S5 g-2 





D, D, 
5896 A 5890 À 


Fig. 6.7. The sodium D lines, g = 2J + 1, give the statistical weight of each level. 
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Figure 6.8 shows the sodium D lines recorded in absorption. This 
spectrum shows interstellar medium (ISM) sodium ground state atoms 
absorbing against light from reddened stars. The example shown is 
strongly saturated as the intensity ratio for the two D lines is almost equal. 

Figure 6.8 also shows two features which are called diffuse inter- 
stellar bands (DIBs). DIBs are ubiquitous ISM absorption features which 
are present in all reddened lines of sight. They are broader than atomic 
lines and therefore are almost certainly molecular in origin. There have 
been many proposed assignments for these features, most of which have 
been proved to be incorrect and none of which are completely accepted. 
DIBs have been observed since the 1920's and the long-running failure to 
resolve the DIBs problem has been described by Patrick Thaddeus as 'the 
scandal of modern astronomy’. 

Besides the sodium D lines, the Na I 3s—4p doublet at 3302 A can also 
be observed in absorption in the ISM. These transitions are much weaker 
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Fig. 6.8. Interstellar absorption by the sodium D lines recorded towards 
HD 97950 which is the central object in compact galactic cluster NGC 3603, 
recorded using the Anglo-Australian Telescope. Note the two unassigned absorp- 
tion features which are due to diffuse interstellar bands. [Reproduced from 
W.B. Somerville and J.C. Blades, Mon. Not. R. Astron. Soc. 192, 719 (1980).] 
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than the D lines, A(3s—4p) = 2.8 x 10°s~!, compared to A(3s-3p) = 6.2 x 
107 s^!. This means that the 3s -4p doublet can be used when D lines are 
saturated. However they lie in a difficult spectral region where there are 
many telluric (i.e. atmospheric) features due to ozone. 

Higher transitions of sodium are often also observed and provide 
important spectral markers in the atmosphere of cool stars. Figure 6.9 shows 
the spectrum of an L-subdwarf star, a cool star with a mass about 8% of our 
Sun. Absorption features due to Na T, K1 and Rbr are clearly visible in the 
spectrum along with much more complicated molecular features. 

All the alkali metal features appear to be doublets. In fact the NaI 
absorption features arise from a triplet 3p —3d transition: 


8183.26 À 3 oed *D3; 
8194.79 A3 "Pj-3 "Dy; 
8194.82 À 3 "PS -3 "D; f 


2x10716 


LSR1610-0040, sdL: 
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Fig. 6.9. Spectrum of the L-subdwarf star LSR 1610-0400 showing clear spectral 
features due to NaI, KI and Rbr as well as diatomic molecules CaH, TiO, CrH and 
FeH, and water. [Reproduced from S. Lépine, R.M. Rich and M.M. Shara, Astrophys. 
J. 591, L49 (2003).] 
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This appears as a doublet since the splitting in the upper 3 ?D term is too 
small to be resolved. 

The K t and Rbr transitions both belong to the resonance line which is 
a doublet in each case: 


KI 7664.91 À 4°S;—4°P3; 
KI 7698.97 À 4*8, - 4*P3 ; 
RbI 7800.27 À 5°S;-4°*P3; 
Rb 1 7947.60 À 5 E a f 


The significantly larger splitting in the Rb 1 resonance line is due to the 
larger spin-orbit effects in this heavy atom. Indeed the transitions of all 
three species give a nice illustration of the variation of spin-orbit effects 
with principal quantum number n, and atomic number Z. 


6.5 Other Alkali Metal-Like Spectra 


There are a number of astronomically important ions whose ground states 
consist of a single s electron outside a closed shell. These ions have spectra 
similar to those of the alkali metals. Examples include: 


Ca or potassium-like calcium 


The transitions 
H line: 475, - 4?P$ at 3968.47 À, 
2 
K line: 4 E — 4?P$ at 3933.66 À, 
2 


are seen in the solar spectrum. Indeed the labels ‘H’ and ‘K’ are due to 
Fraunhofer. These lines absorb strongly in cool stars (see Fig. 6.10). The H 
line is often blended with the hydrogen Balmer line He at 3970.07 Å. High- 
resolution studies of the structure of K line absorption profiles observed in 
the Sun can give detailed information on the vertical distribution of Ca* 
in the solar chromosphere. 

Calcium has an ionisation potential of 6.1 eV, significantly lower than 
hydrogen. This means that Ca* is found in the ISM, where it can be 
observed unblended with He, since He involves a transition between two 
excited states, where such excited state lines are not present in the cold ISM. 
The unsaturated intensity ratio of the K to H line in absorption is 2 to 1, as 
found for the Na D lines, as this is a general property of S-P doublets. 


94 Astronomical Spectroscopy 





3840 9550 "oen aen 35650 


Fig. 6.10. High-resolution spectrum of the cool, hydrogen-rich white dwarf star 
WD 163344333 showing Ca 1 H and K lines, recorded using the Keck Telescope. 
[Adapted from B. Zuckerman et al., Astrophys. J. 596, 477 (2003).] 


Ca Il also has a strong triplet transition in the red, sometimes referred to as 
the Ca T lines. These lines are: 


8498.0 A 4?p$ -3 *D3 ’ 
2 
8542.1 À 4°P$ -3 °D; , 
8662.1 À 4?P5 -3 °D; . 
Absorptions due to the Ca T lines are clearly visible in the atmosphere of 


cool stars where they are used as an important diagnostic test of metallic- 
ity, the proportion of atomic species heavier than helium (see Fig. 6.11). 


Mg tl or sodium-like magnesium 
The transitions: 
2802.7 À3 "S; -3°P9, 
2 
2795.5 Å 3 °S, -37P5, 


lie in the ultraviolet and so have to be observed by satellite. These lines 
can be routinely monitored by the Hubble Space Telescope (see Fig. 6.12). 
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Fig. 6.11. Call triplet red absorption spectra in stars HD 115444 (upper) and 
HD 102870 (lower) recorded at the Observatoire de Haute Provence. HD 102870 
is an F-type star and has a significantly higher abundance of metals than K-star 
HD 115444. These give rise to the other weaker features. [Reproduced from 
T.P. Idiart, F. Thévenin and J.A. de Freitas Pachieco, Astron. J. 113, 1066 (1997).] 
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Fig.6.12. Ultraviolet interstellar Mg II absorption-line measurements for the sight 
line towards O-star u Columbae, obtained with the Hubble Space Telescope. 
[Adapted from J.C. Howk, B.D. Savage and D. Fabian, Astrophys. J. 525, 253 (1999).] 
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Mg II lines are observed as strong absorptions in stellar atmosphere and 
in the ISM, and as emissions from circumstellar shells. If these lines are 
saturated, the weaker transitions 


1240.4 À 3 "S; -4 "P ! 
1239.9 Å 3 E -4 uu , 
can be monitored instead. 
CIV or lithium-like carbon 
The transitions 
1550.8 A2 ^S; -2 A , 
1548.2 À2 "S; -2 "P3 ; 


were also extensively studied by the International Ultraviolet Explorer 
(IUE) (see Fig. 6.13). These lines are very prominent in emission from cir- 
cumstellar shells and quasars, and can also be seen in absorption. 
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Fig. 6.13. CIV emission and absorption lines towards Q 00135.01-4001 recorded 
using the International Ultraviolet Explorer satellite. This figure illustrates three 
types of absorption lines found in quasi-stellar objects (QSOs). Type A, the 
broad trough about 4230; Type B, the sharp pair of lines at z ~ 1.86; and four 
Type C doublets at Z = 1.62, 1.76, 1.78 and 1.83. [Adapted from R.J. Weymann, 

R.F. Carswell and M.G. Smith, Ann. Rev. Astron. Astrophys. 19, 41 (1981).] 


Alkali Atoms 97 


All the above lines show the characteristic fine structure patterns of 
the alkali metals. This fine structure behaviour has itself been the sub- 
ject of extensive recent studies to try and address the question of whether 
physical constants remain constant with time. The fine structure splittings 
all depend on a dimensionless physical constant called the 'fine structure 
constant', whose value, in atomic units, is simply the inverse of the speed 
of light. By measuring fine structure splittings for a large range of red- 
shifts, it is possible to determine the value of the fine structure constant 
at different epochs. The present evidence point towards a gradual change 
in this constant with time suggesting that the physical constants in our 
Universe did not always have today's values. However this result awaits 
confirmation. 


Problems 


6.1 By considering the levels that arise from the following configurations 
and dipole selection rules, determine what spectral lines would be 
produced in the following transitions: 

(a) Nal 3p- 4d, 

(b) Nar3d 5f, 

(c) Nal 4s—4d, 

(d) Kr4s-4p. 

Comment on the possible observability of these transitions in different 
astronomical locations. 

6.2 A series of transitions in atomic potassium consists of emissions 
from several np levels to the 4s atomic ground state. Transitions 
are observed at 12985cm ! when n = 4; 24701 cm ! when n = 5; 
28999 cm! when n = 6. The series limit, which can be assumed to 
be the same as the ionisation potential of the 4s level, is at 35010 cm '!. 
Calculate the quantum defects of the 4p, 5p and 6p levels. Hence esti- 
mate the wavenumber of the corresponding n = 7 transition. 

6.3 The potassium transitions shown in Fig. 6.9 are a doublet linking 
the ground state to the ^P and AP levels. The two lines lie at 
764.494nm and 769.901nm, respectively. Use this information to 
calculate the constant of proportionality, A(L, S), in the expression for 
the spin-orbit interaction energy, Eq. (6.13). 

6.4 Observations are made of atomic sodium emitting from its 
1s*2s*2p°4d configuration to its 1s?2s?2pÓ3p configuration. At high 
resolution, three transitions are observed. Label these transitions using 
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spectroscopic notation. A lower resolution survey spectrum only 
resolves the Na 4d -3p emission as a single line. What other Na emis- 
sion lines must be present in the spectrum? Suggest an astronomical 
location where the Na 4d —3p transitions might be observed. An astro- 
nomical observation resolves all three transitions in absorption but 
finds them to be of approximately equal intensity. Can one determine 
the column density of Na atoms? 

Configurations of triply ionised carbon, C+, can be written 1s?nl!. In 
terms of Rœ and the quantum defect, give an expression for the energy 
levels of this system. Briefly explain the physical significance of the 
quantum defect and how it depends on n and l. 

The ionisation energy of C?* is 520178 cm! in its 1s?2s! ground state 
and is 217329 cm™! in the 1s?3s! state. Estimate the ionisation energy 
of the 1s?4s! state of C?*. 

The 1s?2s?2p? ground state configuration of O* (oxygen-like nitrogen) 
leads to three terms: *S°, ?P? and 2D° . 


(a) What energy order would you expect the terms to be in? 

(b) Give the levels for each term in full spectroscopic notation. 

(c) Assuming a spin-orbit splitting constant A’ = 30cm ^, sketch the 
pattern of level splitting for each of the doublet terms. 


In spectroscopic notation an atom has a level which is designated ?F9. 
Explain the meaning of this symbol and give values for the angular 
momenta it represents. What other levels arise from the same term? 
Suggest a configuration of atomic carbon that could give rise to this 
level. 

The ?F3 level can emit to levels ?P5, ?P$, ! D and ?Ds. Assuming elec- 
tric dipole selection rules, order these transitions by their probable 
strength, giving your reasons. 


CHAPTER SEVEN 


SPECTRA OF NEBULAE 


* Apparelled in celestial light’ 


— William Wordsworth, Intimations of Immortality (1807) 


Nebulae often have particularly rich spectra (see Figs. 7.1 and 7.2 for exam- 
ples). The spectra contain a wealth of information and can allow pictures 
to be built up of the different regions involved. The spectra of nebulae 
largely involve atomic emissions. These emissions are driven by a number 
of different physical mechanisms which are responsible for creating atoms 
and ions in excited states which allow them to emit. Each mechanism has 
a different spectral signature. The primary mechanisms are: 


Electron collisions: electron collisions largely populate low-lying excited 
states. Excitations are essentially thermal in nature and depend on the elec- 
tron temperature T. 

Recombination: the cascade of emissions following the recombination of 
an electron with an ion leads to emissions from more highly-excited states 
than electron collisions. 

Optical pumping (or resonance-fluorescence): a special mechanism which 
depends on the detailed physics of the system being observed. As discussed 
below, optical pumping is characterised by rather specific line emissions. 


7.1 Nebulium 


In 1918, extensive studies of the emission spectra of nebulae found a series 
of lines which had not been observed in the laboratory. Particularly strong 
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Fig. 7.1. Optical spectra of NGC 6153 from 3540 to 7400 À, obtained for a deep 
(10-minute) exposure using the ESO 1.52 m telescope in Chile. The two spectra 
plotted are (a) obtained by uniformly scanning the long-slit across the entire neb- 
ula, and (b) taken with a fixed slit centred on the central star. [Reproduced from 
X.-W. Liu et al., Mon. Not. R. Astron. Soc. 312, 585 (2000).] 


were features at 4959 A and 5007 Á. For a long time this pair could not 
be identified and, inspired by Lockyer's success with helium, these lines 
were attributed to a new element, ‘nebulium’. 

Ten years after their original observation, Ira Bowen (1898-1973) found 
the true explanation. Bowen realised that in the diffuse conditions found 
in nebulae, atoms and ions could survive a long time without undergoing 
collisions. Indeed, under typical nebula conditions the mean time between 
collisions is in the range 10-10000s. This means that there is sufficient 
time for excited, metastable states to decay via weak, forbidden line emis- 
sions. These lines could not be observed in the laboratory where it was not 
possible to produce collision-free conditions over this long timeframe. 

Bowen identified the doublet as transitions within the 1s?2s?2p? 
ground state configuration of O?*: 


5006.84 A [O r1] 1s72s*2p? !D; P;, 
4958.91 A [O iri] 1s*2s*2p* !D; P, . 
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Fig. 7.2. Continuum-subtracted spectra of planetary nebula NGC 6153 from 4000 
to 4960 À obtained using the ESO 1.52 m telescope in Chile. The spectra show the 
rich recombination-line spectra from C, N, O and Ne ions. The upper spectrum 
was obtained by uniformly scanning the entire nebular surface using a narrow 
long-slit, and the lower one was obtained with a fixed slit centred on the central 
star. The spectra were normalised such that the flux of H9 = 100. [Reproduced 
from X.-W. Liu et al., Mon. Not. R. Astron. Soc. 312, 585 (2000).] 
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Such transitions are strongly forbidden for electric dipoles by the Laporte 
rule. They only occur as weak magnetic dipole transitions. Even for 
magnetic dipole lines these transitions are weak as they are also spin- 
changing, intercombination transitions. They have Einstein A coefficients 
of 1.8 x 10-19 s^! and 6.2 x 107" s, respectively. The line 


4931.23 A [O iri] 1s?2s*2p* !D; —3Po 


is even weaker since this is magnetic dipole-forbidden and occurs as an 
electric quadrupole transition with A = 2.4 x 10 ^s 1. 

Other ^nebulium' lines turned out to be due to [O11] and [N r1] forbid- 
den transitions. 

The forbidden spectrum of the O% ion, [O rri], are often strong in neb- 
ulae. There are a number of [O 111] transitions within the ground state con- 
figuration of O?* 1s?2s?2p? (see Fig. 7.3). As shown in Sec. 4.10, 2p? has 
terms ?P, ! D, 'S in that energy order. All transitions between these terms 
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Fig. 7.3. [Orr] transitions within the o2+ ground state configuration 1s22s22p?. 
All transitions are electric dipole-forbidden by the Laporte rule. 
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are strongly dipole-forbidden by the Laporte rule. The observed transi- 
tions are therefore electric quadrupole or magnetic dipole ones. 

The excited terms associated with the ground state configuration are 
populated by electron collisions. The !S term is higher in energy and there- 
fore requires collisions with more energetic electrons to be populated. This 
line decays via an electric quadrupole transition emitting green light famil- 
iar from the earth's aurora 


4363.2 A [O ir] 1s?2s*2p? 185-1 D;. 


The emission rate from !D relative that to !S depends directly on the 
electron temperature T,. To interpret these data it is necessary to know 
from laboratory studies not only the relevant transition probabilities for 
the observed lines but also the branching ratios, which give the propor- 
tion of ions which decay via these transitions as opposed to the alter- 
native decay routes which also exist. With this information the ratio of 
the unsaturated, observed line strengths can be used to give a mea- 
sure of T,. 

[N 11] is also C-like and therefore has a similar structure to [O m]. [N 11] 
can also be used to measure T,. In the [O 111] and [N 11] examples, the upper 
states have different energies but similar lifetimes. Together they provide 
a useful thermometer. Other transitions can provide a measure of electron 
density. 

A good example of this is the C1 doublet at 1907 and 1909 À. 
These two transitions both involve the same configurations but while the 
stronger transition, 


1908.73 A C m] 2s? !S, -2s2p °P? 


is an intercombination line with Einstein A coefficient of 114 s^1, the other 
weaker line, 


1906.68 A [C 1m] 2s? !S, -2s2p ?P2 


is completely electric dipole-forbidden as AJ — 2. This transition occurs as 
a very weak magnetic quadrupole transition with A = 0.005s~!. The huge 
difference in lifetime for the two transitions is reflected in their critical 
density and hence the relative strength of the two lines can be used as a 
sensitive probe of electron density. 
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7.2 The Bowen Mechanism 


In Planetary Nebulae (PN), many emission lines of O III can be observed at 
visible wavelengths. Some of these transitions are between excited states 
of O?* which are far too high in energy to be populated by collisions. In 
particular there is a whole series of lines in the 3100-3800 À region (see 
Fig. 7.5). These correspond to electric dipole allowed transitions of O III 
but between high-lying states. 

The 3s -3p lines which lie in the 3100-3800 A region follow as a result 
of a cascade from 3p —3d transitions. There therefore has to be some non- 
thermal mechanism whereby the 1s?2s?2p3d configuration is preferen- 
tially populated. Bowen realised that the °P§ level of this configuration 
can be populated through an accidental resonance. The He II Lya line is at 
303.78 À; the transition 


303.80 A O 1 2p? ?P; -2p3d ?P$ 


(note that the closed shell 1s?2s? is assumed) lies at almost the same wave- 
length as He 11 Lya. The small difference is covered by thermal (Doppler) 
shifts within nebula. 

Helium can be doubly ionised in hot nebulae and hence the recombi- 
nation spectrum of He II leads to a plentiful supply of Lya photons. These 
then excite O?* , specifically its 2p3d 3P9 level. This level then decays by 
emission. This mechanism is depicted schematically in Fig. 7.4. 

The Bowen mechanism is an example of a physical process known 
as optical pumping or resonance-fluorescence. Although such processes rely 
on an accidental coincidence in level spacing between different species, a 
number of such occurrences are now known to be astronomically impor- 
tant. For example hydrogen Ly photons at 1025.72 À can pump the 
transition: 


1025.77 A O12p* ?P; -2p?(*s?)3d ?D? . 


In this case, unlike O III, the fine structure splitting in the upper level is so 
small that all three ?D? levels are populated. 

In an ordinary H II region, there are fewer high-energy, ultraviolet 
photons so there are very little He**. Under these circumstances, O?* 
ions are still present and the forbidden [O 111] lines are seen. However the 
Bowen mechanism cannot operate and the allowed lines are not observed. 

O?* transitions can also be observed in the infrared, see Fig. 7.6. In 
this figure note also the Balmer transitions H12a and H13a. The O?* lines 


Spectra of Nebulae 105 






He II "A O III rae 

I 

I 

i 

[| 

2p3p 

| 2p 

I 

i 1 4303.80 
[| 

I 


s 


I 
.62 
303.78 1303: 
I 
I ^ I 
|! 2p3p P^ i 
I Aet S l I 
I sO CUN I I 
I X NON I I 
1 Ns oon I i 
I Me SAY ! [| 
1 A374.44^, ^X l 
I LN E NUS ; I 
I $a Ña I 
N N [ 
s es I l 
an 23 
NV 1?s 2p P 


Fig. 7.4. Schematic partial energy level diagram of O r1 and Hell showing the 
coincidence of Hell Lya and ONI 2p? 3P, -2p3d apo transitions. The Bowen 
resonance-fluorescence lines in the visible and ultraviolet are indicated by the 
solid lines. 


are two forbidden fine structure transitions within the 2p? ground state 
configuration: 


51.8 um [O m] 2p? ?P, -?P»; 
88.4 um [O m] 2p? ?P5 —?P,. 


These are weak magnetic dipole transitions. These infrared transitions 
could not be seen from the ground but were observed by the Infrared 
Space Observatory (ISO). 

It should be noted that the spectrum of NGC 7027 (Fig. 7.6) contains 
transitions assigned to various molecules and molecular ions; the spec- 
tra of these species are discussed in Chapter 10. The planetary nebula 
NGC 7027 has a particularly rich chemistry and the spectrum contains a 
number of lines whose identity is either unknown or uncertain. 
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Fig. 7.5. OII emission lines excited by the Bowen mechanism from planetary 
nebula NGC 3242, recorded using the Isaac Newton Telescope at La Palma. 
[Adapted from X.-W. Liu and J. Danziger, Mon. Not. R. Astron. Soc. 261, 465 (1993).] 
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Fig. 7.6. Infrared spectrum of planetary nebula NGC 7027 recorded using the 
Infrared Space Observatory (ISO). The spectrum is observed on a rapidly vary- 
ing background, as shown in the insert, and the detailed line spectrum is obtained 
after subtracting a polynomial fit to the continuum part of the spectrum. [Repro- 
duced from X.-W. Liu et al., Astron. Astrophys. 315, L257 (1996).] 
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Fig. 7.7. Ultraviolet spectra of planetary nebulae in the Large Magellanic Cloud 
recorded using the Hubble Space Telescope. The Lya emissions are geocoronal in 
origin. [Adapted from E. Vassiliadis et al., Astrophys. J. Suppl. Ser. 114, 237 (1998).] 


O?* transitions are also observed in the ultraviolet (see Fig. 7.7). The 
O11] doublet near 1663 A area pair of intercombination lines, as shown 
by the single bracket, meaning that the electric dipole selection rules are 
all satisfied except for spin conservation, i.e. AS Æ 0. These transitions link 
the lowest quintet state to the ground state, which is a triplet: 


1660.81 A O ir] 1s*2s?2p* ?P, -1s?2s2p? 585. 
1666.15 A O ur] 1s*2s?2p* °P, -1s?2s2p? 585. 


The upper level is just a little higher in energy than the ground configura- 
tion !S. It can therefore be populated by collisions with energetic electrons. 
Other terms with the configuration 2s2p? are 395°, ?Do, ?po, 1P° ang !D°. 
These all have higher energy and so the lines linking them to the ground 
state are at still shorter wavelengths. These transitions can be used to mon- 
itor still higher electron temperatures. 
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7.3 Two Valence Electrons 


The previous chapter considered the spectra of alkali metals and alkali 
metal-like ions which all have a single optically-active electron. The pres- 
ence of two outer shell electrons introduces important new physical 
processes. 

The calcium atom will be used to illustrate these processes. The 
ground state of calcium has two 4s electrons outside a closed shell: 


1s?2s?2p63s?3p64s? . 


Astronomical spectra of CaI can be understood entirely in terms of 
these two outer electrons, which are called valence or optically-active elec- 
trons. One can therefore safely assume that the other electrons are inactive 
and remain frozen in their ground state orbitals. 

Excitation of one 4s electron gives the configuration 4sn! which can 
have spin S — 0 or 1, except for the 4s? configuration itself which gives 
rise to 5 — 0 only. The orbital angular momentum possessed by these 
states is L = |. The resulting energy level diagram for the 4sril states of 
Ca is given by Fig. 7.8. For these states, the ionisation limit, represented 
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Fig. 7.8. Main branch energy level diagram of calcium. The diagram gives 4snl 
configurations of the outer two electrons. These converge to Ca? in its 4s ground 
state plus a continuum electron. 


Spectra of Nebulae 109 


by n — oco, is Ca* in its 4s configuration, which is the ground state of the 
ion. The diagram represents what is called the ‘main branch’ of the Ca 1 
spectrum. However it is not sufficient to give all the observed transitions 
of Cal. 

To explain how the other transitions arise it is necessary to con- 
sider bound states of the calcium atom which have both outer elec- 
trons excited. For example, consider the case where one electron is in 
the 3d orbital and the other electron jumps, giving a configuration 3dnl. 
This gives rise to many new electronic states which involve excitation 
of two electrons. These are depicted in Fig. 7.9. The 3dnl states go to 
an ionisation limit as n — oo of Ca* in its first excited 3d configura- 
tion. A complete understanding of the spectrum of Cal also requires 
4pnl configurations going to the Cat 4p configuration at their ionisation 
limit. 

Structure of this sort is important in spectra of nearly all complex 
atoms. An exception is helium for which all states with two excited elec- 
trons are unbound. 


I LUI 


/ / / 


Ca* 3d 






3d6p 


3d5p 


Energy 








3d4p 


4s? 


Fig. 7.9. Energy level diagram of calcium. The diagram gives 4snl and 3dnl con- 
figurations of the outer two electrons. These converge to Cat in its 4s ground state 
and Ca* in its 3d first excited state respectively, plus a continuum electron. 
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7.4 Autoionisation and Recombination 


As illustrated by the various series for calcium, each series of energy levels 
converges towards its relevant continuum level. That is, at the ionisation 
limit for each series, there remains an ion in some state and an electron 
in the continuum. However, series involving two-electron excited states, 
ie. ones which are not on the main branch, result in ion states which are 
also excited states. As illustrated by Fig. 7.9, the higher states of such series 
lie in the continuum of the ground state of the ion. In other words, there 
are states with two excited electrons which lie above the first ionisation 
threshold. Such states have the same energy as a state of the ion plus a free 
electron. Free electrons lie in the continuum since their energy levels are 
not quantised. 

When there is resonance between two states with the same energy, the 
system goes rapidly from one state to the other and back. ‘Bound’ states in 
the continuum are therefore called 'resonances' since they (always) lie at 
the same energy as an unbound, continuum state. The states of radioactive 
nuclei are also resonances; these states of the nucleus are ones which spon- 
taneously decay with emissions of radioactive particles or rays. Similarly, 
electronic resonances can change spontaneously from being a doubly- 
excited state to an ionised state and an free electron: 


A eA Es (7.1) 


This process is called autoionisation and is illustrated schematically 
by Fig. 7.10. It is a radiationless transition in which one electron jumps 
to a lower level and the other electron escapes to the continuum 
simultaneously. 

Transition probabilities between bound states and resonances obey the 
same selection rules as transitions between truly bound states. However 
resonances do not have a precise energy but have some width or spread of 
energies, AE. This is a consequence of Heisenberg’s Uncertainty Principle: 
as the resonance has a finite lifetime, 7, it cannot have a definite energy. 
The Uncertainty Principle gives 


AE.T > k . (7.2) 
Allowed transitions between bound states are generally much 


stronger than bound-free transitions. This leads to the resonant enhance- 
ment of autoionisation (see Fig. 7.11). As allowed transitions to resonances 
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Fig. 7.10. A schematic representation of autoionisation: a radiationless process 
whereby a system with two excited electrons loses one to the continuum while the 
other simultaneously jumps to a lower energy level. 
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Fig. 7.11. Resonance-enhanced ionisation involves an allowed transition from a 
bound to a resonance state followed by a radiationless transition to the continuum. 


are also stronger, the presence of a suitable resonance can increase signifi- 
cantly the overall ionisation rate in a given radiation field. See Fig. 7.12 for 
example. 

Resonances therefore play an important role in photoionisation. Astro- 
nomically they play an even more important role in recombination. As 
discussed in Sec. 3.9, recombination is an important source of emission in 
regions where some or all of the atoms are ionised. For complex atoms 
there are two possible recombination mechanisms: 


(1) ‘Direct’ radiative recombination, as discussed for H-like atoms in Sec. 3.9. 
(2) Dielectronic recombination: this can be thought of as the reverse of res- 
onant photoionisation. A continuum electron with the correct energy 
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Fig. 7.12. Photoionisation cross section, in units of 10-2 m2, for the 2s2p? ?D 
excited term of C* as a function of photon energy calculated by A.R. Davey, 
PJ. Storey and R. Kisielius, Astron. Astrophys. Suppl. Ser. 142, 85 (2000). The rate 
of this process is dominated by the near-threshold resonance. (P.J. Storey, private 
communication.) 


forms a resonance by undergoing a radiationless transition. This reso- 
nance can then decay by autoionisation or it can become a truly bound 
state by the emission of a photon. 


The process radiative recombination can thus be written as 





At+e — A" — A* hv 
sA deer. (7.3) 


As autoionisation is a rapid process, only dipole-allowed radiative tran- 
sitions can compete with it. As the resonance is a two-electron excited 
state, dielectronic recombination often results in different transitions, i.e. 
photons with different wavelengths, than direct radiative recombination 
which only goes via the main branch (see Fig. 7.13). 

Dielectronic recombination relies on the presence of resonances which 
are low-lying enough to be accessible by thermal electrons. Whether such 
resonances exist is a property of the physics of the particular ion in 
question. Dielectronic recombination is particularly important in some 


Spectra of Nebulae 113 


[dll 


dielectronic 






Continuum 


recombination 
line 


| recombination line 


main branch 


Fig. 7.13. Dielectronic and direct (or radiative) recombination. 


cases, such as Ct. C* has a very large and low-lying resonance (see 
Fig. 7.12), formed by the temporary capture of an electron by C?*. CH 
recombination spectra are completely dominated by lines from dielec- 
tronic recombination. 

Snapshot spectra, which only survey the strongest lines, show neb- 
ula emissions dominated by hydrogen recombination lines plus a few 
atomic (ion) transitions. However deep spectra which probe weak tran- 
sitions can show a huge number of lines. See Fig. 7.2 for example. These 
lines generally arise from atomic recombination spectra. But why are 
there so many of them? This question is best answered by considering an 
example. 

The strong C1 line at 4267 À in Fig. 7.2 is due to the recombination 
spectrum of C11. The C?* core is a closed shell with the configuration 
1s?2s?(1S). The resulting recombination lines are due to transitions of the 
form 


Cir 18228? (1S) nl -18?28?(1S) n'l . 


Thus the line at 4267 A is a 3d - 4f transition. This transition is seen as only 
a single line as the (triplet) fine structure is not resolved. 

Figure 7.2 shows a large number of weaker Or II recombination lines. 
These lines also arise from the 3d — 4f recombination transition but occur 
at many wavelengths. The ground state of O?*, which gives rise to the 
main branch recombination lines, is an open shell: 1s*2s*2p? (P). When 
recombined with an extra electron, the open core can couple with the extra 
electron to form a number of different terms as both the core and the extra 
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electron have spin and orbital angular momenta. Thus 


SP plus 4f ?F gives 6 separate terms and 18 distinct levels; 


?P plus 3d ?D also gives 6 terms and 16 distinct levels. 


There are therefore a large number of possible transitions between the dif- 
ferent terms and/or levels. Which transition is the strongest? Since for 
a given atom or ion all 3d — 4f transitions have similar line strengths, 
the relative strength of the transitions are given by considering statistical 
weights. For the OII example above, the strongest emission lines which 
come from the 4f configuration arise from the G° term which has a sta- 
tistical weight of g = (2L + 1) x (25 +1) = 9 x 4 = 36. In contrast the ?D? 
term has a statistical weight of only g — 5 x 2 — 10. 


Problems 


7.1 By considering the levels which arise from the following configura- 
tions and dipole selection rules, determine what spectral lines would 
be produced in the following transitions: 


(a) Ca14s? - 4s4p, 
(b) Ca14s4f — 4s5d? 


Comment on the possible observability of these transitions in different 
astronomical locations. 
7.2 Give the ground state configuration of the N* ion and explain which 
of its terms you would expect to have the lowest energy. 
A series of N I recombination transitions can occur as a result of radia- 
tive recombination of ground state N* ions and an electron. Many 
optical 5f—4d transitions are observed. Explain why so many transi- 
tions are seen. Neglecting fine structure effects, how many allowed 
transitions would you expect? Which one is likely to be the strongest? 
Use Figs. 7.3 and 7.4 to construct a partial energy level diagram for the 
O?* ion, with energies in eV. 
What is the statistical weight of the atomic level *5+1L)? Give the sta- 
tistical weight for the O?* ion states in your diagram. 
Assuming the population of the 2p? ?P, ground state to be 1, use the 
Boltzmann distribution for thermal populations, and a temperature of 
T — 10000K, to give the occupation of the other O?* levels in your 
diagram relative to the ground state. 


7. 


ies) 


Spectra of Nebulae 115 


Discuss the significance of this result in relation to the observed O III 
spectrum in a nebula with T, = 10000 K. 
The following will be useful: 


1cm-! = 1.23985 x 1074 eV, 
Boltzmann's constant, k = 8.62 x 10^ eV. K^. 


CHAPTER EIGHT 


X-RAY SPECTRA 


X-rays are high-energy photons which have shorter wavelengths than 
those in the ultraviolet. X-rays can be produced as a result of transitions 
involving an inner shell in heavy atoms and their name was initially asso- 
ciated with this particular process. Indeed X-rays produced in this fashion 
in X-ray tubes as a result of some heavy metal being bombarded by high- 
energy electrons remain important for routine medical inspections of 
bones inside living people. 

Photons with even higher energies than X-rays are called y-rays. Like 
X-rays, 7-rays were originally associated only with the process that pro- 
duced them. In the case of y-rays this was the radioactive emissions from 
nuclei. However the terms X-ray and y-ray are now just used to label parts 
of the electromagnetic spectrum. Indeed astronomically neither of the two 
original processes are major sources of high-energy photons. 

The previous chapters have considered spectra associated with the 
movements of outer electrons which result in a so-called 'optical' spec- 
trum. For an atom, with (fairly) high Z, transitions which involve move- 
ments of inner shell electrons occur at much shorter wavelengths. In 
principle, so long as only one electron jumps, there is no difference 
between the two cases and much of the theory already outlined still 
applies. In particular, X-ray transitions obey the dipole selection rules 
given in Table 5.1. 

If an inner shell electron is removed from an atom by any pro- 
cess including high-energy electron impact, a vacancy is created. This 
inner shell vacancy can then filled by an outer electron jumping down 
(see Fig. 8.1). Such transitions give rise to characteristic X-ray emission 
lines. 
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Fig. 8.1. X-ray line spectrum caused by the removal of an inner shell electron by 
electron impact followed by the subsequent relaxation of an outer shell electron 
with emission of an X-ray photon. This is how X-rays are usually produced in the 
laboratory. 


Transitions involving K shell (n = 1 electrons, see Table 4.1) involve 
most energy: ‘hard’ X-rays or about 10° eV for heavy atoms. These energies 
correspond to wavelengths, A, of about 0.1 A. Higher shells give progres- 
sively 'softer' X-rays, with wavelengths up to tens of À. Wavelengths 
between about 100 A and 1000 A are often described as the ‘extreme ultra- 
violet' or EUV. 

In fact, most X-rays produced in astronomical sources are not asso- 
ciated with particular transitions but are part of the continuum spec- 
trum produced in high-energy environments such as active galactic nuclei 
(AGN). For example, it is thought that AGN generate X-rays from the 
atoms being accreted by a central massive black hole. These spectra can 
be extremely bright, producing more radiation than the rest of the entire 
galaxy; they give direct information on the distribution of accreted mate- 
rial. The solar corona provides a nearby source of X-ray spectra. Spectral 
lines come largely, not from inner shells of neutral atoms, but from highly 
ionised atoms. See Fig. 8.2 for example. 

Consider the spectrum of FeXXVI as an example. Fe??* is a one- 
electron system with Z — 26. This means that its spectrum is similar to 
hydrogen but with all energies and transition frequencies scaled by Z?. 
For H, the Lya photon has an energy of 10.2eV, while for Fe XXVI, Lya 
has an energy of 10.2 x (26)? eV — 7 keV. Emission of FexxvI Lya can 
be seen in the X-ray emission coming from intergalactic plasma which 
lies within clusters of galaxies. The continuum emission from this plasma 
is thermal emission with a temperature in the region of 108K. At this 
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Fig. 8.2. X-ray spectrum of the entire solar disk recorded using a rocket-based 
experiment. [Adapted from M. Malinovsky and L. Heroux, Astrophys. J. 181, 1009 
(1973).] 


temperature, iron is predominantly ionised to Fe??* (see Fig. 3.10). This 
spectral region also contains lines from other highly-ionised species such 
as Fe XXV (see Fig. 8.3). This figure also shows the spectral signature of 
Fer Ka transitions. These transitions arise from ionisation of a 1s electron 
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Fig. 8.3. X-ray spectrum of Seyfert galaxy NGC 7314 showing various iron tran- 
sitions recorded using satellite Chandra [see T. Yaqoob et al., Astrophys. J. 596, 85 
(2003)]. (T. Yaqoob, private communication.) 


in Fe and then relaxation of a high-lying electron which emits at one 
of these characteristic wavelengths. As discussed above, such transitions 
have long been observed in the laboratory but had been thought to be 
unimportant astronomically until the recent batch of satellite-borne X-ray 
experiments. 

Astronomical X-rays cannot be observed from ground. Early obser- 
vations were made from rockets but now several satellites have been 
deployed. Recent space missions include SOHO and XMM-Newton 
launched by ESA and Chandra launched by NASA. Other space satellites 
are planned. These missions have substantially changed our view of the 
high-energy regions of the Universe as they have opened up many objects 
for detailed observations of X-ray spectra which previously could only be 
performed on the Sun, which has a very large flux of X-rays. 


8.1 The Solar Corona 


The visible spectrum of the solar corona shows a number of strong emis- 
sion lines which for a long time could not be identified. These emissions 
were labelled ‘coronium’, although they were not generally thought to 
belong to a new element. 

Seventy years after their original observation in 1869, the Swedish 
astronomer Bengt Edlen (1906-1993) showed that the lines were due to 
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forbidden transitions in a number of highly-ionised atoms. The three most 
prominent lines are: 


6374 A [Fe x] 3s73p° ?P$ —?P?, 
2 2 
5302 A [Fe xiv] 3s?3p *P§ —?P3, 
2 


3 
2 
5694 A [Ca xv] 2s?2p? °Po -?P . 


These transitions all occur within a single configuration and, indeed, 
within a single term. They are thus all strongly forbidden by electric dipole 
selection rules as they do not satisfy the Laporte rule. The transitions are 
all weak magnetic dipole fine structure lines. The green [Fe XIV] line at 
5302 À is the strongest of these coronal lines. 

Transitions involving highly-ionised atoms such as Fe?*, Fe?* and 
Cal^* will occur at high temperatures. Inspection of Fig. 3.10 suggests 
temperatures in the region of one million Kelvin. This temperature is 
much higher than the assumed temperature of the solar corona so Edlen's 
assignments were at first very controversial. However his analysis is now 
accepted but the mechanism for heating the corona remains not fully 
understood. 

Many other coronal lines from highly-ionised metals can be observed 
in the solar spectrum in the extreme ultraviolet (EUV) and X-ray. Figure 8.4 
shows the X-ray spectrum of the nearby star Capella, which is a G-star like 
the sun. It shows many emission features due to highly-ionised ions. 


8.2 Isotope Effects 


As discussed in Sec. 33, the energy levels of an atom depend on the 
reduced mass of the atom. For a hydrogen-like atom with nuclear mass M, 
the reduced mass is: 
|. Mm, 
eae Me 





(8.1) 


Isotopes have the same atomic number Z and the same electronic struc- 
ture, but have different nuclear masses, as expressed by the atomic 
weight, A. For example, deuterium is an isotope of hydrogen. In standard 
notation, the integer atomic mass (i.e. the total number of protons and 
neutrons) is given as a leading superscript on the atomic symbol. Thus 
hydrogen is 1H and deuterium, often designated D, should more correctly 
be given as 7H. 
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Fig. 8.4. X-ray spectrum of binary star Capella recorded using the XMM-Newton 
satellite [see M. Audard et al., Astron. Astrophys. 365, L329 (2001)]. (M. Audard, 


private communication.) 
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For H-like atoms, the energy levels and frequencies simply scale 
according to the reduced mass factor [see Eq. (3.8)]. For complex atoms 
the situation is not so simple. There is no longer a single scale factor as the 
mass effects of several electrons have to be allowed for simultaneously. 
The effect is also smaller than in hydrogen and immeasurably small for 
atoms with Z beyond about 30. 

The measurement of the isotopic composition of a particular environ- 
ment is important as it varies with the age of the Universe and the nuclear 
reactions occurring in a particular star. Isotopes thus provide an important 
measure of local nuclear reactions. For example, deuterium was formed in 
the early Universe but is burnt early in the life of all but lightest stars. The 
abundance of deuterium thus decreases with the age of the Universe and 
its presence can be used as an indicator that a brown dwarf lies below the 
critical mass necessary to burn D. Similar analyses can be conducted with 
other isotopes, although these are generally made rather than destroyed 
by stellar nuclear synthesis. 

For most elements the best way of determining isotopic abundance is 
via molecular spectra; methods to do this are discussed in Sec. 10.1.1. For 
very heavy atoms, however, the effect of nuclear size becomes important. 

So far we have always treated the nucleus as a point charge. Although 
very small, a factor of about 10? smaller than a typical atom, the nucleus 
does have a finite size. Furthermore the nucleus slightly distorts the elec- 
tronic wavefunction and associated energy levels. This distortion is par- 
ticularly important for electrons in s (l = 0) orbitals as these orbitals have 
nonzero densities at the nucleus (see Fig. 3.1) and for heavy elements 
which have larger nuclei. 

An example is the element mercury, Hg, which has Z — 80, and so is 
heavier than iron and is therefore only present in significant quantities in 
chemically peculiar stars. Mercury can exist in one of seven stable isotopic 
forms with atomic weight A equal to 196, 198, 199, 200, 201, 202 and 204. 
Hg* has a prominent line: 


3984 A Hg 11 5d’6s? - 5d" 6p "D; -?P$ . 
2 


Note that this transition is formally forbidden as it is a two-electron tran- 
sition, but in practice, strong configuration interaction (CI) effects ensure 
that it is actually strong. The change in the occupancy of the 6s orbital in 
this transition makes it sensitive to nuclear size effects, since all s orbitals 
have some probability of the electron overlapping the nucleus. When 
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observed at very high resolution, the 3984 À line can be separated into 
different wavelength components ranging from 3983.77 A for P6Hg* to 
3984.07 for Hg+, allowing the isotopic composition to be determined. 
On earth the commonest isotopes of mercury are 7°°Hg and ?'?Hg; yet 
for no star are the terrestrial abundances of mercury isotopes observed. 
Instead, chemically peculiar stars show a range of isotopic compositions 
but with a tendency for the heaviest isotope ?'"Hg to be the most 
abundant. 

The analysis of mercury spectra discussed above is complicated by the 
fact that the spectrum of each isotope has its own hyperfine structure due 
to nuclear spin effects. Isotopes differ not only in their atomic weights but 
also in their nuclear spin, I. For example hydrogen, 1H, isa single proton 
and therefore has I = $, but deuterium, ?H, has I = 1. Similarly the nucleus 
of *He, which is an alpha particle, has zero spin, I — 0, but ?He has I = 1. 
For heavy elements, I can be larger and values up to 7 or 8 are found. 

As discussed for hydrogen (see Sec. 3.14), coupling the nuclear spin, 
I, with the total (electronic) angular momentum, J, leads to the ‘final’ 
angular momentum quantum number, F: 


E-14J. (8.2) 


Weak splittings of transitions according the value of F are called the hyper- 
fine structure. The hyperfine structure gives isotopic information, but in 
practice, is only fully resolved in ultrahigh-resolution spectra and is there- 
fore not generally astronomically useful. Use of such structure for different 
isotopes also requires very detailed laboratory data which is not always 
readily available. 


Problems 


8.1 Give an expression for the energy levels of the hydrogen-like atom 
of nuclear charge Z in terms of the Rydberg constant Ræ. Obtain 
the wavenumber of the Lya transition of hydrogen-like oxygen, O7*. 
From what astronomical environments would such transitions occur 
and how might they be observed? 

8.2 Mercury has atomic number Z of 80. Give the full configuration of the 
5d"6s? ?D; state of Hg*. The ground state of Hg* has the configura- 
tion 5d'°6s. Obtain the term and level for this state. Why is the 5d?6s? 
?D; state metastable? 


CHAPTER NINE 


MOLECULAR STRUCTURE 


‘Ban dihydrogen monoxide’ 


— Spoof webpage www. circus.com!~nodhmo/ (2004) 


It is easy to record the astrophysical spectra of atoms and ions since 
they are present in the atmospheres of hot stars which shine brightly. 
Molecules, conversely, are found in cooler and less active regions which 
can make them more difficult to observe. Nonetheless, molecules can be 
found in many different astronomical environments. These range from 
planetary atmospheres and comets to cool stars and sunspots; from cold, 
giant interstellar clouds to planetary nebulae and even the tori which form 
round active galactic nuclei. 

The spectra of molecules are studied at wavelengths from the ultra- 
violet to the radio. However it was the development of radio astronomy 
that led to the realisation that substantial portions of our galaxy, and hence 
the Universe, are dominated by molecular processes, indeed complicated 
ones. So far over 120 molecules have been detected in the interstellar 
medium by direct observation of their spectra. These molecules contain 
up to 13 atoms and their spectroscopy is rich and complex. There have 
also been claimed detections of significantly larger molecules, but these 
are harder to verify. 

The introduction to molecular structure and spectroscopy presented 
here will largely be restricted to diatomic (two-atom) molecules as they 
show most of the key features of molecular spectra. For a discussion 
of the spectroscopy of larger (polyatomic) molecules, see Banwell and 
McGrath (1994), or Bernath (1995) in further reading. Both these textbooks 
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give a more advanced and a more general treatment of the laboratory 
spectroscopy of molecules. 


9.1 The Born-Oppenheimer Approximation 


Molecules formed from two atoms are called diatomics while those 
formed from more than two atoms are described as polyatomic. Because 
molecular binding energies are relatively small, i.e. generally less than 
ionisation energies, molecules are only found in cooler, or less active, astro- 
nomical environments. Polyatomic molecules are only a significant com- 
ponent of matter at temperatures below about 4000 K. Diatomic systems 
can survive to somewhat higher temperatures and may be found in envi- 
ronments with temperatures up to about 8000 K. Of course, these figures 
are a rule of thumb with details depending on the molecule in question 
and the physical conditions of the environment. 

The structure, and hence the spectra, of molecules are more compli- 
cated than atoms in two ways: 


(1) There is no single charge centre about which electrons move. The elec- 
tronic wavefunctions therefore have lower symmetry, making them 
harder to calculate and harder to work with. 

(2) The nuclei themselves move, giving rise to both rotational and vibra- 
tional motions of the atoms within the molecule. These motions give 
rise to discrete spectra. 


To a very good approximation, known as the Born—Oppenheimer 
approximation, the motions of the electrons and the nuclei can be com- 
pletely separated. The validity of this approximation relies on the fact that 
the electrons are much lighter than the nuclei. For hydrogen, which is 
the worst case, My ~ 1836 me. This means that the electrons move very 
much faster than the nuclei and can be considered to relax instantly to any 
change in the positions of the nuclei. To use an everyday scale analogy, it 
is like flies buzzing round an elephant — as the elephant moves the flies 
move with it. 

Within the Born-Oppenheimer approximation, one separates the 
wavefunction for the motions of electrons from the wavefunction for the 
motions of the nuclei [see Eq. (9.7)]. One can then consider the electronic 
wavefunction separately for each position of the nuclei, as if the nuclei are 
held fixed. The electronic energy associated with each of these electronic 
wavefunctions gives the familiar potential energy curves which show the 
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interactions between two atoms. Such curves, examples of which are given 
in Fig. 9.2, only exist within the Born-Oppenheimer approximation. 

To a less good approximation, one can also consider separately the 
two types of nuclear motion: vibration and rotation. The energies of these 
motions are such that the energy associated with motion of the electrons 
(the electronic energy) is always very much greater than the energy con- 
tained in vibrational motion which, in turn, is greater than the energy of 
rotation. This ordering is very useful when considering molecular struc- 
ture but it is important to remember that in molecular spectra a particular 
transition can involve changes in more than one type of motion. For exam- 
ple, one gets rotational structure on vibrational transitions, and electronic 
transitions have vibrational structure and also fine structure due to simul- 
taneous changes in rotational motion. The rotational fine structure in elec- 
tronic transitions produces many lines close together. This fine structure 
can only be resolved at very high resolution which is often not available 
for astronomical spectra. If the rotational fine structure is not resolved, 
then the result is characteristic molecular band spectra. 


9.2 Electronic Structure of Diatomics 


Within the Born-Oppenheimer approximation, consider the electronic 
state of a diatomic molecule as a function of the separation, R, between 
the two atoms. Figure 9.1 shows typical coordinates for a diatomic, AB, 
where the atoms have nuclear charges Z4 and Zp respectively. For clar- 
ity only two electrons are shown explicitly, so this figure can be taken to 
represent molecular hydrogen, H;. 

When considering the electronic state of the molecule AB as a function 
of R, the following considerations come into play: 


As R — 0, at very small R the overall interactions are strongly 


repulsive. There is repulsion due to nuclear-nuclear interaction 


e 
1 I15 ez 


Fig. 9.1. Coordinates for a diatomic molecule, AB, with two electrons. 
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whose potential depends on 4448. There is also repulsion due 


to the electron-electron interactions, which also behave approxi- 
mately as R^!. However the most important interaction at small R 
in a many-electron molecule is due to exchange forces. Exchange 
interactions arise as a direct consequence of the Pauli Principle 
(see Sec. 4.3), which strongly forbids any attempt to force the two 
electrons to occupy the same space. This behaviour is similar to 
that found in collapsed stars where degeneracy pressure is a con- 
sequence of the Pauli Principle. It is the exchange forces that cause 
the atoms in molecules to behave almost as hard spheres when the 
molecule is compressed. The repulsion due to exchange forces at 
small R increases approximately as an exponential. It is therefore 
best modelled as e°’, where a is a system-dependent constant. 


As R — co, the molecule is pulled apart and it separates into 
atoms in a process known as dissociation. The energy of the system 
at dissociation is clearly just the sum of the atomic energies. 


At intermediate R, to get binding there must be some region 
of R where the molecular energy is less than the sum of the 
atomic energies. In this case the electronic state is described as 
'attractive' and there is a minimum in the potential energy curve. 
The attraction can arise from interactions such as the increase 
in the electron-nuclear attraction brought about by bringing the 
atoms together. States can also be purely repulsive. Such states are 
referred to as 'dissociative' since they do not support bound states 
and will always lead to dissociation. 


Figure 9.2 gives potential energy curves for the lower-lying electronic 
states of molecular hydrogen. The scheme used to derive labels for these 
curves is discussed in the next section. Note that the electronic potential 
is by convention denoted V. Two important quantities for the curves in 
Fig. 9.2 are the equilibrium bondlength and the dissociation energy. Re, 
the equilibrium internuclear separation is defined by 


dV 
TR p =" (9.1) 


that is Re is the value of R at the minimum of the potential. D, is the dis- 
sociation energy of the molecule. It is the minimum energy required to 
fragment a molecule. For a diatomic molecule it can be defined by 


D, = V(R,) — V(R = oo). (9.2) 
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Fig.9.2. Potential energy curves for the lowest seven electronic states of molecu- 
lar hydrogen, H5. The lowest two curves both dissociate to two ground states, 1s 
hydrogen atoms; the higher curves dissociate to one H atom in its 1s ground state 
and the other into either a 2s or a 2p excited state. [Reproduced from D.T. Stibbe, 
PhD thesis, University of London (1997).] 


9.2.1 Labelling of electronic states 


Electronic states of atoms are labelled according to the so-called spectro- 
scopic notation (see Sec. 4.8). The notation for electronic states of molecules 
follows a similar system but is more complicated because of their lower 
symmetry. In general terms the labelling scheme used is based on the 
symmetry of the molecule in question. Here only diatomic molecules will 
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be considered. Of course, all diatomic molecules are linear molecules and 
only two separate symmetry cases need to be considered: 


Homonuclear diatomics, where both atoms are the same, e.g. H5, N2, O2. 


Heteronuclear diatomics, where the atoms are different, e.g. CO, HF, CH. 


As the electron spin is usually not strongly coupled to the frame of the 
molecule, the treatment of spin follows very much the same principles as 
in the atomic case. Each electron has its individual spin angular momen- 
tum, s;. These can be summed to give a total spin angular momentum, 
S, again remembering that paired electrons in closed shells make zero 
contribution to this sum. The electronic states of molecules are designated 
by their spin multiplicity, 2S + 1, which is given as a leading superscript, 
exactly as in the atomic case. 

In an atom the treatment of the individual orbital angular momenta, 
li, follows along similar lines to the treatment of spin. However molecules 
are not spherical and the orbital angular momentum of the individual 
electrons is no longer a conserved quantity. For diatomic molecules, the 
total orbital angular momentum L is strongly coupled to the nuclear axis. 
It is therefore necessary to consider the components of L, designated A, 
along the diatomic nuclear axis which, by convention, is taken to define 
the z-axis of the system. What this means is that while the value of the 
total orbital angular momentum in a diatomic molecule can change, its 
projection onto the diatomic axis is conserved. As the projection of L onto 
z-axis can be either positive or negative, states with A # 0 are 
twofold degenerate while 2 states, which have A = 0, are singly 
degenerate. 

Electronic states are labelled by their value of A rather than L. Values 
of A are denoted using the Greek letter equivalent of the Latin letter 
used to denote L (see Table 9.1), thus &, IT, A are equivalent to the atomic 


Table 9.1. Letter designations for projected 
total orbital angular momentum quantum 


number, A. 

A= 0 1 2 3 4 .. 
Orbitals o m ô ó qa 
States 2 II A ® D.. 
Degeneracy 1 2 2 2 2: 
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symbols S, P, D. Again, uppercase letters are used to denote many electron 
states and single electron orbitals are labelled using lowercase versions of 
the same symbols. 

This means that: 


15 denotes a state with S = 0 and A — 0, 
?II denotes a state with S = 1 and A = 1, etc. 


For most (stable) diatomics, the electronic ground state is a closed shell, 
meaning that it is 15. Examples include H2, Nz and most other homonu- 
clear diatomics. The exception is O2 which has a 3x ground state. CO and 
many other heteronuclear diatomics with an even number of electrons 
also have !X ground states. Diatomics with an odd number of electrons 
usually have S = 3. For example H7, CH* and CN all have 7 ground 
states. Ground states need not have A — 0 (2), just as some atoms have 
ground states with L > 0. Examples include CH, OH and NO which all 
have ?II ground states. Molecules in electronic states with A > 0, and in 
particular those with A — 1, have extra lines in their spectra due to a pro- 
cess called A-doubling. This is a significant extra complication which will 
be discussed in Sec. 10.4 because of the astronomical importance of the 
CH and OH molecules. Indeed, CH was one of the first molecules to be 
detected outside our solar system. 


9.22 Symmetry 


There are two symmetry properties of certain diatomics which are similar 
to the parity of an atom. These contribute to electronic selection rules and 
hence are included in the electronic state designation. 

Homonuclear diatomics have identical nuclei. This means that all 
properties must be unchanged with respect to swapping the two nuclei. 
The usual way of treating this is to consider the behaviour of the electronic 
wavefunction when the molecule is inverted through its centre-of-gravity. 

Consider an electronic wavefunction, v; (R4, Rg), where R; is the vec- 
tor connecting the molecular centre-of-gravity with nucleus A, and Rg 
connects the centre-of-gravity with nucleus B. As the two nuclei are iden- 
tical the probability distribution of the electronic wavefunction must be 
unchanged with respect to inverting the molecule A - B — B — A, which 
has the effect of interchanging the nuclei: 


|We(Ra, RB)? = [jv RA, —Ra)P. (9.3) 
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This equation has two solutions 





We(Ra, Rs) = ~Ye(—Ra, —Rg) . (9.4) 


Both solutions are acceptable. Electronic states which are positive with 
respect to interchange are known as 'gerade', the German word for 'sym- 
metric’, or more simply ‘g’. Electronic states whose wavefunctions change 
sign upon inversion are known as ‘ungerade’, or ‘u’. The g/u label is 
given as a trailing subscript on the symbol used to designate the elec- 
tronic state. Thus, for example, the ground state of Hp» is of D sym- 
metry while the next two stable states are !£f and ?II, which lie above 
the repulsive ?X7 state. It should be noted that the 2/1 label applies not 
just to diatomic molecules but to all molecules with a centre-of-symmetry 
and that this label is not related to the total angular momentum of the 
system. 

For £X states only, there is an extra symmetry that needs to con- 
sidered. This is the reflection of the wavefunction through any mir- 
ror plane containing the nuclei. Again, the electronic wavefunction V, 
either does or does not change sign with respect to this reflection. 
This parity is denoted by a trailing superscript for X states only, that 
is states are denoted either X* or X . Nearly all states are Z*, which 
are sometimes written simply as X. The one notable exception is the 
oxygen molecule which has a ground state configuration 10210720726; 
30717, 172 which gives the lowest electronic state of symmetry ?X;. 
This reflection symmetry is hard to visualise since it requires two elec- 
trons in an open non-o orbital, such as the 172 orbital in Oz, to geta 
XM state. 

The discussion above has focused on the overall label for the electronic 
state which is of course a many-electron property. Molecular orbitals, 
which are one-electron functions, have symmetry labels similar to those 
used for the states. Thus orbitals are labelled c, 7, ô, ... according to the 
projection, A, of the (single-electron) orbital angular momentum onto 
the molecular axis. As in the atomic case, lowercase letters are used 
for one-electron properties. For homonuclear systems the orbital labels 
are subscripted g or u following the prescription given above. Thus 
the ground state configuration of Hz is 107. It should be noted that 
the c-orbital does not have a + or — superscript. This is because a 
one-electron state can only have symmetric parity and all c-orbitals are 
assumed to be +. 
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9.23 State labels 


Giving an electronic state a symmetry designation such as pra is not 
unique. In fact there is likely to be an infinite number of states with this 
symmetry for a given molecule. In atoms it is common, if not univer- 
sal practice, to precede the spectroscopic term value with n, the principal 
quantum number of the outer electron. For molecules this system would 
still not give a unique set of labels, so a different, rather more ad lioc system 
is used. 

Each electronic state is proceeded by a letter. The following convention 
is used to assign the appropriate letter: 


X labels the ground electronic state; 
A, B, C, ...label states of same spin multiplicity as the ground state; 
a, b, c, ...label states of different spin multiplicity to the ground state. 


In principle, states are labelled alphabetically in ascending energy 
order. In practice there many exceptions. For example, the lowest triplet 
state of Hp is the b ?Z7 with the a ?E7 lying somewhat higher. Sim- 
ilaly there is a higher singlet state, which has a double minimum, 
labelled E, F xs The reasons for these various mislabellings are usually 
historical. 


9.3 Schródinger Equation 


The Schródinger equation for a diatomic molecule with nucleus A, of mass 
Mag and nuclear charge Za, and nucleus B, of mass Mg and nuclear charge 
Zp, and N electrons can be written as 


hi n? 
-Vi -z V$ Vi +Ve— E) (Ra, Rsr)-0, (95 
(i t AVH (RA, Rp, r;) (9.5) 
where the first two terms are the kinetic energy operators for the motions 
of nuclei A and B respectively, the third term gives the kinetic energy oper- 
ator for the electrons, and V, is the potential. The potential is given by the 
various Coulomb interactions within the molecule: 


—(-3 Zi iy m a yl AR). (9.6) 


Areo iz1 “Ai ii Bi i=2 j=1 lij 








Molecular Structure 133 


These terms represent respectively the attraction of the electrons by 
nucleus A, the attraction of the electrons by nucleus B, the electron- 
electron repulsion and the nuclear-nuclear repulsion. The coordinates 
used in this expression are defined in Fig. 9.1. 

Within the Born-Oppenheimer approximation, one writes the wave- 
function as a product of electronic, v, and nuclear, Yn, wavefunctions: 


(Ra, Rs, Fi) = Pelt) bu (RA, Re). (9.7) 


In this case the electronic wavefunction satisfies the simpler Schrédinger 
equation 





h2 N " 
=m; 2N + Ve — Ee | Pelr) =0. (9.8) 
This equation is solved separately for each value of the internuclear sep- 
aration, R. The resulting eigenvalue, E,, is the electronic energy at each R 
and gives the electronic potential V(R) upon which the nuclei move. 

The nuclear wavefunction satisfies the Schródinger equation 


2 2 
(uu VA uu HVO- E) WRR) 09) 
where E, the eigenvalue, is the total energy of the system. 

Solving the Schrödinger equation (9.9) to obtain the total energy E is 
actually not very useful. This is because the equation deals with all the 
motions of the nuclei: vibrations, rotations and the translation of the whole 
system through space. The energy of the translational motion gives a con- 
tinuum and including it in the total energy would obscure the remaining 
discrete spectrum of energy levels. However Eq. (9.9) is a two-body equa- 
tion which can be separated — just as in the H-atom problem — into an 
equation for the translational motion of the centre-of-mass of the system 
plus an equation for the internal motion of the molecule. For a diatomic, 
this latter equation represents the motion of one body in a ‘central’ poten- 
tial, a potential which depends on the distance between the particles but 
not their orientation. The effective mass of this one-body problem is the 
reduced mass, u, where 

Ma Meg 


ERIS 9.10 
PT IMS (9.10) 
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Unlike the H-atom problem (see Sec. 3.3), the reduced mass differs signifi- 
cantly from the mass of the individual nuclei. For example, for a homonu- 
clear system, Ma = Ms and u = Ma, 

If one neglects the translational motion, the Schrödinger equation for 
nuclear motion becomes 





I vt evay- g| d(R) =0, (9.11) 
2u 


where coordinate R = (R, 0, ¢) is a vector giving the internuclear sepa- 
ration R plus the orientation (0, ¢) of the molecular axis relative to the 
laboratory z-axis. 

The vibrational and rotational motion cannot be separated rigorously. 
However such a separation, which gives a good first approximation, can 
be written as 


V&(R) = v(R)v«(0, $). (9.12) 


Considering first the angular or rotational motion, one obtains 


i 1 df. d 1 d 


where the differential, kinetic energy operator is the angular part of the 
Laplacian operator V?, which has been written using spherical polar 
coordinates. 

Equation (9.13) is a Schródinger equation with a zero-potential. It is 
actually the same equation as the angular part of the H-atom problem. 
The solutions are therefore the same ones satisfied by the angular wave- 
functions of the H atom: 


Y,(0, 6) = Yim (0, $), (9.14) 


where Yj, is a spherical harmonic (see Fig. 3.2), and 


h2 
E, = zug +1), (9.15) 


where the rotational angular momentum quantum number J takes inte- 
ger values, J = 0,1,2,.... The projection of J onto the molecular axis is 
denoted MJ. It takes values M; = —J, — J+ 1, ...0,...] — 1, J. There are 
a total of 2J +1 possible values of Mj, which represents the degeneracy 
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of the rotational state. As in atoms, the M}; states split in the presence of a 
magnetic field. 

If the molecule is taken to rotate as a rigid body, which is the simplest 
approximation to molecular rotation, then R can be fixed at some value 
such as Ro. Within this rigid rotor model, the moment of inertia of the 
molecule is Ij = j,R2 and one can write the rotational energy expression as 


2 
E= F-J] += BJ +1), (9.16) 
0 


where Bo is known as the rotational constant of the molecule. Note the 
subscript denotes a vibrational level [see Eq. (9.30)]. In practice, molecules 
are not rigid: atoms move apart as a molecule rotates faster (J increases), 
giving rise to centrifugal distortion terms which can be used to correct the 
expression (9.16). 

It should be noted that if J = 0 then E, = 0. This means that there is no 
rotational zero point energy. 

Measurement of the rotational spectrum of a molecule gives a value 
for Bo. The rotational constant Bo can in turn be used to give an accu- 
rate determination of the average bondlength of the molecule. Laboratory 
spectroscopy, particularly of pure rotational spectra, is the main source of 
accurate data on the geometric structure of molecules. 

The radial equation arising from separating Eq. (9.11) is 


Bog 
B EE SI 47 


where, for simplicity, the dependence on the rotational energy E, has been 
neglected. Within this assumption the energy of the system is given by 


E=E,+E,. (9.18) 


In Eq. (9.17), V(R) is not a simple function of R, so the equation has no 
general algebraic solution. However, for low values of E,, V(R) can be 
approximated by a parabola 


Viso(R) = Vo + KR — RJ? (9.19) 


where the curvature of the potential about its minimum is given by 


_&V 


k E L4 
AR? |k R, 


(9.20) 
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where k is known as the force constant. This parabolic form of the potential 
is the one which gives simple harmonic motion. Figure 9.4 compares a 
harmonic potential with a more realistic one. 

Setting the zero of energy at the minimum of the potential well, Vo = 
V(R = Re) = 0, Eq. (9.17) becomes 


Bg 1 " 
-arge + aA Ro? - Eo] WR =O. (9.21) 


This is the quantum mechanical equation of the harmonic oscillator which 
is solved in most introductory texts on quantum mechanics. See Rae (2002) 
in further reading for example. 

The quantum mechanical solution for energy levels of the harmonic 
oscillator model is 


E, = hw (: af 5) ; (9.22) 


where the vibrational quantum number v takes integer values v — 
0, 1,2, ..., meaning that the vibrational levels supported by a harmonic 
potential are evenly spaced (see Fig. 9.3). In Eq. (9.22) w is the angular 
frequency. It is standard to quote hw in wavenumber units of cm !. The 

















Fig. 9.83. Evenly-spaced energy levels supported by a harmonic potential. Note 
the zero point energy. 
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frequency is related to the constants in Eq. (9.21) by the expression 


w= (=) (9.23) 
H 


The dependence on the reduced mass p of this expression is important 
for studies of molecules with different isotopic composition. This will be 
discussed in Sec. 10.1.1. 

It should be noticed that the vibrating molecule cannot exist at the 
bottom of the potential well. The minimum energy is called the zero point 
energy (zpe). For a diatomic harmonic oscillator the zero point energy 
equals 4ħw. This zero point energy is a direct consequence of Heisenberg's 
Uncertainty Principle relating the position and the momentum. Because 
the potential has a minimum, a state with zero energy would be located 
precisely at Re with no energy and hence no momentum; this would 
violate the uncertainty condition 


AXApx = k ] (9.24) 


The zero point energy also means that the actual dissociation energy, Do, 
is less than D, defined by Eq. (9.2) above. The difference is given by 


Do = D, — zpe. (9.25) 


In other words the amount of energy taken to pull a diatomic molecule 
apart is reduced by the amount of the zero point energy. 


9.4 Fractionation 


The zero point energy of molecules has an important influence on the 
chemistry of the cold interstellar medium (ISM). Consider the simple deu- 
terium substitution reaction: 


H +D > HD +H. (9.26) 


At high temperatures, the ratio between n(HD) and n(H2) will be similar 
to that of n(D) to n(H), where n(X) is the number density of species X. At 
the very low temperatures, below 20K, found in giant molecular clouds in 
the ISM, this is not true due to zero point energy effects. 
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Consider the above example. The harmonic frequency of Hy is 
4396 cm, implying a zpe of 2198cm '!. Mass scaling using Eq. (9.23) 
gives 


w(HD) _ ( um, V 
w(Hy) |— (55) ‘ ne 


if one assumes that the force constant k does not change with isotopic 
substitution, which is true within the Born-Oppenheimer approximation. 
However 


Mu MaMp | 


Assuming Mp ~ 2My and substituting into Eq. (9.27) gives 


1 
My. 3M B o O E 
WHD — WH, (= x zu) = 79 VH = 3807 cm `. (9.28) 


zpe(HD) = z X 3807 = 1903 cm™!. 


A more accurate treatment gives wHp = 3817 cm7! anda zpe of 1909 cm t. 
The energy change for reaction (9.26) is given by 


AE = zpe(HD) — zpe(H;) ~ —289cm™! ~ 420K. 


The lower zero point energy of HD compared to Hz means that at low 
temperature, this reaction, and other similar ones, strongly favour HD for- 
mation over H». This effect, which can lead to between a hundred and 
10000 times more D being observed in molecules than would be expected 
on abundance grounds, is termed 'deuterium fractionation'. It is observed 
in many species, including H?/HD, HCO* /DCO* and Hj /H2D*. 

One can get fractionation of other isotopes such as °C over ^C. How- 
ever, the effect is less extreme as the mass differences are smaller. 





9.5 Vibration-Rotation Energy Levels 


Within the harmonic oscillator approximation all energy levels are evenly 
spaced. Real molecules are not harmonic oscillators. Their potential is elas- 
tic but the repulsion at short bondlengths is stronger than the attraction at 
long bondlengths, so a parabola is only the first approximation to V(R). 
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Fig. 9.4. Comparison of a harmonic potential (dashed line) with a more realistic 
one (full line) for a diatomic molecule with a fairly shallow potential well. 


This approximation is reasonable near R, but becomes increasingly poor 
at larger displacements which correspond to higher energies. 

Figure 9.4 gives a qualitative comparison between the levels in a 
harmonic oscillator potential and those supported by a more realistic 
potential, one which actually leads to dissociation. Potentials which are 
not harmonic are referred to as anharmonic. Unlike the evenly-spaced har- 
monic energy levels, the energy levels get closer together as v increases. 
Also, the harmonic oscillator potential supports an infinite number of 
energy levels, whereas in reality the discrete levels stop at dissociation. 
Thus, for example, H supports 14 vibrational levels in its ground elec- 
tronic state. The precise number of levels depends on the well depth, cur- 
vature and the reduced mass involved. 

A more complete expression for the energy of a vibration—rotation 
state above the bottom of the well for a !E (S = 0, A = 0) diatomic molecule 


is given by 
1\2 
): BeJ(J +1) - were (: * ) 


DPJ +1)? IC ! ;) eme (9.29) 
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where the signs are chosen so that each constant, defined below, is 
positive. In Eq. (9.29), the expansion constants, known as spectroscopic 
constants, are: 


We, the harmonic frequency; 

Be, the rotational constant; 

WeXe, the leading anharmonic correction; 

De, the centrifugal distortion (not the dissociation energy); 
Qe, which gives variation of B, with vibrational state. 


Thus the rotational constant for vibrational state v is given by 


By =B.—a.(v+5) t. (9.30) 
The constants are all usually positive and given in cm !. They are usually 
determined empirically by analysing observed laboratory spectra and are 
different for each electronic state of a given molecule. Huber and Herzberg 
(1979) provide an extensive tabulation of these constants for the observed 
electronic states of diatomic molecules. 

Values of spectroscopic constants for Hz, CO and the molecular ion 
CH" are given in Table 9.2. It is notable that even for hydrogen which is a 
relatively non-rigid system, the energy expressions, at least for low levels 
of excitation, are dominated by the harmonic oscillator, rigid rotor terms. 
In other words, the use of only we and B, gives quite a good representation 
of the actual energy levels. It should also again be noted that the rotational 
energies are much smaller than vibrational energies. 


Table 9.2. Spectroscopic constants for the 
vibration-rotation levels of the ground 
electronic state of molecular hydrogen, car- 
bon monoxide, and the CH* molecular 
ion. Values, which are all in cm ^ 1, are taken 
from Huber and Herzberg (1979). 


Constant H5 CO CH* 


We 4401.21 2170 2740 
WeXe 121.33 13.5 64.0 
Be 60.853 1.93 14.177 
De 0.0471 6x107 0.0014 


Qe 3.062 0.017 0.492 
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9.6 Temperature Effects 


Before considering temperature effects in detail it is useful to give a sum- 
mary of the energy scales found in molecules. Figure 9.5 gives the energies 
in the H2/Hj system. Note that the rotational energy differences are too 
small to be visible on the scale given. 

As molecules have many closely-spaced energy levels, they act as a 
useful thermometer in a variety of environments. To understand this it 
is necessary to consider the population of the various levels in a typical 
molecule as a function of temperature. 


9.6.1 Rotational state populations 


Consider the Boltzmann distribution for a diatomic molecule rotating as a 
rigid rotor. For simplicity take the J — 0 level as the zero of energy. Define 


Ht + H(1s) 





H5* 


14 + 


12 + 
I(H*) 
(01 [ræ 





H(1s)  H(1s) 


Energy/ eV 


D.(H2) 














H2 





-— zpe 


Fig. 9.5. Energy scales in hydrogen molecule. Zero is taken as the H5(v = 0) 
ground state. I represents the ionisation potential and D the dissociation energy. 
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P(J) as the probability of finding the molecule in level J. To avoid having 
to evaluate the partition function, let us consider the ratio of P(J) to the 
probability of finding the molecule in its rotational ground state, P(0): 


PU 8 -Ej 
BO a a 


where 2; is the degeneracy factor which is (2J + 1) for a heteronuclear 
diatomic. Note that for homonuclear diatomics there is an additional fac- 
tor caused by nuclear spin effects and ortho/para states, but this is beyond 
the scope of the present book. The energy levels of a rigid rotor are 
E; = BJ(J + 1), therefore 

PJ) -BJU +1) 

PO > (2J + 1)exp Em (9.32) 
For a typical diatomic, B ~ 1 cm '!, and in these units it useful to remember 
that Boltzmann’s constant, k, equals 0.695 cm~!-K~!. For a temperature of 
T — 300K, these values give 


= -1x2 _ 2,-0.0096 _ 
nj = See ( ran) = 3e = 2.97, 
which means that there are more molecules in the J = 1 level than in the 
J = 0 level. In fact P(0) ~ P(28) for this example. 

For a diatomic rigid rotor with rotational constant B, the rotational 
level with the largest population is given by the approximate relationship: 


kT\? 1 
= (5) -5 (9.33) 


At room temperature many rotational levels of a molecule are occu- 
pied. Even at cold interstellar medium temperatures of 10-20K, most 
molecules have several levels occupied. Indeed, one of the earliest 
molecules observed in the interstellar medium, CN, was found in the 
early 1940s. The Nobel Prize-winning molecular spectroscopist Gerhard 
Herzberg (1904-1999) commented on this work that ‘from the inten- 
sity ratio of the lines... a rotational temperature of 2.3K follows, which 
has of course a very restricted meaning’. It is not as restricted as he 
thought, as Arno Penzias and Robert Wilson received a Nobel Prize 
for their 1963 demonstration of the existence of the cosmic microwave 
background radiation with a temperature of 2.7K. The CN rotational 
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temperature measurement was actually the first observational evidence of 
the cosmic microwave background radiation. The spectrum of CN is now 
used as one means of measuring this background temperature in different 
locations. 

The rotational level populations of molecules such as CO and CN can 
be used to determine temperature in a particular region. However in many 
diffuse environments, the rotational levels of heteronuclear diatomics, or 
indeed any molecule with a permanent dipole moment, are usually not 
thermally occupied as the radiative lifetime, given by Ay is less than the 
mean time between collisions. That means the molecules exist in an envi- 
ronment which is below their critical density. The molecular emissions in 
planetary nebula NGC 7027 (Fig. 7.6) are like this. 


9.6.2 Vibrational state populations 


Vibrational energy separations are larger than rotational ones and are 
therefore only sensitive to higher temperatures. To illustrate this, consider 
the probability, P(v), of finding a diatomic molecule in vibrational level v. 
For simplicity, assume that the diatomic is a harmonic oscillator and set 
the zero of energy at the v = 0 level. Again taking the ratio of the vth state 
to the vibrational ground state gives 


PRO)? - —hwv 
PO = exp (== ) : (9.34) 


It should be noted that there are no degeneracy factors associated with the 
vibrational motion. 

Diatomic vibrational frequencies range from 500 cm“! for a heavy sys- 
tem to 4000 cm ! for H2. 1000 cm '! is a typical value so taking this for a 
temperature of 300K gives 


P(1) —1000 
P(0) ( 


P(0) 0.695 x x) ey 
which means that less than 1% of the molecules are in their v = 1 vibra- 
tional level. 

So at low temperatures, molecules are vibrationally cold with only 
the vibrational ground state significantly occupied. However, hot environ- 
ments can lead to thermal emissions from vibrationally excited molecules. 
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This is quite often found for Hz such as in Figs. 7.6 and 10.10 for 
example. 


9.6.3 Electronic state populations 


For most molecules the dissociation energy is less than the lowest elec- 
tronic excitation energy. This means that electronically-excited states are 
not thermally occupied as the molecules are destroyed instead. Emission 
from electronically-excited states are still observed but these emissions 
are usually not thermal. For example, emissions are observed in plane- 
tary aurora, from electronically-excited H» in Jupiter and Saturn, and from 
other diatomics on Earth. These emissions are the direct consequence of 
excitation by collisions with electrons. 

Astronomically there is one important exception to this rule. The 
atmospheres of cool stars, such as M-dwarfs, which have temperatures in 
the 2200-3700K range, contain molecules. In these stars the most impor- 
tant diatomic molecules are H) and CO. However, a number of molecules, 
such as TiO, FeH, ZrO and VO, are formed which contain a transition 
metal atom. These molecules have open shells which result in many low- 
lying electronic states. These states are thermally occupied in the atmo- 
spheres of cool stars. These molecules have dense spectra which are often 
a significant source of opacity as they absorb radiation from the interior 
of the star over a wide range of wavelengths. The spectra of the molecules 
are very complicated. They are often poorly understood and in many cases 
remain completely uncharacterised. Figure 6.9 gives an example which 
shows that the absorptions in small, cool stars are dominated by molecular 
absorptions. The spectra of these open shell species will not be considered 
further here. 


Problems 


9.1 A sample of a diatomic molecule with rotational constant B, is in ther- 
modynamic equilibrium at temperature T. Show that the ratio of the 
number of molecules in rotational level J to the number in rotational 
level zero is a maximum for the level with 


1 
kT\? 1 
1- (55) E. 


where k is Boltzmann's constant, which has a value 0.695 cm! - K- 1. 


9.2 


9.3 
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For the molecule CO, B, = 1.93cm !. Which are the most-occupied 
rotational states of CO at 


(a) the ISM temperature of T — 20K, 
(b) room temperature, T = 300K, 
(c) the temperature of a typical M-dwarf star of 3000 K? 


Assuming a harmonic frequency w(CO) = 2170 cm~}, what proportion 
of the molecules are in the v = 1 vibrational state compared to v = 0 at 
each of these temperatures? 

The fundamental vibrational frequency of the "CH* molecular ion 
is w = 2075.5cm . Estimate the zero point energy of "CH" and, 
assuming integer values for the atomic masses, estimate the zero 
point energy of the less-abundant PCH* ion. Why might one want to 
observe spectra of the less-abundant PCH* species instead of "CH? 
CH’ ions are found in giant molecular clouds. How would you expect 
the ratio n( DCH^) : n("CH^) to compare to the ratio n(?C) : n(?C), 
where n(X) represents the number density of species X? 


CHAPTER TEN 


MOLECULAR SPECTRA 


The astronomical spectra of molecules give rise to three distinct types 
of transitions, which have to be considered separately. 

Pure rotational transitions lie at long wavelengths ranging from radio 
frequencies for heavy polyatomic molecules, through the microwave, to 
the far-infrared for light hydrogen containing diatomics. 

Vibrational transitions are important at mid-infrared wavelengths. 

Electronic transitions lie at similar wavelengths to the allowed transi- 
tions of neutral atoms: the visible and ultraviolet. The wavelengths quoted 
are only typical ones. In each case the precise spectral region depends 
on the spectroscopic constants of the molecule in question. Each of the 
transition types, and their associated selection rules are considered in 
turn below. 


10.1 Selection Rules: Pure Rotational Transitions 


For a diatomic to undergo a pure rotation transition, it must have a 
permanent dipole moment, u. Actually for any molecule to have a dipole- 
allowed rotational spectrum it must have an asymmetric charge dis- 
tribution giving a permanent dipole moment. Heteronuclear diatomics 
possess a permanent dipole moment but homonuclears, such as H2, do 
not. Molecular hydrogen therefore does not have a dipole-allowed rota- 
tional spectrum. 

The strength of a rotational transition depends on p°, therefore 
molecules with large permanent dipoles have intense transitions. This 
means that molecules such as sodium chloride, Na* - CI-, which have 
very large dipole moments as a result of the almost complete charge 
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separations, have particularly intense rotational transitions. This has 
greatly aided the astronomical detection of such species even at very small 
concentrations or column densities (see Fig. 10.1). 

For diatomic molecules in € symmetry electronic states, dipole- 
allowed transitions obey the following selection rule: 


AJ 2343. (10.1) 





Within the rigid rotor approximation it is therefore straightforward to 
obtain the energy of the rotational transitions: 


AErot = Bo o + 1) m TU T 1)| , (10.2) 


where the standard convention for molecular spectra has been followed — 
that’ is used to denote the upper state and " means the lower state. 
As the selection rule gives J’ = J” + 1, 


AE rot = 2BoJ' . (10.3) 


This means that the rigid rotor approximation leads directly to very reg- 
ular pure rotational spectra where the transitions are regularly spaced by 
2B. Thus the transition J'- J” = 1—0 is at 2B, 2-1 at 4B, 3-2 at 6B, and 
so forth. It should be noted that, in contrast to atomic transitions, molecu- 
lar transitions are generally denoted upper — lower. Sometimes an arrow is 
inserted so that J’ — J” denotes absorption and J’ — J” denotes emission. 

Centrifugal distortion leads to a slight narrowing of the gap between 
neighbouring transitions as J increases, however the regular progression 
is still clearly seen. Figure 10.2 gives a sample pure rotational spectrum 
of CO recorded in the laboratory. The regular spacing of the spectrum 
for both the most-abundant species, 1*C!°O, and isotopically-substituted 
species is easily seen. Astronomical pure rotational spectra are, of course, 
also regularly spaced when plotted using a frequency scale. However it 
unusual for a single astronomical spectrum in the radio or far-infrared 
to span a wide range. The more usual representation for series of astro- 
nomical rotational lines are given in Fig. 10.1 for NaCl and Fig. 10.3 
for CO. 

Carbon monoxide, CO, is a particularly important species for astro- 
nomical observations. CO is the most stable diatomic molecule. It has 
a dissociation energy Do of 11.1eV, which is more than double the Do 
value found for most other diatomics. As a result, in astronomical envi- 
ronments where molecules form, C and O usually combine to form CO, 
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Fig. 10.1. Five spectra observed towards IRC+10216 with the IRAM 30m 
telescope showing the J = 7-6 to J = 12-11 rotational transitions of NaCl. 
[Reproduced from J. Cernicharo and M. Guélin, Astron. Astrophys. 183, L10 (1987).] 
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Fig. 10.2. Laboratory far-infrared absorption spectrum of CO. [Reproduced from 
P.F. Bernath, Spectra of Atoms and Molecules (Oxford University Press 1995).] 


which is very stable and long-lived. Any surplus C or O is then available 
for other chemistries. CO is thus the second most abundant molecule in 
the Universe after H2. 

For CO, the rotational constant of the ground vibrational state, Bo, 
is 1.93cm !. The wavelengths of the first few rotational transitions are 
1-0 at A = 2.60mm, 2-1 at 1.30 mm, and 3-2 at 0.87 mm. Examples of 
low J transitions are shown in Fig. 10.3. Much higher rotational levels can 
also be observed in warm objects such as the planetary nebula NGC 7027 
(see Fig. 7.6). 

The J = 1-0 transition of CO is the second most important spectral 
line in radio astronomy after the hydrogen 21 cm line. CO is widely dis- 
tributed in the interstellar medium and maps of the CO J = 1-0 transition 
are a standard tool for investigating the ISM. One reason for this is that 
cold H» is very difficult to observe directly because its pure rotational tran- 
sitions are not only very weak but lie in the near-infrared where ground- 
based observations are not possible. The abundance of CO is therefore 
often used to estimate the total amount of molecular gas present in a 
given environment. It is generally assumed that the number density of 
CO, n(CO), is approximately 10~* n(H;). 
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Fig. 10.3. Comparison of profiles from CO 2-1 and 1-0 absorption against 
extragalactic sources observed with the Plateau de Bure interferometer. The 2-1 
spectra are displaced vertically by +1 for clarity. CO temperatures can be obtained 
by comparing the two absorption profiles. [Reproduced from R. Lucas and H.S. 
Liszt, in Molecules in Astrophysics: Probes and Processes, ed. E.F. van Dishoeck 
(Kluwer, 1997).] 


If, as often happens, the CO 1-0 line is optically thick, one can use 
higher transitions such as the CO 2- 1 line instead. These higher lines, 
at millimetre wavelengths, can be observed using telescopes such as the 
James Clerk Maxwell Telescope on Mauna Kea, Hawaii. Figure 10.3 gives 
an example of monitoring both the 1-0 and 2-1 lines together, which 
can be used for temperature determination. Another option to avoid the 
effects of optical thickness is to observe isotopically-substituted CO, which 
is present with much lower densities and whose transitions are therefore 
much less optically thick. As discussed below, isotopic substitution has 
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a particularly marked effect on molecular spectra, making such observa- 
tions fairly straightforward. 

It should be noted that CO actually has a rather small permanent 
dipole of only 0.12 Debye. This means that even though the rotational 
spectrum is driven by this electric dipole, the Einstein A coefficient is 
actually fairly small for the pure rotational transitions of CO. For exam- 
ple Ai_9 = 7 x 10 ?s^!. Note that the factor v? in the definition of the 
Einstein A coefficient [see Eq. (2.4)], means that the radiative lifetime of 
rotationally-excited states is in any case significantly longer than those 
generally found for electronically-excited states of atoms or molecules 
whose transitions lie at much higher frequencies. 


10.1.1 Isotope effects 


Isotopic substitution means the replacement of an atom in a molecule by 
another isotope of that atom. Thus in CO, whose normal form is C160, 
one can replace *C by PC, giving ?C!6O. This change leaves the chem- 
istry unaltered, except for the fractionation effects discussed in Sec. 9.4, 
but leads to changes in the spectrum due to mass effects. 

The rotational constant for CO is 


n? McMo 


=-——, So 10.4 
2ure K Mc + Mo ( 


Bo 


Replacing ?C by PC leads to an increase in the reduced mass u and 
hence an increase in the moment of inertia, I, assuming that the effective 
bondlength rp is unchanged upon substitution. This leads to decreased 
rotational constant, B, for C160 and hence a smaller separation between 
the rotational transitions. This change in transition frequency is significant 
and can easily be resolved even with moderate resolution. Thus the 1-0 
transition is at 2.60 mm for !2C!°O, at 2.67 mm for CO, and at 2.72 mm 
for 8C!°O. 

The spectrum of ?C!5O can easily be observed. See Fig. 10.3 for exam- 
ple. In fact, astronomical spectra of all versions of isotopically-substituted 
CO have been recorded. 


10.1.2 Rotational spectra of other molecules 


There are presently over 120 different molecules whose spectra have been 
observed in the interstellar medium; this number is growing steadily by 
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about 3 species per year. Most of these molecules have been detected 
by the observation of rotational transitions at radio frequencies. Radio 
spectra of giant molecular clouds are rich in features, particularly if one 
probes deeply by using long integration times and sensitive detectors. 
Figures 10.4 and 10.5 show a survey and a close-up spectrum recorded 
looking at the nearby massive star-forming region in Orion. Such studies 
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Fig. 10.4. Radio spectrum spanning the 138.3 to 150.7 GHz region obtained 
towards Orion-KL. Strong lines are designated their molecular identifications. 
‘U’ marks unidentified lines. [Reproduced from C.W. Lee, S.-H. Cho and S.-M. Lee, 
Astrophys. J. 551, 333 (2001).] 
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Fig. 10.5. Detailed radio spectrum spanning obtained towards Orion-KL. Lines 
are designated by their molecular identifications. [Reproduced from C.W. Lee, 
S.-H. Cho and S.-M. Lee, Astrophys. J., 551, 333 (2001).] 
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highlight the large difference in intensity scales of the various transitions 
that can be observed. That there are molecules still waiting to be found 
is beyond doubt because there are many unidentified spectral features. 
See for example the strong lines marked ‘U’ in Fig. 10.4. Higher-resolution 
studies reveal many, many more such features. Most often such features 
are assigned as a result of detective work based on spectroscopic studies 
performed in the laboratory. 

Interstellar medium molecules are most often observed in emission, 
but can also be sometimes found in absorption. Indeed, there are a num- 
ber of cases where molecules have been observed in absorption against the 
cosmic microwave background; these observations imply that the absorb- 
ing molecule is at a temperature below 2.7K as measured by its rota- 
tional population. As in the atomic case, observation of high-resolution 
line profiles can yield further information: Doppler shifts can give rela- 
tive motion, Doppler broadening can give the (translation) temperature, 
and so forth. Figure 10.6 gives an example of self-absorption by the carbon 
monosuphide molecule, CS. Light from further away is absorbed by the 
molecules nearer to us. In this case the nearer molecules are cooler as their 
absorption profiles are significantly narrower. These particular observa- 
tions are interpreted as evidence for the collapse of the dark molecular 
cloud Barnard 335 as a precursor to star formation. 


Temperature (Ty; K) 





Velocity (km s") 


Fig. 10.6. Rotational transitions of the CS molecule observed towards the cen- 
tre of massive dark cloud Barnard 335. The dashed lines are from a model fit. 
The J = 2-1 and J = 3-2 show self-absorbed profiles with stronger blue peaks. 
This behaviour is typical of in-falling gas. [Reproduced from S. Zhou, Molecules in 
Astrophysics: Probes and Processes, ed. E.F. van Dishoeck (Kluwer, 1997).] 
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10.1.3 Rotational spectra of molecular hydrogen 


Hz is a homonuclear diatomic and therefore does not have a permanent 
electric dipole moment or an allowed rotational spectrum. It should be 
noted that the molecule HD does have a dipole but only a very small one 
of 0.0008 Debye. This leads to an Einstein A coefficient of only 3 x 10 5s! 
for the lowest J = 1-0 rotational transition. 

However the abundance of Hz in many locations in the Universe 
means that electric quadrupole transitions can be observed. These transi- 
tions are weak, like the forbidden quadrupole transitions found in atomic 
spectra. For diatomic molecules, the selection rule for pure rotational spec- 
tra driven electric quadrupoles are 





AJ =+2. (10.5) 


The H; molecule is very light so that By = 62cm™!, meaning that the rota- 
tional transitions of H; lie at infrared wavelengths. For example, the low- 
est 2-0 line is at à = 28 um. These transitions are about a billion times 
weaker than dipole-allowed rotational transitions. The Einstein A value 
for the 2-0 transition is very small, 3 x 10 ! s^, which corresponds to a 
radiative lifetime against spontaneous emission of about 1000 years. 

H» quadrupole rotational transitions can be seen in the warm ISM in 
emission (see Fig. 10.7). Such transitions have been extensively studied in 
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Fig. 10.7. Quadrupole rotational emission spectra of molecular hydrogen 
towards the warm photon-dominated region (PDR) $140 recorded using the 
Infrared Space Observatory (ISO). The observations probe the warm, dense 
molecular gas adjacent to the $140 HII region. The flux density (vertical scale) 
is given in Jansky (— 10 a: ergs 1.cm-?. Hz- 1). Dashed lines give the laboratory 
wavelengths for the transition indicated. [Adapted from R. Timmermann et al., 
Astron. Astrophys. 315, L281 (1996).] 
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Orion; indeed higher rotational transitions have been observed in Orion, 
up to S(17) or J = 19-17, than those detected in the laboratory. In the lab- 
oratory, getting a sufficient pathlength of light through a sample of H; to 
observe these very weak transitions is a major challenge. 


10.2 Vibrational Transitions 


Within the harmonic oscillator approximation, there is a simple and rigor- 
ous electric dipole selection rule: 


Av = +1. (10.6) 





This leads directly to 





AE, = hoo +14 5) «(s ' 5) =w, (10.7) 


where w is known as the fundamental frequency. 

For anharmonic molecules any change in Av is allowed in principle, 
but in practice, Av = +1 always leads to much stronger transitions. The 
intensity of individual transitions falls off rapidly with increasing Av. 
Vibrational transitions which change v by more than one quantum are gen- 
erally called overtones. 

Vibrational transitions are usually accompanied by a change in rota- 
tional state. For electric dipole transitions within a X electronic state, the 
additional selection rule on rotational motion is 





AJ 2343. (10.8) 





For non-X states and polyatomic molecules, transitions with AJ — 0 can 
also occur. To get dipole-allowed vibrational excitation, one must have a 
change in the dipole moment 


du 
do #0, (10.9) 


where Q is the vibrational coordinate undergoing excitation. 

For homonuclear diatomics, 4 = 0 for all internuclear separations, 
so vibrational transitions are all dipole-forbidden. This is not true, how- 
ever, for symmetric polyatomic linear molecules such as carbon dioxide. 
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CO, does not have a dipole moment in its equilibrium geometry but 
distortions from this can give an instantaneous dipole moment. Hence 
some CO; vibrational transitions are dipole-allowed. It is these vibrational 
transitions which contribute to the greenhouse effect both here on earth 
and on Venus, where the super-abundance of CO; in the atmosphere is 
responsible for maintaining a huge greenhouse effect. 

For heteronuclear diatomics, such as CO, the derivative of the dipole 
with respect to bondlength, ae is only zero by accident. So heteronuclear 
diatomic molecules all have allowed vibrational spectra. 


10.2.1 Structure of the spectrum 


For simplicity the discussion will be confined to rotation-vibration spec- 
tra in E symmetry electronic states. Consider the so-called fundamental 
transition given by o' = 1 — v” = 0 centred on the the fundamental fre- 
quency, w. The rotational selection rules are AJ = +1, meaning that the 
vibrational band has two branches. Assuming the rigid rotor approxima- 
tion this gives: 





R-branch: AJ = +1, ie. J/ = J” -1 





AEvg = hw + J'(J --1)B4 — JJ" +1)Bo ~ hw --2B(J" +1) (10.10) 





for J” = 0,1,2,.... It has been assumed that B4 ~ Bo, which is generally 
a good approximation, and is about as reliable as neglecting the effects of 
centrifugal distortion. 


P-branch: AJ = -1,ie.]/—J"—1 


AEqg = hu + J'(J!' +1)Bi — "(" +1)Bo ~ hw — 2B]" (10.11) 





for J” = 1,2,3,... but not J” = 0 as it would imply J’ = —1. In standard 
notation the transitions for the two branches are labelled R(J") and P(J”), 
where J” is the rotational quantum number associated with the lower 
vibrational state. 

This series of transitions give a characteristic spectral pattern as 
depicted in Fig. 10.8. The rotational selection rules lead to a series of lines 
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Fig. 10.8. Stick spectrum showing the structure of a vibrational band for a 
diatomic molecule with a £ symmetry electronic state. 


which are evenly spaced, 2B apart. The exception is at the centre of the 
band where no line is observed. 

Figure 10.9 shows a CO absorption spectrum recorded looking 
towards the massive young stellar object AFGL 4176, which is embed- 
ded in a dense molecular cloud. The symmetric structure of the P- and 
R-branches, and the band centre, which can be pinpointed by the missing 
line, are all clearly visible. The shorter wavelength feature due to solid 
carbon dioxide shows how different these largely structureless condensed 
phase spectra are. The 4.7 um region is heavily obscured by the earth’s 
atmosphere. This spectrum was recorded by a satellite, the Infrared Space 
Observatory (ISO). 

Besides the structure of the lines, the vibration-rotation spectra of 
diatomics also have a characteristic shape. The intensity of individual lines 
is directly proportional to the rotational population in the initial state of 
the molecule. This means that this intensity distribution can be used to 
measure the (rotational) temperature of the molecule. The temperature can 
then be estimated from a spectrum using formula (9.33), although a more 
reliable result can usually be obtained by fitting the whole spectrum to a 
temperature-dependent rotational population. 

For diatomic molecules not in a È electronic state, that is, ones with 
A > 0, and for polyatomic molecules, the rotational selection rule (10.8) is 
generalised to 


AJ =0,+1, notJ=000. (10.12) 





158 Astronomical Spectroscopy 








> 
+ 
a 
& 
Ya 
z 
o 
2 
E Gaseous CO 
[2 
6 
Solid CO; 
AFGL 4176 
ISO-SWS 
4 l à. 
42 4.4 46 4B 5 5.2 


Wavelength (um) 


Fig. 10.9. Spectra obtained towards the massive young stellar object AFGL 4176 
using the Infrared Space Observatory (ISO). The strong, largely structureless band 
absorption at 4.27 um is due to solid CO», whereas the characteristic vibration- 
rotation P- and R-branch structure between 4.4 and 4.9 um is due to the pres- 
ence of warm, gaseous CO along the line of sight. [Reproduced from E.F. van 
Dishoeck, in The Molecular Astrophysics of Stars and Galaxies, eds. T.W. Hartquist 
and D.A. Williams (Clarendon Press, Oxford, 1998).] 


The extra AJ — 0 transitions form what is known as a Q-branch. Within the 
rigid rotor approximation, all Q-branch transitions in a vibration-rotation 
spectrum lie at the same frequency, giving the spectrum a characteristic 
sharp peak at the band origin. Effects beyond the rigid rotor model can 
introduce slight shifts in these transitions but it is common for the indi- 
vidual Q-branch transitions to not be resolved even with high resolution. 

Vibrational spectra can be seen in emission from hotter regions of the 
ISM and of planetary atmospheres particularly from auroral regions. They 
are seen in absorption in the atmospheres of cool stars and brown dwarfs. 
Absorption spectra of cold ISM molecules can be seen against the light of a 
suitable, more distant infrared source. Such observations, which are often 
necessary to identify interstellar molecules which do not have a perma- 
nent dipole, can be difficult. 
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10.2.2 Isotope effects 


Isotopic substitution alters vibrational frequencies as well as rota- 
tional ones. However, whereas the rotational transition frequency shifts 
according to the inverse of the reduced mass u (see Sec. 10.1.1), the vibra- 
tional band origin shifts as p i. This can be seen from Eq. (9.23) which is 
based on the harmonic approximation. A further assumption, valid within 
the Born-Oppenheimer approximation, is that the force constant, k, does 
not change with isotopic substitution. Thus for C!^O the band origin is 
shifted to 2096 cm! from 2143 cm '! found for C'°O. Such shifts are suf- 
ficient to be detectable with a spectrograph of moderate resolution. 


10.2.3 Hydrogen molecule vibrational spectra 


Like its rotational spectrum, Hz only undergoes very weak electric 
quadrupole-induced transitions for its vibrational spectrum. For these 
transitions, vibrational state changes with Av = +1 are still dominant, 
although transitions with Av > 1 are observed astronomically. The rota- 
tional selection rules for these quadrupole transitions are: 





AJ =0,+2, not] 20.0. (10.13) 





This leads to a new set of labels for the different branches. Transitions 
with AJ — 2 form the S-branch; AJ — 0 form the Q-branch; AJ — —2 form 
the O-branch. Transitions are labelled S(J”), Q(J") and O(J”), where J” 
is the rotational quantum number of the lower rotational state involved 
in the transition. 

Despite the extreme weakness of the transitions — they are approxi- 
mately 10? times weaker than electric dipole-allowed vibration-rotation 
transitions — emissions from vibrationally excited H» are bright from 
a number of astronomical sources. These include planetary atmospheres 
(Fig. 10.10), planetary and reflection nebulae (Fig. 10.11) and hot regions 
of the ISM. H» emissions can be particularly bright from merging galaxies, 
which contain huge quantities of hot hydrogen molecules. 

Both Figs. 10.10 and 10.11 show effects due to the earth's atmosphere, 
the so-called telluric effects. This is a common problem with observing in 
the infrared. In particular, the gap between 1.8 and 2.0 um in Fig. 10.11 due 
to atmospheric effects separates the H-band (1.4-1.8 um) from the K-band 
(2.0—2.4 um). 
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Fig. 10.10. Continuum-subtracted, infrared K-band spectrum of Uranus recorded 
using the UK Infrared Telescope (UKIRT) by Trafton and co-workers [see 
Astrophys. J. 405, 761 (1993)]. The Hz quadrupole emission lines are prominent, 
especially in the Q-branch. Weak emission features from the overtone band of the 
hydrogenic molecular ion Hr can also be seen. The structure in the spectrum near 


2.01 um and long-wards of 2.45 um is caused by telluric features. 
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Fig. 10.11. Infrared spectra of reflection nebula NGC 7023. The wavelengths of 
molecular hydrogen lines have been labelled above the spectrum. The intensity 
has been normalised to the peak 1-0 S(1) intensity. [Reproduced from P. Martini, 
K. Sellgren and J.L. Hora, Astrophys. J. 484, 296 (1997).] 
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10.3 Electronic Transitions 


Electronic transitions involve a jump from one potential energy curve to 
another curve (see Fig. 10.12). The most common situation is that the equi- 
librium bondlength of the excited electronic state, Rj, is bigger than that of 
the ground state, R}, as electronic excitation usually weakens the bonding 
in a molecule. The potential energy curves of Hp (see Fig. 9.2), illustrate 
this behaviour. 

Electronic spectra are significantly more complicated than the pure 
rotational and vibration-rotation spectra considered above. This is 
because it is necessary to consider changes in the electronic, vibrational 
and rotational state of the molecule for each possible change in electronic 
state. 


10.3.1 Selection rules 


Nearly all electronically-excited states of molecules contain enough energy 
to dissociate the molecule. For this reason it is usually only necessary 
to consider electronic transitions which obey electric dipole selection 
rules. Table 10.1 gives the selection rules for electric dipole transitions 


Energy 








Fig. 10.12. Electronic transitions involve a jump from one potential to another. 
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Table 10.1. Selection rules for spectra of diatomic molecules 
undergoing allowed electric dipole transitions. 











Rotations AJ =+1 for A= 0-0, 
AJ —0,-X1  not/-— 0-0, for other AA. 
Vibrations Av any. 
Spin AS — 0. 
Orbital AA — 0, +1. 
X states xbox, 
Le Se: 
Symmetry geu Homonuclear molecules only. 


for diatomic molecules. In principle these selection rules apply to all 
transitions of a diatomic molecule, but in practice the table is of most 
use for considering electronic transitions as these have the most compli- 
cated selection rules. The simplified versions of the selection rules given 
above for pure rotation and vibration-rotation transitions are sufficient for 
these cases. 

For all molecules, allowed electric dipole transitions obey the spin con- 
serving selection rule 


AS — 0. (10.14) 


For diatomic molecules the projection of the total electron orbital angular 
momentum along the diatomic axis, A, obeys the rule: 





AA 20, +1, (10.15) 


which means that Z—-X, X-II and II- A transitions are allowed but a X-A 
transition is forbidden. 

Molecules with a centre-of-symmetry, which include homonuclear 
diatomics, must change symmetry when undergoing an electric dipole 
transition. Thus transitions must link g and u electronic states; transitions 
which are u + u and g — gare forbidden. This rule is the molecular equiv- 
alent of the Laporte rule encountered in atomic spectra. Its application 
explains why there are no dipole-allowed rotation or vibration-rotation 
spectra for homonuclear diatomics since such transitions all occur within 
a single electronic state and therefore leave the symmetry of the state 
unchanged. 


Molecular Spectra 163 





tie] 


in 


erg s'em 
o 


E" 
E- 
E 
‘| =g 

Fs 
d 


Flux (10 * 











M Mem eee lip Pri gym NP 


1003 1020 1040 1080 

















Tlux (10 erg s 


E 
ə 


1105 Hr 1140 


Fig. 10.13. Spectrum towards the star DI 1388 recorded with the Far-Ultraviolet 
Space Explorer (FUSE). The interstellar atomic and molecular lines arising in the 
Galaxy (MW) and the Magellanic Bridge (MB), a region between the large and 
small Magellanic clouds, are identified. [Reproduced from N. Lehner, Astrophys. J. 
578, 126, (2002).] 


Figures 10.13 and 10.14 show ultraviolet spectra of molecular hydro- 
gen recorded by observing starlight passing through the interstellar 
medium. Prominent electronic transitions shown for H, are the 


Werner band: C !I,- X p» at about 1100 Á; 
Lyman band: B !E; -X !Ef at about 1010 À. 


These transitions are described as bands and do not occur at a precise 
wavelength because within each electronic transition there are a series of 
vibrational and rotational transitions. These are considered in turn below. 


10.3.2 Vibrational selection rules 


There are no rigorous selection rules which govern the change in vibra- 
tional quantum numbers during an electronic transition. The chance of a 
particular vibrational transition occurring depends on the squared overlap 
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Fig. 10.14. Spectrum recorded towards O-star ¢ Oph using the satellite 
Copernicus. [Reproduced from L. Spitzer Jr. and E.B. Jenkins, Ann. Rev. Astron. 
Astrophys. 13, 133 (1975).] 


of the vibrational wavefunctions from the two potentials involved: 


| f polo)" molar) (10.16) 


This factor, which takes a value between zero and one, is called a 
Franck-Condon factor. It is derived following the Franck-Condon prin- 
ciple, which can be thought of as an extension of the Born-Oppenheimer 
approximation as it assumes that electronic transitions occur too rapidly 
for the nuclei to move during them. See Bransden and Joachain (2003) in 
further reading for a full discussion of the Franck-Condon principle. 

It is common practice to denote vibrational bands within an electronic 
transition by (v’, v”). So (4, 0) means a transition between vibrational state 
v' = 4 in the upper state and v" = 0 in the lower state. 


10.3.3 Rotational selection rules 


The rotational selection rule for an electronic transition is 





AJ =0,+1, notJ20—0, (10.17) 


except that AJ = 0 transitions are not allowed for transitions between 
Y. electronic states. These selection rules lead to P-branch (J' - ]" = —1), 
Q-branch (J’ = J”) if allowed, and R-branch (J' - J” = +1) transitions. How- 
ever, unlike vibrational spectra, these branches do not display a regular 
structure. 
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Table 10.2. Spectroscopic constants, in cm 1, for the electronic states of 
molecular hydrogen involved in the Lyman and Werner bands. 


State Te We WeXe Be Qe De 


CUTE, 100089.8 2443.77 69.524 31.362 1.644 0.0223 
Blxt 91700.0 1358.09 20.888 20.015 1.185 0.0163 
X ps 0 4401.21 121.33 60.853 3.062 0.0471 





The rotational constant associated with each electronic state involved 
in a transition is different, i.e. B’ A B"; in fact they are often substantially 
different. This is because although the reduced mass of the system stays 
constant during a transition, the equilibrium bondlength, Re, does not. The 
most common situation is that R} > R7 so that B’ < B". Table 10.2 gives 
the molecular constants which can be used to determine the wavelength 
of individual transitions in the Werner and Lyman bands. This shows that 
B, for the X Qu ground state is substantially bigger than B, for the two 
electronically-excited states. 

The difference in rotational constants means that the rotational 'fine 
structure' lines are not evenly spaced. Furthermore the transitions often 
pile up on top of each other leading to a band head where many rotational 
transitions occur very close together. Band heads occur at the red end of 
the band, if B' « B", which is the most common situation, or at the blue 
end of the band, if B’ > B". 


10.3.4 Transition frequencies 


The energy difference, AE, for an electronic transition contains contribu- 
tions from each of the electronic, vibrational and rotational energies. The 
difference is most conveniently written 


AE — T, + Elg — Eve, (10.18) 


where T, is the adiabatic electronic excitation which is defined as the 
energy difference between the minima in the potential energy curves of 
the two states involved (see Fig. 10.12). By convention the minimum of 
the ground state curve is taken as the energy zero. The vibration-rotation 
energy levels within each curve, Eyr, are usually obtained in terms of stan- 
dard molecular constants. Expressions for them are given in Eq. (9.29). 
The relevant constants for the Lyman and Werner bands of H» are given 
in Table 10.2; a general tabulation of spectroscopic constants for diatomic 
molecules can be found in Huber and Herzberg (1979). 
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Worked Example: A cold interstellar column of H; absorbs light in the 
B'Zj-X!Xi Lyman band. This band does not have a Q-branch. Assum- 
ing that the H} molecule is entirely in its vibrational ground state and 
J” = Qand 1 rotational states, what wavenumber of light will be absorbed to 
populate vibrational states with v' — 0 and 1 in the upper electronic state? 


Solution: Denote the energy levels E x(v", J”) and Eg(v', J’) for the ground 
and excited electronic states respectively. Then, using spectroscopic con- 
stants for the appropriate electronic state taken from Table 10.2, 


i PA 


w 


wW, X 








Ex(0, 0) = DA = 1 e —21703cm !, 
Ex(0, 1) = Ex(0, 0) + 2B/ — 4D” — a! = 2288.8cm ! , 

/ LEVIA 
Est 0j e T. T = Lh = 92690.1 cm}, 
Eg(0, 1) = Eg(0, 0) + 2B} — 4D! — 2a!, = 92726.5 cm 1 , 
Eg(0, 2) = Eg(0, 0) + 6B! — 36D! — 6a’, = 92798.9 cm! , 
Eg(1, 0) = Eg(0, 0) + w! — 2w!x!, = 93721.5cm ^! , 
Eg(1, 1) = Eg(0, 1) + w, — 2w!x!, — a! = 93762.0 cm™}, 
Eg(1, 2) = Eg(0, 2) + w! — 2w!x!, — 3o], = 93842.7 cm ! . 





The conditions specified only give six transitions. They are: 


(0,0) R(0) at w = Eg(0, 1) — Ex(0, 0) = 90242.3 cm! , 
(0, 0) P(1) at w = Eg(0, 0) — Ex(0, 1) = 90085.0 cm! , 
(0,0) R(1) at w = Eg(0, 2) — Ex(0, 1) = 90200.9 cm! , 
(1,0) R(0) at w = Eg(1, 1) — Ex(0, 0) = 91156.2 cm ^! , 
(1,0) P(1) at w = Ep(1, 0) — Ex(0, 1) = 91401.3 cm! , 
(1,0) R(1) at w = Eg(1, 2) — Ex(0, 1) = 91510.1 cm! . 


10.3.5 Astronomical spectra 


Most astronomical electronic spectra of molecules are observed in absorp- 
tion. Recent observations by the Far-Ultraviolet Space Explorer (FUSE) 
satellite have extensively monitored the Werner and Lyman bands of H32 
in absorption in the ISM. Such observations provide a direct observational 
handle on cold molecular hydrogen in the ISM. These bands are observed 
in absorption against starlight. See Fig. 10.13 for example. 
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Figure 10.14 shows part of the Lyman band of Hp. This transition is 
BIXi-X Ds and is electric dipole-allowed. Absorption for the vibrational 
ground state starts at about 1100 À and runs to shorter wavelengths. The 
abundance of H» is so great in interstellar clouds that absorption from the 
lowest rotational level, J — 0, is optically thick. Despite the large B con- 
stant of Hz, and the relatively small Einstein A coefficient for the transition, 
absorption from several rotationally excited levels can be seen. The spec- 
trum also shows a line due to HD, which is shifted due to mass effects that 
influence the vibration-rotation energies in both electronic states. Deu- 
terium fractionation effects at cold temperatures mean that although a 
ratio of HD to H» can be obtained from this spectrum, this does not give 
directly the ratio of D to H abundances. 

Electronic spectra can be observed in emission from both comets and 
planetary aurorae. In both these cases the emissions are the result of non- 
thermal excitations. For example, Jupiter's aurorae emit strongly in both 
the Lyman and Werner bands as a result of electrons travelling down the 
magnetic field lines in the polar regions and exciting the Hz gas which 
makes up most of Jupiter's atmosphere. 

Except for Hz, whose rotational fine structure transitions are widely 
spaced, the rotational structure of the electronic transitions are often not 
fully resolved, leading to band structures. However the bands have char- 
acteristic profiles, for example containing band heads, which allow them 
to be identified as molecular in origin. 

Molecules in giant molecular clouds and other cold interstellar envi- 
ronments do not possess sufficient energy for their electronic spectra to be 
observed in emission. However, molecular absorptions can be seen against 
suitable light from a more distant star. As these molecules exist at low tem- 
peratures, only the ground vibrational state and a few rotational levels are 
occupied. For light molecules only one rotational state may be occupied. 
These species give rise to fairly simple spectra with transitions to several 
vibrational states of the upper electronic state, each with a few, maybe 
even only one, rotational fine structure line. 

The first three molecules observed in the ISM were CH, CH* and 
CN. These somewhat unusual species were all detected by their electronic 
spectra which lie in the visible. Figure 10.15 compares term diagrams, 
which can be thought of as the molecular equivalent of Grotrian diagrams, 
for these three species. CH and CH'* are light species so only transitions 
from their rotational ground states are usually observed. The CN ‘violet’ 
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Fig. 10.15. Term diagrams giving the rotational structure only of the strongest 
electronic interstellar bands of the CN, CH and CH* molecules. Note the 
more complicated structure of CH is caused by A-doubling in both the ground 
and excited electronic states. [Reproduced from P. Thaddeus, Ann. Rev. Astron. 
Astrophys. 10, 305 (1972).] 


band has a larger transition dipole and CN is heavier and so has smaller 
rotational constants. Two or three rotational levels may be occupied, which 
leads to either three or five rotational fine structure transitions for each 
vibrational band. Figure 10.16 gives an example where only the J” = 0 
and 1 levels of CN are occupied. The population of CN rotational levels 
are sensitive to the local temperature of the environment and can be used 
to measure the local microwave background, as discussed in Sec. 9.6.1. 
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Fig. 10.16. Optical spectrum recorded towards ¢ Oph showing rotationally- 
resolved absorption due to the CN ‘violet’ band. Absorption by the R(2) line is not 
observed. Adapted from P. Thaddeus [Ann. Rev. Astron. Astrophys. 10, 305 (1972).] 


10.4 Non-!x Electronic States 


Molecules which do not have !X (i.e. S = 0, A = 0) electronic states have 
more complicated spectra. In this case it is necessary to consider the cou- 
pling of the rotational motion to the spin and/or orbital angular momenta. 

This can get fairly complicated as, in particular, there is more than one 
way to do this. The choice of coupling scheme depends on the molecu- 
lar species and even the electronic state in question. These different cou- 
pling schemes are known as Hund's cases. Here only one astronomically, 
and atmospherically important species will be considered: the hydroxyl 
molecule, OH. 

The ground state of OH is X?II which means that S = 4 and A = 1. This 
system is treated using Hund’s case (a) which starts by coupling the spin 
projection along the molecular axis, X, to the electron orbital projection on 
this axis. The resulting projection of the total electron angular momentum, 
Q, can take two values as this is not vector addition but the addition of two 
projections which either point in the same direction or opposite directions. 
The total projection can thus take the values 


Q — |^ — EL, A -X. (10.19) 


For the case of ?II electronic state, therefore, Q equals i and 3. 

The X °H electronic state of OH is split into two states which are des- 
ignated "Tl; and “Tl, , or indeed sometimes ?II, and *II_. The separation 
between these two states is due to relativistic interactions and therefore 
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completely equivalent to the fine structure found in atomic levels. For OH 
the ?IT; is the lower of the two states. 

The extra complication for molecular spectra arises from the addi- 
tional coupling to the rotational motion which is strictly represented by 
the quantum number N. Note that N equals J in the usual case when 
S = A = 0. However for states with A 7 0, the electronic angular momen- 
tum couples with N to give J. This coupling to the rotational motion splits 
each of the energy levels into a ‘+’ or ‘—’ component. The result is known 
as A-doubling. A-doubling effects are particularly strong for molecules such 
as OH which have A — 1. The structure of the low-lying energy levels of 
OH, and the resulting transitions between them, is given in Fig. 10.17. 
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Fig. 10.17. Low-lying energy levels of the OH molecule showing the wavelengths 
of allowed radio and near-infrared transition. The expanded portion shows the 
hyperfine splitting of the lowest levels and the frequencies of the associated hyper- 
fine transitions. 
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Transitions between the A-doublets are allowed. The lowest of these 
lies at radio frequencies with a wavelength of 18cm. Observation of this 
transition in 1963 led to OH becoming the fourth molecule observed in the 
interstellar medium and the second species observed at radio frequencies. 

To add one more layer of angular momentum coupling, the A-doublets 
are themselves split by hyperfine coupling, which couples the rotational 
motion quantum number J to the total nuclear spin quantum number I. In 
the case of OH, !6O has a zero nuclear spin and the hyperfine effects arise 
from coupling J to the i = 5 of the H nucleus. Hyperfine effects cause the 
18-cm line to split into four components (see the inset in Fig. 10.17). These 
lines can be observed in high-resolution studies. 


10.5 Maser Emissions 


Microwave amplification stimulated emission of radiation (maser) action 
is observed from at least 36 molecules including SiO, OH and water, usu- 
ally at infrared or microwave frequencies. The population inversion nec- 
essary to cause maser action can be created by a number of mechanisms 
including optical pumping, radiation trapping in certain long-lived levels 
and selective collisional excitation of the masing molecule. 

Figure 10.18 depicts a simplified case of a maser driven by collisional 
excitation. In this figure the levels are vibration-rotation levels of the elec- 
tronic ground state. If level A is excited in some collision process, such as 


A 
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excitation 


maser action 


B 





Rapid 
emission 


Ground state 


Fig. 10.18. Typical scheme for a three-level masing system driven by selective 
collisional excitation. 
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scattering with H3, and level B is not, then the population of level A can be 
greater than the population of level B. This situation gives a non-thermal 
population and can lead to maser action which, by its nature, is highly 
directional. 

The detailed mechanism which causes masing in each system can be 
regarded as an accident of the physics for that species. Molecules such 
as SiO have been observed masing in over 400 stars. OH and water are 
Observed to mase in the ISM. Such maser action is often very time variable 
and — since masing is very dependent on the detailed physics — carries a 
lot of information about the local environment. Indeed the original obser- 
vations of what turned out to be OH masers were so puzzling that they 
were dubbed ‘mysterium’. Instead of being a possible new element, it was 
considered possible that these emission lines with strong, wildly fluctuat- 
ing intensities were due to an alien lifeform. However it was demonstrated 
that simple physical phenomena giving maser action could give rise to 
such effects. The masing hyperfine lines of OH are now widely studied. 
Figure 10.19 gives an example of maser emission from a variable late-type 
supergiant star. Water and SiO maser lines can both be observed from this 
object. The velocity structure of the masing line displayed in Fig. 10.19 
indicates that the emissions are coming from an expanding shell of gas at 
constant velocity. 
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Fig. 10.19. Velocity profile of OH 1612 MHz maser emission from the variable, 
late-type supergiant star VX Sagitarii, observed using the Lovell 76m antenna. 
(A.M.S. Richards, private communication.) 
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Problems 


10.1 


10.2 


10.3 


10.7 


Interstellar CO is observed absorbing infrared radiation at 2134.63, 
2138.46, 2146.16, 2150.03 and 2153.94cm '!. Identify the transitions 
observed inthis spectrum and explain what astrophysical information 
might be obtained from its study. Where might one expect to observe 
an infrared spectrum of CO in (a) absorption and (b) emission? 

One wishes to observe the molecules ?C!6O and C160. How do 
(a) the J = 1-0 rotational transition frequency and (b) the v = 1-0 
fundamental vibrational frequency depend on the reduced mass of 
the system? For ?C!60, the J = 1-0 transition lies at 3.86cm ! and 
the fundamental frequency lies at 2170 cm !. Assuming integer val- 
ues for the atomic masses, estimate these quantities for C160. Why 
might it be necessary to observe different transitions for 1?C!6O than 
13C160 to determine their relative abundances? 

Observations are planned of CO pure rotational emission spectra 
using the J = 5-4 and J = 20 -19 transitions. Assuming CO is a rigid 
rotor with rotational constant B = 1.93cm7!, estimate the transition 
frequency of the two transitions. Explain which of these estimates 
will be more accurate. 

The ?CH* molecular ion has a rotational constant B = 11.94cm ^! 
and a fundamental vibrational frequency of w = 2075.5 cm 1. A cold, 
interstellar column of ^CH* in its J = 1 rotational state is observed 
in (a) the far-infrared and (b) the mid-infrared. Give the quantum 
numbers and wavenumber (in cm !) of the transitions you would 
expect to observe. 

Assign quantum numbers transitions labelled CO in Fig. 7.6. You can 
assume these are all pure rotational transitions. 

Figure 10.9 shows a vibrational spectrum of CO. Estimate the wave- 
length of the band origin of CO from this spectrum. Starting at the 
band origin, assign rotational quantum numbers as far as the maxima 
in the two branches of the spectra. Hence obtain an estimate of the tem- 
perature of the CO gas involved in the absorption. You may assume 
that CO is a rigid rotor with a rotational constant B = 1.93cm"!. 
Estimate the wavenumber in cm~! at which cold interstellar H, will 
absorb in the Werner band. Assume H; is a rigid rotor and a harmonic 
oscillator, and use the constants B,, we and T, given in Table 10.2. Only 
consider absorption to the three lowest vibrational states, v' — 0,1,2, 
in the excited, C !II,, electronic state. 


SOLUTIONS TO MODEL PROBLEMS 


Several solutions below give precise numerical answers to problems for 
which an estimate was requested. For these problems, indicated * below, 
correct working may yield answers differing slightly from those given. 


Chapter 1 

11v-— ch. 
Line 1: v = 955: 797 x 2.998 x 105 = 100.1 kms ^; 
Line 2:0 = 830-4 x 2.998 x 108 = 100.1 kms; 
Line 3: v = 5772772619 x 2.998 x 105 = 99.9 km s71; 
Line 4: v = $8:57—0625 x 2.998 x 108 = 100.0 kms !; 


All lines show a Doppler shift of 100 kms ! which is the speed the star 
is moving away from earth as the shift is positive. 





Chapter 2 
2.1* Doppler shift, measuring v, as approximately —600 km s“!, is 
v, 6 x 10° 
Av = T TTC PA 99.023 — 0.198 GHz. 


Hence v = 99.023 + 0.198 = 99.221 GHz. The laboratory frequency for 
H509 is 99.223 GHz. 


Chapter 3 


3.1 E, = — t. Lya is n = 2-1 so 


E, —E,= ns (a — i) = 109677.58 x i = 82258.19 cm™!. 
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Rə» is defined for a one-electron atom with an infinite nuclear mass. 
Iron is a heavy atom, so the reduced mass for an electron plus an iron 
nucleus is very close to the mass of an electron, me, which is the mass 
used in defining Ræ. Lya for Fe** therefore lies close to 26? x 0.75 x 
109737.31 = 5.564 x 107 cm 1. Such transitions will only be observed 
from very hot environments such as intergalactic plasma. 

Selection rules: Al = +1, An any. Observed emissions: 


4-3 Pa at 18751 A; 
4-2 Hf at 4861 A; 
4-1 Ly* at 972 À; 

3-2 Ha at 6563 A; 
2-1 Lya at 1215 À. 


Starting from the 4s level, 3-1 Ly? at 1026 Å will be observed, but 
not Lyy. 
The cascade of emission from the H 4p level can be depicted as: 


3. 


N 





S p d f 
n=4 —— —— =— => 
n=3 
n-2 — 
n=1 


3.3 AX 1= Efa = +) for Ha nı = 2, m = 3. Only the Rydberg constant R 
2 


ny 
changes between H and D, via the reduced mass, so: 


Ap Ru jn Myme Mp + me 


eee S — = ———— x ——__, 
AH Rp HD My + me Mpme 
Ap 18361 36714 


2D — T X LL = 0.999728. 
Au 1837.1 536704 — 299728 
Therefore Ap(H/9) = 6561.00 A. 
À AH Ra 1 


Rm SS SE SESS SS EE 
AX Ay—AD Rọ -Rọ 1-Rgu/Rp 
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so the resolving power R is independent of nı or np. 


1 
| 1— 0.999728 


Note that telescope resolving powers are usually given to one digit 
and a power of ten. It is meaningless to specify the resolving power 
to significantly more digits than this. 

Since R is independent of nı and nz, the same resolving power is 
required for Lyo. The main problem with observing both H and D 
transitions together is that if there is approximately 10? more H than 
D, then if the D transition is observable, the H one is likely to be 
optically thick. 


R ~ 3700. 


From Sec. 3.8.1, the volume of an H atom with principal quantum 
number n is 


4 
V,a 37" x 1.48 x 10?! mè. 


(a) Star with N = 1016 cm ? = 10? m ?, so each atom can occupy up 
to about 10 7? m°. Equating gives: 


4 
37" x 148 x 10?! ~ 10772 


which gives n ~ 23. 

H II region N = 10* cm ? = 101? m ?,so each atom can occupy up 
to about 101? mê. Equating with V,, gives n ~ 2300. 

The partition function sum is finite since in any environment the 
H atom is confined in a finite volume. This truncates the sum at 
some maximum value of n, e.g. 23 for the star and 2300 for HII 
region considered above. 


S 


ADD R(4 — i) For 80a, nı = 80 and n} = 81 which gives 
T 2 

A = 2.37668 using Rx and à = 2.37797 using Ry. The resolving 

power 

A 2.377 

AA 2.37797 — 2.37668 


A one-digit estimate of the resolving power gives R — 2000. 


R= z~ 1800. 


Chapter 4 


4.1 


Be has 4 electrons; its ground state configuration is 1s*2s*. Both orbitals 
are fully occupied, giving a term 1S and a level 'So. 


4.2 


4.3 


4.4 


4.5 
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The excited state has configuration 1s?2s!3d! giving l(2s) = 0, s(2s) = 
and (3d) = 2, s(3d) = 1 5. Adding these gives L = 2 and S = 0 or 1. Thes 
correspond to terms 1D, 3D, and levels ! D», ?D4, ?D5, ?Ds. 

4d? means 1; = 2, sı = 1 and I; = 2, s2 = $. Adding these L = lı + h = 
0,1,2,3,4 and S = s, + s2 = 0,1. Without the Pauli Principle the terms 
1S, ! P, ID, !E, !G,?sS,?P, 3D, ?F and °G would be possible, but the Pauli 
Principle imposes the rule that for indistinguishable electrons the sum 
of L + S must be even. Applying this rule gives 'S, ?P, ! D, °F and !G 
only. If Hund's rules are obeyed for this configuration then the energy 
order would be ?F, ĉP, !G, ! D, 'S, going from lowest to highest. 


1 
2 
e 


Tem L S J=]|L-S|,... L+S 
1S 0 0 0 

oP 1 1 0,1,2 

iD 2 @ 2 

3F 3 1 2,3,4 

1G 4 0 4 


giving 35, 3ES, 3F4, 3Po, 3P}, 3P5, IG; ID», 1S. 3B, should be lowest in 
energy. All terms are even parity since (—1)^*^ = (—1)?? = 1. 
^F? means S = $, L = 3, J = Z. It will have 2J + 1 = 8 states (which 


L 35722 
S = 5/3/33 





are Mr sublevels). The other levels are given by J = L4 
which gives S, "FS and "FS in addition to FS. 

As S = 3, the system must contain at least 3 electrons. A possible three- 
electron configuration could be 1s2p3d. Other configurations are also 
possible, but not 2p? which cannot give ^F? by the Pauli Principle. 

1D has L = 2 and S = 0 so it can exist only for J = 2 i.e. 'D2, not !Dy; 
°D; would mean 2S + 1 = 0, ie. S = —4 but S must be > 0; 


?P means S = L = 1 giving integer J values of 0, 1 or 2 but not 5. Half- 

integer J’s mean odd numbers of electrons and 25 + 1 even. 

(a) 1s?2p: one open shell electron which has | = 1, s = 3 and term ?P°. 

(b) 1s2s3s:3 open shell electrons, all with / = 0 so L = 0. All have s = L 
adding them sequentially S12 = 00r1,s05 = S12 - s = L L 3 (note 
the two separate values of S = 4). The terms are 7S, 7S and 4S. 

(c) 1s2p3p also has 3 open ipe electrons. The spin coupling is the 
same as (b), giving S = i, L 3. Coupling lı = 0, h = 1 and l; = 1 
gives L = 0, 1 and 2 since lı does not contribute. So the resulting 
terms are 25, 2S, ^S, °P ?p ^P 2D, 2D and ^D. 
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4.6 Only the 3d5p electrons need to be considered as the others are all 
in closed shells. These have l4 = 1, s, = i and h = 2, s2 = i giving 
L = 1,2,3 and S = 0,1. The terms are therefore !P, ! D, | E, 3P ?D and °F. 
The lowest is ?F which has levels ?F;, ?F4 and Fy. 


Chapter 5 


5.1 Ha is n = 2-3. Hell has Z = 2 which means the transition will be 
scaled by a factor Z? compared to H I. This means that He ll n = 4-6 
lies close to Ho and that this transition satisfies: 


where Rye = ( Mne Ra Assuming Mye = 4Mg gives Rue = 


Mute + Me 

109722.4cm~! and a wavenumber of 15239.22cm~! or \ = 6562.0 A. 
5.2 1640 À corresponds to a Hell Balmer a transition with n = 3-2. If 

n = 3-2 emissions are observed then Lyman f, n = 3-1 at 256 À and 

Lyman o, n = 2-1 at 304 Å. 
5.3 Shell A: H rand Hell recombination spectra. 

Shell B: HI and He 1I recombination spectra. 

Shell C: HI recombination spectra only. 

Shell D: no recombination spectra will be observed as there are no ions. 
5.4 (a) 1s2s has]; =], = L = 0 and sı = s? = L giving 1S, and 35}. 

(b) 1s2p has |; = 0, h = 1, so L = 1 and sı = s2 = 3,80 S = 0 or 1. 

This gives TP, and 3Po, ?P, and 3P>. 

Decays in approximate order of strength (strongest first): 


1s2p - 1s? 'P?—'So, resonance line, 
1s2p- 1s? apo —16,, intercombination line, 
1s2p - 1s? ?P? -!S,, magnetic quadrupole line, 
1s2s - 1s? 3S; - 15, spin-forbidden magnetic transition. 
The two J = 0 levels cannot decay to the ground state since transitions 
with J = 0—0 do not occur for any selection rules. 
5.5 (a)-(b) 1s?2p - 1s2s3s is a forbidden transition: parity changes 
(odd to even) but two electrons jump. 
(a)-(c) 1s?2p - 1s2p3p is allowed. 
(b)-(c) 1s2p3p - 1s2s3s is rigorously forbidden by the Laporte rule. 


Chapter 6 
6.1 (a) Na 3p is °P, "P$; 4d is D;, ?D;; 
2 2 


Allowed transitions are ?P9 EAD ?po -?D; ,?P9 -?Ds. 
2 2 2 


6.3 


6.4 
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(b) Na 3d is *Ds, *Ds; 5f is T, FS; 
Allowed transitions are ?D; =F, "D; E C, ?D; -EFE. 
(c) Na 4s is "84; 4d is ?D;, *Ds; 
Both are even parity so no electric dipole transitions possible. 
Electric quadrupole transitions have AJ = 0, + 1, + 2, parity con- 
served. Both *S; -*D3 and °S, -*D; are quadrupole-allowed. 
(d) K 4s is ^81; 4p is "PS, "PS; 
Allowed transitions are 2S i= 2Po 25 1 -Tp. 
All the allowed transitions should be observable in absorption in stel- 
lar atmospheres of suitable temperature. Only (d) K 4s - Ap is a possible 
interstellar absorption feature as all others are between excited states. 
E(4s) = —35010cm !; E(4p)= —35010 + 12985 = —22025cm 5; 
E(5p) = —35010 + 24701 = —10309 cm !; E(6p) = —35010 + 28999 = 
—6011 cm !. 
The quantum defect, un is given by: 











(n ior Ent)? : 
so (4p) = 1.768; u(5p) = 1.737; u(6p) = 1.727. 
Assuming (7p) = 1.72 gives E(7p) = —3936cm' ! and 4s-7p at 
31074 cm ''. These transitions correspond to the 4 °S, -n PS series; 


for this series the 4s -7p transition is observed at 31070 cm™!. 


E(nl) = 





Need the splitting AE = 13080.5 — 12988.7 = 91.8 cm™!. 
Level L S J i[JJ*1-L(L-1)- S(S4-1)] 
ap 1d 3 + 
2 
Ap? "E E —1 


So 3A’ —91.8cm !; A = 61.2cm7!. 

4d has L = 2, S = $ so J = 3 or 3; 3p has L=1,S5=4s0J=3 
or 3. The selection rule is AJ = 0, +1 so the allowed transitions are 
*Ds — "PS, ?D; — "Pe and ?D; => "Pj (see Fig. 6.5). 

If 4d emits to 3p, then the following cascade emissions must also fol- 
low: 4p —4s, 4s —3p, 3p -3s and 4p -3s. 

The triplet emissions should have different intensities according 
to their statistical weight. If all three transitions have the same 
intensity, they are optically thick and column density of Na cannot 
be determined from this observation. 
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where Zeff = 4. The quantum defect gives the departure of the energy 
levels of the outer electron from the pure hydrogenic case. It is largely 
determined by penetration effect, with the values generally decreas- 
ing with l; for a given / it usually only depends weakly on n. 

Using E = 520178cm ! and n = 2 gives (2s) = 0.163; E = 
217329 cm ! and n = 3 gives (3s) = 0.158. So using u(4s) = 0.155 
and n = 4 gives E(4s) = 118763 cm '!. Observed value is 118830 cm '!. 
(a) Hund's rules give *S? lowest, then ?D? and ?P? highest. 

(b) Levels are ‘SS, "D$, "D$, FS and T. 








Level L S J àÀiJ-c-10-L(.-c1)-S(S41) 
"D$ 2.1 3 —1.5 

"D? 24 3 +1 

"T; 38 3 2 

^F? 3 1£ 41.5 





?F9 means S = 1, L = 3, J = 2 and odd parity. ?F? term also gives 

?F$ and ?F?. The term could arise from the configuration 1s?2s?2p3d 

(other answers are possible). 

In order of diminishing strength: 

?F? — 3D; is dipole-allowed and will be strong. 

3F9—- 1D; has AS Æ 0 and is an intercombination line. 

?F2 -?P, has AL = 2 and is a forbidden transition. 

?F?2-?P$ does not change parity so is completely dipole- 
forbidden by the Laporte rule. 
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Chapter 7 


7.1 (a) 4s? gives 1S); 4s4p gives MSS dus apo and Spo: Decays in approxi- 
mate order of strength (strongest first): 


4s4p — 4s? 1P? -!S,, resonance line, 

4s4p — 4s? ?P? —1So, intercombination line, 

4s4p — 4s? ?p9 - 155, magnetic quadrupole line. 

?P8 does not have an allowed decay to the 4s? 'Sp ground state. 


(b) 4s4f gives 1F9, SES, ?pe and RES 4s5d gives 1D», ?D4, ?D» and 3D3. 
Allowed dipole transitions: 


1F9 - 1D, ?F2 -?Ds, ?Fg -?D», Fg -?Ds, °F -5D», ?F9 -?Ds, 
3po mM 3D 
4 3: 


All the allowed transitions could be observed in emission as part of the 
Ca I recombination spectra in HII regions. The transitions (a) could be 
seen in absorption from the ISM. 

7.2 N* is C-like: 1s?2s?2p?. Gives terms °P, ! D and !S as C. Hund's rule (1) 
means ?P will be the ground state. 
Transitions 1s?2s?2p?(?P)5f — 4d. Need to consider coupling between 
?P core and the extra electron. 


Upper state: ?P gives L = 1, S = 1; 5f gives] = 3, s = 4}. So L = 2,3,4; 
S = 5,3 giving ?D^, *D°, ?F?, 4F°, ?G? and 4G°. 

Lower state: ?P gives L = 1, 5 = 1; 4d gives l = 2, s = $. So L = 1,2,3; 
S = 5,3 giving P, ^P 7D, *D, °F and *F. 

Allowed transitions, 12 in total: 


2pe Iy 2P 2pe = 2D 2pe cas 2F 2pe — ?D 2pe > 2B 260 = 2B 
4De =y 4p, 4De = 1D, 4De M 4F, 4po y 4D, 4po ET 4F, 4G9 E 4k 


(In fact inclusion of fine structure effects leads to 71 distinct transi- 

tions.) The strongest transition is *G? — ^F as: 

(a) *G? has the highest statistical weight and therefore occupancy; 

(b) *G? only emits to ^F. Terms D? and F° have competing decay 
routes. 


7.3 The transitions lead to the following energies relative to the 2p? °Po 
ground level: 
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Level Energy (eV) g Population, T = 10000K 
2p? 3Po 0.000 1 1.000 

2p? ?P, 0.014 3 2.958 

2p? ?P; 0.038 B 4.814 

2p? !D; 2.514 5 0.405 

2p? 'So 5.356 1 0.0047 

2p3p?P; 40.849 5 9.1 x 10718 


Note that several of these levels can be calculated in different ways 
which provides a consistency and accuracy check. The level 75+1L; has 
statistical weight g = 2J + 1; individual degeneracies are given in the 
table. Population relative to the ground state is given by Boltzmann's 
distribution g exp(— 7), where it has been assumed that for the ground 
state E = 0 and g = 1. Results for T = 10000K are given in the 
table. At a typical nebular temperature of 10000K the fine struc- 
ture levels of ?P are populated almost in proportion to their sta- 
tistical weights, so the populations of these levels are not a good 
indicator of temperature. The ! D» and !S; levels have reasonable ther- 
mal populations, which are temperature sensitive, so these levels will 
give reasonably strong transitions which can be used to determine 
temperature. 

The thermal population of 2p3p 3P, is negligible; this level can only be 
populated by some other mechanism. 


Chapter 8 
8.1 


Ovi Lya is at 4032846cm !, which corresponds to 500eV. This 
lies in the X-ray and such transitions can only be observed using 
satellites. 

8.2 1s?2s?2p3s?3p63d ^ 4s?4pSAd' 4f ^ 5s?5p65g?6s?. Hg * 5d!^6s gives term 
?S and level ?S;. Both this configuration and the 5d?6s? configuration 
are even parity so dipole transitions between them are forbidden by the 
Laporte rule. The 5d*6s? °D; state is metastable since it cannot decay 
by a dipole transition. 
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Chapter 9 
9.1 Boltzmann’s distribution gives the probability, Pj, of being in level J 


P BJ +1) 


For a maximum: 


dB 
ri = 2exp|- 2002 | — QJ +1) = exp - ar) 2 zm: 


a kT\? 1 
e= =h) o 


9.2 Using the formula derived in answer 9.1: 
(a) J = 1; (b) J = 7; (c) J = 22 or 23. 
Using Eq. (9.34): 
(a) Effectively zero; E 0.00396; (c) 3596. 
9.3 Zero point energy = ^? = 1037.8cm '!. Within the harmonic approxi- 


$ Z 
mation w = (£) , SO: 


1 1 
z H12 z V 12 13 + 16 z E —1 
(713 = W12 (=) = 2075.5( eu 143 — 2069.4cm 


Zero point energy is = ^? = 1014.7 cm !. This assumes that k is the 
same for both "CH* and PCH'*. 

13CH* may be optically thin when ^CH' is optically thick. 

At low temperature, fractionation effects will increase the ratio 
n(PCH") to n("CH*) compared to n(?C) to n(?C) since “CHT has 
the lower zero point energy. 


Chapter 10 


10.1 The gaps between the lines are 3.83, 7.70, 3.87 and 3.91 cm 1, which 
correspond (approximately) to 2B, 4B, 2B, 2B. The 4B gap is the band 
centre and the transitions are therefore P(2), P(1), R(0), R(1) and R(2), 
respectively. CO will absorb in the infrared both in the ISM and in 
the atmosphere of cool stars. CO infrared emissions are less common 
but can occur in the warm ISM (e.g. shocked regions or nebulae). 

10.2 (a) Assuming a rigid rotor, the J = 1-0 transition at 2B and 

Bap 
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10.5* 
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(b) Assuming a harmonic oscillator, the v = 1—0 transition is at hw 

and w x p 3. 

. MmnMo , nm MisMo . 

H12 Mp 4- Mo Mo’ 13 Mis + Mo X Mo' 

u2 _ Mi(Mis + Mo) 

LM = 0.9560, 

H3 — Mi3(Ma2 + Mo) 
assuming My) = 12, Mj3 = 13 and Mo = 16u. J = 1-0 for 8C!°O 
will be at 3.86 x 0.9560 = 3.69 cm. 
hw for ’C!O will be at 2170 x 0.9560 = 2122cm-!. 
The abundance of ?C is much lower than that of C so observable 
lines for C!^O may well be optically thick for C160. 
Rigid rotor transitions at 2BJ' so J = 5-4 will be at 10B = 19.3cm ! 
and J = 20-19 will be at 40B = 77.2cm7!. The estimate for J = 5-4 
will be more accurate since the contribution from centrifugal distor- 
tion will be lower. 
(a) Far-infrared, rotational spectrum: 

In emission: J = 1-0 at 2B = 23.88cm7!; 

In absorption: J = 2-1 at 4B = 47.76cm"!. 

(b) Mid-infrared, vibration-rotation transitions, observed only in 
absorption: 

P(1) at w — 2B = 2051.6cm-!; R(1) at w + 4B = 2123.3 cm™t. 
Even for J ~ 30, CO behaves as a near rigid rotor with transitions 
spaced by 2B, ~ 3.8cm ''. Measuring approximate wavelengths 
from Fig. 7.6 gives: 





A(um) hwu(cm?) J'—J" 


186 53.8 14-13 
174 57.5 15-14 
162 61.7 16-15 
153 65.4 17-16 
145 68.9 18-17 
138 72.5 19-18 
131 76.3 20-19 
123 81.3 21-20 
118 84.7 22-21 
113 88.5 23-22 


108 92.6 24 —23 
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10.6* Figure 10.9: central feature (where a line is missing) at Ag ~ 4.66. The 
R-branch lies at A < Ao with lines up to the R(0), R(1), R(2), R(3), R(4), 
R(5) up to the maximum at R(6). The P-branch lies at A > Ao with 
lines P(1), P(2), P(3), P(4) up to the maximum at P(5). R(6) means 
J'-J" =7-6 ie. J” = 6; P(5) means J' - J” = 4-5 so J” = 5. So the 
highest occupied level in the absorbing gas is J" — 5 or 6. Using the 
formula derived in answer 9.1 gives T ~ 168K for Jmax = 5 or 235K 
for [max = 6. 
Note: this is very approximate. A reliable temperature estimate can 
be obtained by a fit to the intensity of all the transitions in the band. 
10.7 Denote the energy levels Ex(v", J”) and Eg(v', J’) for the ground and 
excited electronic states only. As the H» is cold, only J” = 0 or 1 will 
be occupied. Then, using constants for the appropriate electronic 
state taken from Table 10.2, 


wi 


Ex(0,0) = = 2200.6cm ^1 , 

Ex(0,1) — EQ, 9 + 2B” = 2322.3cm™}, 
Ec(0,0) = T, + e = 101311.7cm ! , 
Ec(0,1) = Ec(0,0) + 2B] = 101374.4cm ! , 
Ec(0,2) = Ec(0,0) + 6B! = 101499.9cm ! , 
Ec(1,0) — s 0) + w! = 103755.5cm !, 
Ec(1,1) = Ec(0, 1) + w, = 103818.2cm  , 
Ec(1,2) = Ec(0,2) + w, = 103943.6cm ! , 
Ec(2,0) = Ec(1,0) + w, = 106199.2cm  , 
Ec(2,1) = Ec(1,1) + w, = 106261.9cm 1, 
Ec(2,2) = Ec(1,2) + w, = 106387.4cm^! 





The conditions specified give twelve transitions. These occur at: 
(0,0) Q(1) at o = Ec(0,1) — Ex(0,1) = 99052.1cm ^! , 


(0,0) R(0) at w = Ec(0,1) — Ex(0,0) = 99173.8cm ! , 
(0,0) P(1) at w = Ec(0,0) — Ex(0,1) = 98989.4 cm! , 
(0,0) R(1) at w = Ec(0,2) — Ex(0,1) = 99177.6cm^! , 


) ) ) 
) ) ) 
) ) ) 
(1,0) Q(1) at w = Ec(1,1) — Ex(0,1) = 101495.9 cm 
(1,0) R(0) at w = Ec(1,1) — Ex(0,0) = 101617.6 cm7 
(1,0) P(1) at w = Ec(1,0) — Ex(0,1) = 101433.2cm ! , 
) ) ) 
) ) ) 
) ) ) 
) ) ) 
) ) ) 


, 


$ 


(2,0) Q(1) at w = Ec(2,1) — Ex(0,1) = 103939.6 cm 
— Ex(0,0) = 104061.3 cm7 
— Ex(0,1) = 103876.9 cm7 
— Ex(0,1) = 104065.1 cm7 


(2,0) R(0) at w = Ec(2,1 
(2,0) P(1) at w = Ec(2,0 
(2,0) R(1) at w = Ec(2,2 


, 


, 


1 
1 
1 
(1,0) R(1) at w = Ec(1,2) — Ex(0,1) = 101621.3cm ! , 
1 
1 
1 
1 
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